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Abstract 



High energy reactions necessarily involve large momentum transfer to the target's con- 
stituents, which thus rebound near light speed. The experimentally observed correlation 
between initial and final states measures the system's response along the light cone. A con- 
venient theoretical description of such reactions is light-front dynamics where the properties 
of physical states are described along the advance of a wavefront of light. 

In this work we report on using light-front dynamics to describe generalized parton 
distributions, which are correlation functions encountered in virtual Compton scattering 
at large momentum transfer. This two photon process requires pair annihilation contri- 
butions which are subtle to handle in light-front dynamics. We are careful to derive such 
contributions in the light-front Bethe-Salpeter formalism. This derivation highlights the 
connection to the Fock space picture as well as provides a toy model for generalized parton 
distributions. 

Ultimately we use light-front dynamics to build a phenomenological parametrization 
for the proton's generalized parton distributions. This task is non-trivial due to the field 
theoretic constraints required of these distributions. To meet the constraints imposed by 
Lorentz covariance, we approach generalized parton distributions from the perspective of 
double distributions. Despite avid phenomenological use, surprisingly little work has been 
done to calculate double distributions. We first focus on a number of simple pedagogical 
examples and demonstrate how double distributions can be correctly determined. The 
double distributions we derive, including those for the proton, satisfy known constraints 
and are adequate for phenomenological estimates of cross sections. The models used allow 
one to study the interplay between the light-front formulation and Lorentz covariance. 
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Glossary 



BSE: Bethe-Salpeter equation. 

Covariant: Lorentz covariant, i.e., transforming properly under the Lorentz group. 
DD: Double distribution. 
DIS: Deep-inelastic scattering. 

DSE: Dyson-Schwinger equation, a coupled field theory equation of motion for a Green's 
function. 

DVCS: Deeply virtual Compton scattering. 

GDA: Generalized distribution amplitude, or two pion distribution amplitude. 
GPD: Generalized parton distribution. 
IMF: Infinite momentum frame. 

LFTOPT: Light-front time-ordered perturbation theory. 

Light-cone coordinates: The coordinates that describe the world sheets of light, namely 
in our conventions x"^ = {x^ it x^)/^/2. 

Light-front energy: The Fourier conjugate to light-front time. 
Light-front time: The hypersurface defined by x'^ = constant. 
PDF: Parton distribution function, or quark distribution. 

Polynomiality: Property of the moments of generalized parton distributions required by 
Lorentz invariance. 

Positivity: Bounds that are required of generalized parton distributions due to the posi- 
tivity of the norm on Hilbert space, properly referred to as positivity bounds. 

QCD: Quantum chromodynamics, the quantum field theory of colored quarks and gluons. 

QED: Quantum electrodynamics, the quantum field theory of electrically charged particles 
and photons. 
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Chapter 1 

Introduction 



The description of the bound states of strongly interacting particles calls for a synthesis of 
relativity and quantum mechanics. While this synthesis is afforded by quantum field theory, 
the direct investigation of bound states from a covariant framework has been limited. Other 
complimentary methods are desirable for both varied intuition and differing calculational 
strategies. 

In a Hamiltonian approach, one attempts to calculate wavefunctions for the physical 
states of the theory. More than a half century ago, Dirac's paper on the forms of relativistic 
dynamics [1] introduced the front-form Hamiltonian approach, which differs from conven- 
tional Hamiltonian dynamics by a relativistic boost to infinite momentum. Applications 
of this form of dynamics to quantum mechanics and field theory were overlooked at the 
time due to the appearance of covariant perturbation theory. The reemergence of front- 
form dynamics was largely motivated by simplicity as well as physicality. The light-front 
approach has the largest kinematic subgroup of operators [2] of any Hamiltonian theory. 
Today the physical utility of light-front dynamics is transparent: hard scattering processes 
probe a light-cone correlation of bound states. Not surprisingly, then, many perturbative 
QCD applications can be treated on the light front, see e.g. [3]. Outside this realm, physics 
on the light cone has been extensively developed for non-perturbative QCD [4-7] as well as 
applied to nuclear physics [8,9]. 

In this thesis we consider applications of light-front dynamics to generalized parton 
distributions (GPDs). As GPDs enter into amplitudes for hard exclusive reactions, light- 
cone quantized fields make a natural appearance. GPDs are theoretically challenging objects 
to calculate. Even simple phenomenological parameterizations cannot be written down at 
will [as is the case of conventional parton distributions functions (PDFs)] due to the field 
theoretic constraints GPDs must satisfy. Ultimately we will generate a phenomenological 
parametrization for the proton's GPDs. Before doing so, however, we investigate from a few 
perspectives the difficulties inherent in describing these new distribution functions. We first 
motivate the physicality of the light-cone framework by considering deep-inelastic scattering 
and then discuss bound states in field theory and the front form of Hamiltonian dynamics 

1.1 Deep- inelastic scattering 

Deep-inelastic scattering (DIS) is the canonical hard scattering process. We shall therefore 
review DIS in some detail in order to introduce the relevance of light-front dynamics in pro- 
cesses amenable to treatment in perturbative QCD. Our treatment parallels the analysis 
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Figure 1.1: Diagrammatic depiction of the DIS process. An electron of momentum k is 
incident on a proton of momentum P. The leading-order electromagnetic process is virtual 
photon exchange between the electron and proton. As a result, the electron's momentum is 
altered to k' and due to the high virtuality of the photon, the proton breaks up into various 
reaction products collectively denoted by X. 

of [10,11]. We show how light-cone coordinates and the collinear approximation enable sim- 
plification of the DIS amplitude into a product of hard scattering coefficients and universal 
functions. The latter are the conventional parton distribution functions (PDFs). These 
distribution functions are momentum probabilities with respect to the quark wavefunction 
of the target defined along the advance of a wavefront of light. The analysis presented 
below will be extended later in Chapter |31 to deeply virtual Compton scattering (DVCS). 
Nonetheless, this Section gives practical motivation for the need to investigate bound states 
on the light cone. 

1.1.1 Bjorken limit 

In DIS off the proton, the process under consideration is electron^ scattering at large mo- 
mentum transfer, namely 



Here k, k' label the momenta of the initial and final electrons respectively, while P is the 
momentum of the proton. The reaction happens at such large momentum transfer that the 
proton breaks up, as it were, and subsequently forms new hadrons. We have used X to 
label any of the myriad of possible final states. 

We define the four-momentum transfered to the proton as q = k — k' . In the language 
of QED, this momentum is transfered via the virtual exchange of a photon between the 
electron and proton, see Figure ITTTl The offshellness of the photon is given by q'^ = —Q^, 
where is assumed to be very large and positive in the deep-inelastic limit. In accordance 
with the Heisenberg uncertainty principle, the larger the momentum transfer, the smaller 
the region probed. Hence one can view the electromagnetic process in DIS as resolving the 
structure deep within the proton. From the electron's perspective, however, the proton is 
severely Lorentz contracted and thus the virtual photon is only capable of resolving small 
sizes transverse to the beam direction. 

In evaluating the cross section, we have the freedom to choose both the beam direction 
and the overall frame in which to view the scattering reaction. With a suitable orientation of 
the beam direction and a convenient choice of reference frame, we can work simultaneously 
with P-*" = {P^, P"^) = and q"*" = 0. Now as is conventional, we introduce two vectors 

^One can also consider muon, positron, neutrino, etc. scattering off the proton (or other) targets. 
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Figure 1.2: The forward virtual Compton amplitude T^'^ (q, P). The dashed line denotes 
the cut that yields the imaginary part needed for the DIS cross section. 

and n'^ into which we decompose the remaining momenta. Explicitly these vectors are 

^^^ = ^(1, 0,0,1) (1.2) 

and 

n'^ = ^(1,0,0,-1), (1.3) 

where A is an unfixed parameter associated with the unexploited freedom to boost along 
the z-direction. Notice both vectors are lightlike: = = 0. Additionally their product 
\s p ■ n = \. 

In terms of these p and n vectors, the proton's momentum appears as 

P^^=p^^ + —n'^, (1.4) 

so that = M^. In the proton's rest frame, A = M/^/2. In this frame we define the 
energy of the virtual photon = u and hence P ■ q = Mu. The photon's momentum in the 
proton rest frame is thus 

q^^ = {iy,0,0,-VW+^) (1.5) 

because the virtuality is q'^ = —Q^. 

Lastly we define the Lorentz invariant Bjorken variable x 

(1.6) 



2P-q 



In the proton rest frame, we have x = Q"^ /2Mv. In the limit of elastic proton-electron 
scattering x is fixed by momentum conservation, x = 1. For a general inelastic process, x 
is kinematically restricted to S (0,1]. Notice in particular this implies u'^ ^ Q^. In DIS, 
one refers to the Bjorken limit in which Q'^ ^ oo while x remains fixed by the kinematics. 
In this limit the doubly differential cross section for the reaction in Eq. (|1.1|) has the form 

^':L^'^ik,q)W,A<l,P), (1-7) 



dE'dQ' MQ^ E 

where we have included only the leading-order electromagnetic contribution from one photon 
exchange. Here E' is the final electron energy and dD.' its differential solid angle, while E is 
the initial electron energy. L'^'^{k,q) is the leptonic tensor representing the QED amplitude 
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squared for virtual photon emission off the electron line and a^; is the electromagnetic 
fine-structure constant. The unknown quantity in the cross section is the hadronic tensor 
W^'^ {q, P) which denotes the amplitude squared for the virtual photon-proton interaction. 
The hadronic reaction products X, however, are not observed in inclusive DIS and are 
accordingly summed over in W^'^{q, P). Thus we utilize the optical theorem 

W'^''{q,P) = -^Qm[T>^-'{q,P)] (1.8) 

to relate W^'^{q, P) to the forward virtual Compton amplitude T^'^{q,P) that is depicted 
in Figure E21 In terms of the Heisenberg current operator J^{y), T^'^(q,P) can be written 
as the diagonal proton matrix element of the time-ordered product of two currents, namely 

T^'•'{q,P) = i y Ae^^'^(-P|r{j^(y)J''(0)}|P). (1.9) 
1.1.2 Collinear approximation 

Above we have defined the kinematics of DIS and written the experimentally observed cross- 
section in terms of the forward virtual Compton amplitude T^'^ (q, P). We now proceed to 
isolate the leading contribution to T^'^{q, P) in perturbative QCD.^ Large momentum q flows 
through the diagram depicted in Figure E21 Since T^"{q,P) is a Fourier transform, most 
contributions are averaged out as q^ — oo. Indeed the only non- vanishing contributions 
that remain arise from singularities in the product of current operators as y"^ 0. 

Due to the nature of the QCD running coupling, at large momentum transfer the highly 
virtual photon couples to asymptotically free quarks. Asymptotic freedom justifies the 
use of the impulse approximation for the product of currents and the so called handbag 
dominance emerges (see Figure II. 3() . More rigorously, the handbag mechanism picks out 
the leading singularity as ^ in Eq. (|1.9j) . The handbag contribution to W^'^{q,P) is 
thus^ 

along with an additive similar contribution from the crossed handbag. The matrix element 
M{k, P) in the above expression is the forward free-quark proton scattering amplitude given 
by 

M^p{k,p) = I dSe"'-y{p\T\lPf,{o)My)}\P), (1-11) 

where a and /3 are Dirac indices. 

Lastly we need to isolate the dominant contribution from the matrix element M(k, P) 
in the Bjorken limit. To do this we decompose the active quark momentum k^^ in terms of 
the light-cone vectors and n^, viz. 

k^ = {k-n)p^ + {k-p)n^ + (1.12) 

^Throughout this work, we shall work only to leading order in the strong coupling. Important points 
concerning 0{a.s) corrections are contained in footnotes such as this. For a comprehensive discussion of 
perturbative QCD corrections to DIS, see e.g. [11]. 

^ Quark masses have been neglected in the propagator with large momentum flow. The restriction 
m^/Q^ <^ 1 automatically follows from /Q^ <g 1, which, as we shall see, is the natural scale sepa- 
ration condition in DIS. 



ik + . 



+ ie 



-f^Mik.P) 



(1.10) 
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Figure 1.3: Graphical depiction of the handbag mechanism in DIS. Due to asymptotic 
freedom, the forward virtual Compton amplitude is dominated by free quark scattering. 
There are two contributions: the handbag shown on the left, and the crossed bag shown on 
the right. Figure adapted from [12]. 



and analyze the situation in each momentum channel. The situation is the simplest with 
respect to the transverse component. There, the photon carries no transverse momentum 
and so the active quark's transverse momentum is unchanged after interaction. In the p'^ 
direction, the photon adds q ■ n = —x to the active quark's momentum forcing the final 
momentum to he k ■ n — x. Lastly in the direction, the photon adds q - p k, oo and the 
rebounding quark has nearly infinite momentum in the channel. Thus relative to this 
direction the struck quark has effectively zero momentum, i.e. k"*" ~ and k ■ n — x k, {). 
The latter we enforce by inserting 1 = / dx8{x — k ■ n). Our analysis leads to the collinear 
approximation to the active quark's momentum 



^{k- n) pi", 



(1.13) 



because the transverse momentum is effectively zero and after interaction with the photon 
the n'^ component is infinite, which leads to suppression of the matrix element. In the 
collinear approximation, we have the resulting contribution to W^'^{q, P) 



r 



x^ + i 



{xp + g)^ + ie 



YM^pix) 



(1.14) 



Above we have reduced the free-quark proton Green's function M(k, P) into a function 
of only^ Bjorken x 

/d^k 
^S{x-k-n)M^^{k,P). (1.15) 



In this form, the integral over y in Eq. (|1.11|) can be performed trivially. We arrive at 



Ma^i^) = I ^e^^^{P\i;f,mU^n)\P), 



(1.16) 



^This statement and the factorization in Eq. H1.14^ are modified by 0(qs) corrections. These arise from 
perturbative gluon interactions which lead to renormalization scale dependence, i.e. Ma/s = Mai3{x, /i'^). 
Physically jj, is the scale at which the partons are resolved. Consequently the notion of a parton is renormal- 
ization scale and scheme dependent. The Compton amplitude W^" [q, P) remains scale independent because 
the 0{as) corrections to the hard scattering kernel in Eq. 11.1411 exactly compensate for the /i dependence 
in the quark-proton Green's function. 
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here we have used An to denote evaluation at the position = An^. Thus the DIS 
scattering cross-section depends upon a correlation function containing two quark fields 
with a lightlike separation (above the resulting separation is = A/\/2A, = 0, and 
y-*- = 0). This means that the relevant part of the quark propagator in DIS is only that 
along the light cone. Consequently the amplitude depends upon a correlation of proton 
states along the light cone and inherent to DIS is the system's response along the advance 
of a wavefront of light. Furthermore, as a result of the collinear approximation Eq. 
the active quark's plus momentum is a fraction of the proton's, i.e., fc"*" = xP~^ . 

For unpolarized DIS, we can decompose the Dirac structure of the light-cone correlator 
M{x) into and vectors 

Map{x) = \i>ap fl{x) + ^MV«/3 U{x). (1.17) 

Taking the trace of M{x) with |^ and ^ separately, we deduce 

J ^e'^-{P\-^{0)r4^{Xn)\P) = 2 [/i(x)p^ + M^Mx) n^^] . (1.18) 

Now since both and are zero, and n-p = 1, contributions to the DIS cross section from 
the /4(x) structure function must come out proportional to ^rf^ix) at the end of the day. 
One can neglect this structure due to the power-law suppression in the Bjorken limit. ^ This 
simplification is a peculiarity of the light-cone power counting^ because the overall matrix 
element has one mass dimension but two irreducible structure functions. 

To summarize, we have reduced the virtual forward Compton amplitude in the Bjorken 
limit to one unknown structure function fi{x). In terms of proton matrix elements, this 
structure function is 

fiix) = \j ^e^^-(P|V^(0)7+^(An)|P). (1.19) 

Conventionally /i(x) is referred to as a parton distribution function (PDF). This is tech- 
nically only true in light-cone quantization in light-cone gauge (see Appendix [HI) • Never- 
theless, to adequately model and physically interpret the PDF as a distribution one needs 
to know properties of the bound state (here the proton) on the light front. Such properties 
are encoded in the light-cone wavefunctions. We shall give a preliminary introduction to 
these light-cone wavefunctions in Section [1.21 

1.1.3 Power counting with a twist 

Above we have reviewed DIS in a manner intrinsically tied to the light cone. This will 
guide us throughout the first few Chapters as we investigate properties of bound states 
and wavefunctions on the light cone. There exists an alternate, covariant way to view 
DIS structure functions that will be required in the later Chapters when we take up the 
calculation of double distributions (DDs). Moreover this covariant representation will allow 
us some insight as to the peculiarities of light-cone power counting. Instead of re-deriving 

^This statement is not altered by perturbative QCD corrections. Indeed Eq. 11.181 is modified to include 
scale dependence, but the structure functions fi{x,fi^) and f4,{x,^^) do not mix under evolution. 

^Notice the vector has only a minus component and hence projects out the plus component of a vector 
dotted into it, while the n'' vector has only a plus component and projects out the minus component. In 
older terminology, p'^ gives one the good components, while gives one the bad components. 
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DIS amplitudes in the new framework, we shall take our result above Eq. (|1.19|) and merely 
extract from it a new perspective. 

Consider the n}^ moment of the PDF fi{x). We can write it in the form 



T \ n 



(P|V(0)7+V'(An)|P). (1.20) 



Performing n-integrations by parts enables evaluation of the j;-integral, which finally ap- 
pears as f l^e*'^^ = (5(A). Thus we have 



dxx"/i(x) = -n^n^,...n^„(P|V^(0)7^(i5^i)---(i9^")V'(0)|P). (1-21) 



This suggests we consider the tower of local operators O^") defined by 

0m.../.n ^ ^(o)^{/^iB^i . . . iD^''^}^p{0), (1.22) 

where the action of ^'"^ on Lorentz indices produces the symmetric traceless part of the 
tensor. The derivative D^^ = d^ + igA^^ is the QCD gauge covariant derivative. In extracting 
the form of O^") in Eq. (TT^ from Eq. (fTTTT) . we made a number of trivial modifications 

for free. Firstly we made the partial derivative act symmetrically d = d — d since the ip 
field has no A dependence in Eq. (|1.19|) . Next we upgraded the partial derivatives to gauge 
covariant derivatives to render the tower of operators C^") gauge invariant. Symmetrization 
with respect to all Lorentz indices was also done at no cost since the prefactor of {n + l)-n^ 
vectors in Eq. H1.21() filters out the symmetric part of the tensor operators. Notice that 
symmetrization allows us to identify the operator O^^^ with the quark part of the QCD 
energy-momentum tensor. Lastly tracelessness, which allows us to classify O^") irreducibly 
as spin-(n -|- 1) operators, is a free choice since any trace terms are proportional to g^^'■^^i 
and must vanish when contracted with n^.n^^ since = 0. 

Now we merely observe that because of Lorentz invariance, proton matrix elements of 
the O^") operators have the form 

{P\0^'^'^■■■^'^\P) = 2anP^^ P^^ ■ ■ ■ P^''^ , (1.23) 

where a„ is a Lorentz scalar.'' Recall that n ■ P = n ■ p = 1. Thus we have 

dxx"/i(x) = a„. (1.24) 



Moments of PDFs are proton matrix elements of local operators 0^"'\ A natural question 
to ask, however, is why there are infinitely many such operators. Power counting seems to 
suggest that more P'^'s in Eq. (|1.23|) means more suppression in the Bjorken limit. 

The resolution of this paradox (and also the peculiarities of light-cone power counting) 
involves fully appreciating the nature of the Bjorken limit. Let us denote an arbitrary 
operator of mass dimension d and spin s as 0{d,s). Unpolarized proton matrix elements 
of this operator will pick up s-powers of P due to Lorentz covariance. But the net effect on 

^In general a„ = an{^^). As an alternative to our above discussion of DIS, we can view the scale 
dependence as arising from the renormalization of the composite operators in Eq. I|1.22|l . Notice ^^oo = 

since O'"' is an operator with vanishing anomalous dimension. 



7 



V 



Figure 1.4: Diagrammatic representation of the Bethe-Salpeter equation. The blob repre- 
sents the vertex function V. 

the Compton amplitude T^'^(q,P) must come from contracting these P's with q^s. Hence 
contributions to the physical cross section from 0{d,s) are proportional to 

where the equality utilizes the definition of Bjorken x, Eq. (|1.6|) . Thus the relevant power 
counting in DIS comes with a twist 

T = d-s. (1.26) 

The operators with the lowest twist r give the leading contribution to DIS.*^ The opera- 
tors are twist-2 contributions and are thus all relevant for DIS. The function f^ix) 
encountered above is a twist-4 contribution. In summary, the calculation of model quark 
distributions from their moments is an interesting study of the interplay between the light- 
cone formalism and covariance. We shall take this much further in Chapters |1] and |SJ where 
we obtain model double distribution functions. 



1.2 Bound states 

Above we have reviewed DIS from two perspectives: a covariant development formulated 
in terms of matrix elements of local twist-2 operators, and also from the perspective of 
dynamics along the light cone. For applications throughout the rest of this work, we need 
to introduce light-cone wavefunctions and investigate how the features of covariant field 
theory manifest themselves in the many-body light-cone wavefunctions. To this end, we 
first discuss in Section 11.2.11 two-body bound states in covariant field theory. Next in 
Section [1.2.21 we provide an exactly soluble model wherein one can compare the covariant, 
light-cone and familiar instant form wavefunctions. This involves solving for the model's 
wavefunction in the rest frame and then boosting to an arbitrary frame. 

1.2.1 Bethe-Salpeter equation 

In terms of fully covariant operators, the Lippmann-Schwinger equation for the two-particle 
transition matrix T appears as 

T = V + VGT. (1.27) 

*The renormalization of these operators does not modify the twist expansion because operators of different 
twist do not mix under evolution. One has the additional feature, however, that operators of the same twist 
(and quantum numbers) do mix under evolution, e.g. to 0[ois) one must consider mixing of O'"' above with 
the twist-two gluon operators 
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Figure 1.5: Bethe-Salpeter equation for a point interaction. The state is bound by the 
infinite chain of bubbles. 

Above, V is the irreducible two-particle scattering kernel and G is the completely discon- 
nected two-particle propagator (which is merely the product of two single-particle propaga- 
tors). A pole in the T-matrix (at some = M^, say) corresponds to a two-particle bound 
state. Investigation of the pole's residue gives an equation for the bound state vertex T 



The amplitude ^ is defined as GT and hence satisfies a similar equation, the Bethe-Salpeter 
equation [13-15] 



Following [16], it is convenient to denote quantities able to be rendered in position or 
momentum space with bras and kets. We will employ this notation only for quantities 
that have been stripped of their overall momentum-conserving delta functions, for example 
T{k,R) is defined by T(k,R) = {k\Tji). Here we have used R as a label for the bound 
state for which R^ = M^, and k labels the four momentum of the first particle. The same 
is analogously true for ^{k,R). For the disconnected two-particle propagator, we define 
first G{R) which is the disconnected propagator of total momentum R defined through 
the relation {R'\G\R) = (27r)^5(^) (i?' - R)G{R). Additionally we remove the momentum 
conserving delta function between initial and final states of particle one in the discon- 



nected Green's function to arrive at G{k,R), i.e. {k\G{R)\p) = {2TT)'^5'-^\k - p)G{k, R). 
Using {R'\V\R) = {2Tr)^6^^\R — R')V{R), we see the momentum-space kernel V{k,p;R) is 
V{k,p;R) = {k\V{R)\p). Thus rendered in momentum space, the Bethe-Salpeter equation 
(jl 29() reads (see Figure 



Armed with the Bethe-Salpeter amplitude ^, one can calculate field-theoretic bound- 
state matrix elements by taking the appropriate residues of four-point Green's functions. 
These matrix elements may ultimately require knowledge of higher-point functions which 
then must be solved for consistently in the same dynamics. The Bethe-Salpeter amplitude ^ 
is in some ways the covariant analogue of the Schrodinger wavefunction. While the features 
of relativistic field theory (in particular: particle creation and annihilation, retardation 
effects. . .) make the exact analogy impossible, in the non-relativistic limit, one can show 
that the BSE reduces to the Schrodinger equation. 

1.2.2 Toy wavefunction and boost 

Above we have derived the covariant BSE for two-body bound states. In this Section, we 
consider a toy model for the BSE that is exactly soluble. We shall for simplicity consider a 
scalar bound state composed of two scalar particles. The solution will enable us to compare 



r = VGT. 



(1.28) 



^ = GV^. 



(1.29) 




(1.30) 



9 



and contrast instant form dynamics and light-front dynamics all while maintaining exact 
covariance. The instant and front form of dynamics are reviewed in Appendix IHl 

One can obtain the simplest soluble BSE equation by choosing a point-like interaction 
for V{k,p;R) in Eq. ()1.3U|) . namely V{k,p;R) = ig, where g is a coupling constant. The 
two scalar particles that make up the scalar bound state thus interact infinitely many 
times according to the BSE to bind the state. For a point-like interaction, a bubble chain is 
generated by the BSE and is shown in Figure [T31 With this choice of interaction, the bound 
state equation simplifies tremendously. Since the kernel is independent of momentum, the 
only p dependence that remains in Eq. (|1.3U|1 is in "^{p, R) and this quantity is subsequently 
integrated over all p. The integration merely produces some number that can be absorbed 
into the overall normalization of the wavefunction (along with the coupling constant in the 
kernel). Thus we are left with the solution 

■^{k,R) =iG{k,R), (1.31) 

where the overall constant is a matter of taste. Comparing with Eq. (|1.28jl . we can also 
trivially solve for the bound state vertex function: T{k,R) = i. We shall see later that 
a point-like vertex is a convenient choice for phenomenological applications because the 
crossing properties are trivial. The Bethe-Salpeter equation for r(A;, R) determines the 
mass = R^ of the bound state via the consistency equation 

1 = ^9/ ■0^G{k,R). (1.32) 

Now for simplicity we shall neglect self-energy contributions in the propagators and 
choose the single particle propagator to have the basic Klein-Gordon form. The two-particle 
disconnected propagator is a product of these Klein-Gordon propagators and hence by virtue 
of Eq. (|1.31|) the covariant wavefunction is 

^{k, R) = -i[k^ -n? + ie\-^[{k - Rf - + ie]'^ . (1.33) 

Here we have labeled the constituent mass by m. If we imagine taking the Fourier transform 
of this wavefunction, we have a four-dimensional analogue of the Schrodinger wave function. 
There is an important distinction, however. We also know the time dependence of the 
wave function — the time evolution governed by the Hamiltonian operator is automatically 
included because of the necessity of covariance. Moreover, we know from the Poincare 
algebra that there are other dynamical operators besides the energy^. As to which operators 
are Hamiltonians and thus which are kinematical depends upon the form of dynamics chosen. 
The instant and front forms are reviewed in AppendixEl We shall consider first the familiar 
instant form. 

In the instant form of dynamics, the energy and the Lorentz boosts are the dynamical 
operators. The initial conditions are specified on the boundary = 0, as in classical and 
quantum mechanics. The time evolution of the system is governed by the energy operator 
and involves non-trivial effects due to the system's dynamics. Moreover the relativistic 
properties of the system also are complicated by the dynamics in the instant form. This is 
reflected in the fact that the generators of Lorentz boosts are Hamiltonians. If we write out 
the Fourier transform of the Bethe-Salpeter wavefunction 

^{x,R) = [ —e"-''^{k,R), (1.34) 
J 27r 

^Dirac advocated the use of "Hamiltonian" for all operators that are dynamical, of which the energy 

operator is one. 
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we see that the instant form wavefunction can be obtained by integration over the energy 
k^. Thus at time = and in the rest frame of the bound state, we have 

— ^ik,R) . (1.35) 

Given our solution to the BSE Eq. TTTil . we can carry out this projection onto the 
initial surface. The integration can be done using the residue theorem bearing in mind the 
four poles of the integrand 

A:° jj = ± + m2 =f ie (1.36) 
k^-^ = M± \/k2 + m2 =F ie. (1.37) 

The instant form wavefunction is 2iTi[Kes{k^) + Res(/i;^)]. Explicitly we have 

V^O(k, 0) = ^ } (1.38) 

and overall multiplicative constants have been neglected. Notice the wavefunction is mani- 
festly rotationally invariant. This is indicative of the kinematic nature of the generators of 
rotations in the instant form. 

From the point of view of relativistic quantum mechanics, one might wonder what the 
wavefunction is in an arbitrary frame. For example, in considering electromagnetic form 
factors, the final state is boosted relative to the initial one and thus to model the form 
factors phenomenologically, one requires knowledge of the boost. To deduce the boosted 
wavefunction, one must solve the complicated dynamical equation 

|^0(k,R))=e™|V°(k,O)), (1.39) 

where K are the generators of Lorentz boosts. For any three-vector k, let us define -E(k) = 
\/k2 -|- m2. The boosted wavefunction which solves Eq. ()1.39|) can be compactly written as 



^°(k,R) 



1 



^(R) + ^(k) - E{R - k) I ^(k) [E(R) + ^(k) + E(R - k)] 

^ E{-R - k) [E{K) - E{k) - E{-R - k)] } ' ^"^'^^^ 

To obtain the boosted wavefunction, we chose not solve Eq. (|1.39j) . We merely returned 
to the covariant Bethe-Salpeter wavefunction and evaluated it at x*^ = using Eq. H1.35() 
in an arbitrary frame. The non-transparent relation between the rest frame wavefunction 
Eq. H1.38() and the boosted wavefunction Eq. H1.4U() indicates the dynamical nature of the 
boost. 

In the front form of dynamics, one is interested in the properties of physical states along 
the advance of a wavefront of light. The objects of front form dynamics are the light-cone 
wave functions which are projections onto the initial surface = 0. In analogy with the 
instant form, one refers to as light-cone time, and its Fourier conjugate k~ as light-front 
energy. In the front form, the energy is a Hamiltonian along with two rotation operators 
corresponding to two independent rotations of the wavefront of light. In contrast with the 
instant form, light-front Lorentz boosts are kinematical. 
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By virtue of Eq. H1.34[) , the hght-cone wavefunction is 

/rlk~ 
-^^{k,R), (1.41) 

where x is the fraction of the plus momentum carried by the first particle, x = k^/R'^ 
and the relative transverse momentum is k;^j = k"*" — xR"*". Evaluation of the integral is 
standard and actually simpler than in the instant form because, on the light-cone, there are 
only two energy poles. Taking the appropriate residue, the light-cone wave function for our 
toy model is 

V'(x,k^J = |^^)f^ Irr-^. (1.42) 

x{l—x) 

Notice the wave function maintains only a cylindrical symmetry and further that the relation 
of the rest frame wave function to the boosted wavefunction is all compactly written above in 
terms of k^j. These facts are of course indicative of the kinematic subgroup of the Poincare 
generators (the subgroup includes rotation within the initial surface and light-front boosts). 
Accordingly the full rotational symmetry of the rest frame wavefunction is not manifest, c/ 
Eq. (OrI) and Eq. (OTT) . 

One should now inquire as to how the two wavefunctions ip^ and ip are related to each 
other. We can try to turn the light-cone wavefunction into the rest frame wavefunction. In 
the literature, this is accomplished by introducing an auxiliary variable . So that k^ has 
a physical interpretation in terms of the z-component of momentum, one uses the two-body 
center of mass relation 

k+ I k^ , , 

^ = = ^ + , =. (1.43) 

^ ^ 2^\i\ + (A;3)2 + m2 

Here R'^ is taken to be the sum of the free constituents' plus momentum. There is no 
binding effect in the total plus momentum because longitudinal translations are kinematic 
in the front form of dynamics. Inverted this relation between x and k'^ reads [17] 

Simple algebra yields the relation 




from which we deduce 



This bears a resemblance to the instant form wavefunction in the rest frame Eq. H1.38() . If 
we include the Jacobian 

9fc3 2S(k)3 ^ ' 

as a multiplicative prefactor, then the wavefunctions agree. There is, however, no reason 
why this should work: only a factor of \ ■§§^ is justified to preserve the norm. When 
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one considers not the wavefunctions themselves but the relevant amplitudes, enforcing the 
constraint of Poincare invariance on the physical amplitude can lead to the proper factors 

to include, e.g. for form factors in the Drell-Yan frame the proper factor is The 
variable transformation Eq. H1.44() is additionally convenient for the evaluation of integrals. 
Notice, the light-cone wavefunction itself can be evaluated in the rest frame, however, it is 
not connected to the rest frame wavefunction in the instant form by a non-singular boost. 
For amplitudes involving a fixed number of particles, the above correspondence Eq. (|1.44|) 
can be used along with constraints from Poincare covariance to determine the relation to 
the instant form rest-frame wavefunction. In general, however, there is no relation between 
field theoretic amplitudes and Poincare covariant representations of relativistic dynamics. 

While there is no general way to unboost the light-cone wavefunctions, the instant form 
wavefunction ^'^(k, R) can be boosted to the infinite momentum frame (IMF). As a result of 
this procedure, one recovers the light-cone wavefunction i/>(x,k-'"). Originally field theories 
were re-formulated in the infinite momentum frame [18-21] as a way to simplify covariant 
perturbation theory. The connection to light-cone quantization was made later. 

To boost the wavefunction to infinite momentum, we first simplify matters by choosing 
R-L = 0. In the IMF, ^ oo while i?^ = is kept fixed. Thus one has 



^ ^ 2(ii3)2 + • ' 



(1.48) 



Notice in this hmit, we have i?+ ^ V2R^ oo while ^ M'^/2V2R^ 0. 
Working with the constituent's momentum, we use the parametrization 

k = 'k^ +xR^z. (1.49) 

The energy factors then simplify in the IMF as 

£(k) = rf + ^t-^ + ,,, (1.50) 
£(R-K) = + + (1.51) 

Carrying out the boost to infinite momentum on the wavefunction in Eq. ()1.4U|) yields the 
IMF wavefunction which we denote V''^(k, oo). Using Eqs. (|1.5())1 and (|1.51j) we find 

x{l—x) 

which since i?^ oc R^ is just the light-cone wavefunction evaluated for R-*" = 0. Unboosting 
the wavefunctions from the IMF is generally impossible but can be done in the context of a 
few specific amplitudes for which a correspondence between relativistic quantum mechanics 
and field theory can be made. 

In summary, this simple covariant model for the BSE Eq. H1.31() has allowed us to 
explore both the instant and front form wavefunctions. We saw how the structure of these 
wavefunctions is related to the respective kinematic subgroups of the Poincare algebra. 
Moreover, a fully covariant starting point allowed us a simple way to correctly formulate 
the three-dimensional dynamics. We shall employ this philosophy throughout the work. 
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Chapter 2 

Light-Front Bethe-Salpeter 
Equation 

Above we have seen the utiHty of hght-front dynamics in processes at large momentum 
transfer. Additionally we investigated the properties of light-cone wavefunctions in a sim- 
ple exactly soluble model for the covariant BSE. This Chapter concerns bound states of two 
particles in the light-front formalism, specifically of interest are current matrix elements be- 
tween bound states. We approach the topic, however, from covariant perturbation theory in 
order to dispel rampant misconceptions about bound states on the light front. As demon- 
strated by the tremendous undertaking of [22,23], one can derive light-front perturbation 
theory for scattering states by projecting covariant perturbation theory onto the light cone, 
thereby demonstrating their equivalence — including the delicate issue of renormalization. 
As to the issue of light-front bound states, a reduction scheme for the Bethe-Salpeter equa- 
tion recently appeared [16,24] that produces a kernel calculated in light-front perturbation 
theory. This reduction scheme makes formal earlier observations about the connection be- 
tween the Bethe-Salpeter equation and the light-front bound state equation [25]. For the 
purpose of simplicity, in this Chapter we consider only bound states of two scalar particles 
interacting via the exchange of a massive scalar in the (1 -|- l)-dimensional ladder model. 

Our main consideration is to extend the reduction to current matrix elements to in- 
vestigate carefully valence and non-valence contributions in the light-front Bethe-Salpeter 
formalism. To do this, we calculate our model's form factor. As calculations in the light- 
front reduction are covariant only when summed to all orders, our results violate Lorentz 
symmetry and we choose to extract the form factor from the plus-component of the electro- 
magnetic current in order to make contact with the Fock space representation. Moreover, 
in calculating (1 -|- l)-dimensional form factors we cannot choose a frame of reference where 
Z-graphs vanish. This enables us to investigate their contribution, which in (3 -|- 1) dimen- 
sions has a variety of applications such as to generalized parton distributions, which we will 
detail in Chapter |3J Z-graph contributions haunt light-front dynamics since non- valence 
properties of the bound state are involved, so that valence wavefunction models cannot be 
utilized directly. On one hand, the light-cone Fock representation provides expressions for 
the Z-graph contributions in terms of Fock component overlaps which are non-diagonal in 
particle number [26,27]. While on the other hand, vertices which cannot be related to the 
valence wavefunction (coined as non-wavefunction vertices in [28]) appear in the light-front 
Bethe-Salpeter formalism. A variety of ways have been proposed for dealing with these 
non-wavefunction vertices [29,30] including attempts to model the covariant vertex [31,32], 
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or (when possible) estimating the contribution from higher Fock states [33]. In a Poincare 
covariant framework, Z-graph contributions involve a dynamical light-front rotation and 
represent a similar formidable problem. 

Below we show that non-wavefunction vertices are supplanted by contributions from 
higher Fock states in light-front time-ordered perturbation theory (provided the interaction 
has light-cone time dependence). In essence contributions from non-wavefunction vertices 
are reducible and should only be used when the interaction is (or is approximately) in- 
stantaneous^ . This constitutes a replacement theorem for non-wavefunction vertices which 
trivially extends to (3-1-1) dimensions. When one works from covariant perturbation theory, 
the coupled tower of Dyson-Schwinger equations gives rise to the light-cone Fock compo- 
nents [34,35]. Thus when the Green's functions are treated properly and the light-front 
BSE is derived, instead of postulated, non-wavefunction contributions are absent. 

The organization of this Chapter is as follows. First in Section 12.1.11 we present the 
issue of non-wavefunction vertices and energy poles of the Bethe-Salpeter vertex focusing 
on a typical light-front constituent quark model as an example. For such models, non- 
wavefunction vertices are required to express the form factor. We find the commonly used 
assumptions in quark models necessitate vertices not only without energy poles but without 
energy dependence. Hence covariance is generally lost. Next in Section [2.1.21 we review the 
reduction of the Bethe-Salpeter equation presented in [16] focusing on the energy poles of 
the vertex. We derive an interpretation of the reduction as a procedure for approximating 
the poles of the vertex. Additionally in Section f2. 21 we show how the normalization of the 
covariant Bethe-Salpeter equation turns into a familiar many-body normalization in the 
light-front reduction. The normalization condition resembles a diagonal matrix element of 
a pseudo current and has contributions from higher Fock states. We calculate the explicit 
normalization condition for the ladder model at next-to-leading order in perturbation theory. 
The connection to the familiar many-body Fock state normalization is made in Appendix 
IeI In the following Section 1)2. 3() . we construct the current to be used with the reduced 
formalism. In Section f2. 41 the (1 -|- l)-dimensional ladder model is presented and in Section 
12.4.21 we compare the calculation of the form factor for the model using two different paths 
to the reduction. The comparison allows us to see when non-wavefunction vertices can be 
efficiently used. Lastly we summarize our findings in Section [2. 51 

2.1 Light-front reduction 

To derive a bound state equation on the light front, we must investigate the reduction of the 
covariant BSE Eq. (|1.29|) to the hypersurface = 0. Following the toy wavefunction ex- 
ample presented in Section [1.2. 21 we must integrate over the light-cone energy to obtain the 
proper three-dimensional boundary value. The clear first step to handling such a reduction 
approximately is to make simplifying assumptions about the dependence on the light-cone 
energy to enable the projection. Near independence from light-front energy corresponds to 
an instantaneous approximation with respect to light-front time. This is where we begin 
our discussion. 

^We shall often refer to interactions and vertices merely as instantaneous if they are independent of 
light-cone time, or equivalently light-cone energy. Similarly we are not careful about referring to light-cone 
time and light-cone energy as time and energy, respectively. 
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2.1.1 Poles of the Bethe-Salpeter vertex 

To introduce the reader to non-wavefunction vertices and instantaneous approximations, 
we focus on light-front constituent quark models. We start by writing down the covariant 
equation for the meson^ vertex function T. In (1 + 1) dimensions, it satisfies a simple BSE 
(see Figure [O]) : 

r{k,R)=z I ^V{k,p)^Bs{p,R), (2.1) 

in which we have defined the Bethe-Salpeter wavefunction ^bs as 

^Bs{k,R)=G{k,R)r{k,R), (2.2) 

with the two-particle disconnected propagator G{k, R) = d{k)d{R — k). For scalars of mass 
m, the renormalized, single-particle propagator d has a Klein-Gordon form 

"^^^^ " (F - m2)[l + {k^ - m2)/(A;2)] + ie ' ^^'^^ 

where the residue is i at the physical mass pole and the function f{k'^) characterizes the 
renormalized, one-particle irreducible self- interactions. To simplify the comparison carried 
out in Section [2.41 we shall ignore /(A;^). Above V is the irreducible two-to-two scattering 
kernel which we shall refer to inelegantly as the interaction potential. 

Now we imagine the initial conditions of our system are specified on the hypersurface 
= 0. As we saw for DIS, correlations between field operators evaluated at equal light- 
cone time turn up in hard processes for which this choice of initial surface is natural. In 
order to project out the initial conditions (wavefunctions, etc.) of our system, we must 
perform the integration over the Fourier conjugate to x"*", namely k~ . For instance, our 
concern is with the light-front wavefunction ip{x) defined as the projection of the covariant 
Bethe-Salpeter wavefunction onto = 0, 

/dk~ 
—^Bs{KR), (2.4) 

with X = k^ / R'^ . The overall multiplicative factor in the definition of ip{x) is an arbitrary 
choice. 

Looking at Eq. ()2.2|) . in order to project the wavefunction exactly, we must know the 
analytic structure of the bound-state vertex function. If the vertex function r(A;, R) had 
no poles in /c~, then our task would be simple: the light-front projection of ^ bs would 
pick up contributions only from the poles of the propagator G{k,R). Next we observe 
from Eq. I\'2.1\i . that the k~ dependence of the interaction V{k,p) must give rise to the k~ 
poles of the vertex function T(k,R). Hence an instantaneous interaction gives rise to an 
instantaneous vertex and a simple light-front projection (see, e.g., [29]). 

On the other hand, constituent quark models often assume a less restrictive simplifica- 
tion of the analytic structure of the vertex (see, e.g., [32]) in order to permit the light-cone 
projection. We shall show that this assumption along with the presumed covariance of 
the quark model implies instantaneous vertices. For any momentum p, let us denote the 
(1 + l)-dimensional on-shell energy p~^ = m'^/2p~^ . The propagator G{k, R) has two poles 

{h- — k- — H 
-"-^ - . (2.5) 

k, =R- + {k-R)-^-i^^. 



^Here we use meson to denote a general two-particle bound state. Mesons just happen to be convenient 
examples. 
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p (b) P' 

Figure 2.1: Covariant diagram for the electromagnetic form factor. Labels a, b and c denote 
subscripts used for on-shell energy poles. 



Notice that although we work in (1 + 1) dimensions, the results are trivial to extend to 
(3 + 1) dimensions because the imaginary parts of poles have precisely the same dependence 
on (only) the plus-momenta. This remark applies not only to this Section but to the rest 
of this Chapter. 

In constituent quark models that attempt to find a field theoretic basis, F is assumed to 
have no poles in the upper-half /c~ -plane for < x < 1. Since the light-front wavefunction 
is proportional to 6[x{l — x)], in light of Eq. ()2.5|) we further require any possible poles of 
F to lie in the upper-half plane for x < and in the lower-half plane for x > 1. With these 
restrictions Eq. (|2.4|) dictates the form of the constituent quark wavefunction 

^(x) = D^Jv{x\M^)T{h, R)9[x{l - x)], (2.6) 

where we have defined the Weinberg propagator as 

D^N{x\M'^) = i — (2.7) 

x(\—x) 

and used the abbreviation T{kii,R) to denote evaluation at the energy pole k~ = ap- 
pearing in Eq. 1)2. 5(1 . 

When we calculate the (elastic) electromagnetic form factor for these constituent quark 
models (see Figure ITT]) . we are confronted with more poles. Let A"*" = — C-P^ < be the 
plus-component of the momentum transfer between initial and final state mesons. The form 
factor is then 

(2x - C)F(A:, P)F* (A; + A, P + A) d^k 



Fit)^J [^2 _^2][(p_^)2 _^2]p + ^)2 _^2] (2-8) 
The /c~-poles from the propagators are k~ , k'^ defined in ()2.5() and k^ = — A~ + (/c + A)^^ — 

ie 
x-C- 

Given this pole structure, the contributions to F{t) are proportional to 9[x{l — x)]. In 
the region x < C,, closing the contour in the upper- (lower-) half plane will enclose possible 
poles of the final (initial) vertex. As one can see from considering the region x > the 
form factor can be determined solely from Res(A;^). In (3 + 1) dimensions, where we are 
free to choose frames in which C = 0) this becomes the only contribution to the form factor. 
But if ^ 7^ 0, the additional poles from vertices in the region x < C are required by Lorentz 
invar iance. 

If one advocates no modification to the pole structure of Eq. ()2.8|) due to the vertices, 
then by closing the contour in the lower-half plane, one picks up the residue at k~ without 
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any contribution from poles of the initial-state vertex. Such poles cannot lie in the upper- 
half plane for x < since the form of ip would be both frame dependent and contrary to that 
of Eq. H2.6|) [the (3 -|- l)-dimensional version of which is employed by all light-front quark 
models]. Thus one is actually assuming there are no poles of the Bethe-Salpeter vertex, 
which we originally pointed out in [30]. 

Returning to the definition of the wavefunction ()'2.4() under the premise of a vertex 
devoid of poles, we now find T(ki,,P) = T{ka,P), which we shall call pole symmetry. This 
symmetry is essential for making contact with the Drell-Yan formula [36], since for x > C, 
both initial- and final- state vertices may be expressed in terms of wavef unctions. When 
X < C) however, the final-state vertex becomes r*(A;a + A, P + A) which cannot be expressed 
in terms of 'tp*{x') (where x' = jE^) even with pole symmetry. Such a vertex we refer to as a 
non-wavefunction vertex. Understanding and dealing with such objects from the perspective 
of time-ordered perturbation theory is the primary goal of this Chapter. 

Naturally constituent quark models presume the integral in Eq. I|2.8() converges. This 
in turn enables us to relate an initial-state non-wavefunction vertex to the final-state non- 
wavefunction vertex encountered above via integrating around a circle at infinity. Equating 
the sum of residues Res(fc~) -|- Res{k'j^) + Res(A;~) = 0, we find 

£(^.±14) (t- -t-)- ipn^,- _ = _ (2.9) 

r*(A;6-h A,P') T{kb,P) ciJ a c \ J 

This relation holds for the class of models for which the pole structure of Eq. (|2.8j) is not 
modified by the vertices. The ratio structure of the vertices does not allow for a common 
factor k~ — k~ in the two terms in Eq. ()2.9|) . The equality then depends on delicate 
cancellations between initial- and final- state vertices, which in general are unrelated.^ The 
philosophy of constituent quark model phenomenology is to choose the form of and hence 
the form of T. This rules out any such delicate cancellation. Treating T as free to be chosen, 
the equality Eq. 1)2. 9|) can only hold if both ratios are one. This yields the restriction 

^r(A;„P)=0. (2.10) 

But only enters T{kc,P) through k~ = k^ . Hence T is independent of k~ , which is a 
hidden assumption in light-front quark models.^ Clearly this assumption is problematic for 
the calculation of form factors in frames where C 7^ 0- Furthermore, covariance cannot be 
borne in unless the vertex r(A;, R) is completely independent of k^. 

In general, of course, the vertex T not only has light-front energy dependence but poles 
as well. As we shall see below, the light-front reduction of the Bethe-Salpeter equation is 
a procedure for approximating the poles of the vertex function. Moreover when applied to 
current matrix elements, which require further contributions from higher Green's functions, 
these poles generate higher Fock state contributions. 

2.1.2 The reduction scheme 

To reduce the Lippmann-Schwinger equation (or any such Dyson-Schwinger equation) to 
a light-front version, we must introduce an auxiliary Green's function G in place of G (as 

^In particular, one could consider the final state to be an excited state. While the process is not elastic, 
the analysis goes through unchanged and the initial and final state vertices would not permit the delicate 
cancellations required by Eq. 12.91 . 

^There is another way to reveal this hidden assumption. The vertex F cannot be both covariant and 
devoid of k~ poles for the integral in Eq. 12.811 to exist. 
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in [37]). Thus we have 
provided that 



T = W + WGT, 



W = V + V{G - G)W. 



(2.11) 



(2.12) 



Taking residues of Eq. ()2.11|) gives us an alternate way to express the bound state vertex 
function 

r = WGT. (2.13) 

To choose a hght front reduction, G must inherently be related to projection onto the 
initial surface = 0. For simplicity, we denote the integration J ^^{k~\0{R) = 0{R). 
With this notation, we will always work in (1 + l)-dimensional momentum space for which 
the only sensible matrix elements of 0{R) are of the form 1 0{R) |p" , p^) • The operator 

0{R) is defined similarly. The generalization of both operators to (3 + 1) dimensions 
is obvious. For a useful reduction scheme (one that preserves unitarity), we must have 
G{R) = G{R) . The simplest choice of G that results in time-ordered perturbation 
theory requires 

G{R) = G{R) g-\R) G{R), (2.14) 
where the reduced disconnected propagator g{R) is defined by the matrix elements 

{xR+\ G{R) \yR+) = {xR+\g{R)\yR+). (2.15) 

Explicitly this forces 



{xR+\g{R)\yR+) = 2tt6{xR+ - yR+)e[x{l - x)] 



2TTi 



2R+x(l - X 



-Dy^{x\R^). 



(2.16) 



The inverse propagator g (R) can then be constructed, bearing in mind g (R) only exists 
in the subspace where g{R) is non-zero. Explicitly 

{xR+\g-\R)\yR+\ 



9 R+tM — t) 

2tt5{xR+ - yR+)e[x{l - x)]— — ^ '-Dyv~\x\R^). (2.17) 

27rz 



This forces 

{xR+\g-^{R)g{R)\yR+) = 27r6{xR+ - yR+)e[x{l - x)], (2.18) 

which is unity restricted to the subspace where the operators g{R) and g~^{R) are de- 
fined. We are more careful about this point than the authors [16] since the consequences 
of Eq. (|2.17j) are essential for dealing with instantaneous interactions. Notice G{R) defined 
in Eq. ()2.14|1 is consequently non-zero only for plus-momentum fractions between zero and 
one. 

The reduced transition matrix t{R) is 



tiR) = g'\R) G{R)T{R)G{R) g-\R) 



(2.19) 



Taking the residue of Eq. (fTT^ at R^ = M'^ , gives a homogeneous equation for the reduced 
vertex function 7 

\^R)=w{R)g{R)\^R), (2.20) 
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where the reduced auxihary kernel is 



w{R) = g-\R) G{R)W{R)G{R) g-'{R). 



(2.21) 



Given this structure, the reduced kernel w{x^y\R^) = {xR'^\w{R)\yR^) will always be 
oc 9[x{l — x)\6[y{l — y)]. Moreover the reduced vertex 7(x|M^) = {xR'^\^r) oc 9[x{\ — x)] 
as a result of Eq. ^T^. 

From 1)2. 2U() we can define the light-front wavefunction = g{R)\^ji), notice this 

too restricts the momentum fraction x: ^/'(x) oc 9[x(l — x)]. By iterating the Lippmann- 
Schwinger equation for T twice, it is possible to relate T to t and thereby construct T 
given t, which is clearly not possible from the definition (|2.19|) . Taking the residue of this 
relation between T and t yields the reduced-to-covariant conversion between bound-state 
vertex functions, namely 

\Tr) = W{R)GiR)\ \ jr). (2.22) 
Finally, we can manipulate the covariant Bethe-Salpeter amplitude into the form 



1^' 



R) 



1 + iGiR) 



G{R)jW{R) \G{R) 



(2.23) 



which justifies the interpretation of \iPr) as the light-front wavefunction since it is easy to 
demonstrate that {"^r) = IV'i?)- 

While all light-front reduction schemes when summed to all orders yield the x'^ = 
projection of the Bethe-Salpeter equation, the choice of G{R) in Eq. 



(|2.14)) generates a 

kernel calculated in light-front time-ordered perturbation theory. The normalization of the 
covariant and reduced wavefunction will be discussed below. 



2.1.3 An interpretation for the reduction 

The heart of our intuition about the light cone lies in integrating out the minus-momentum 
dependence of the covariant wavefunction. So we merely cast the formal reduction in a way 
which highlights the contributions from poles of the covariant vertex function. 
Utihzing Eqs. (|TT^ and we can show 



GiR)\TR) = G{R) 



hR) 



(2.24) 



Thus the appearance of G{R)\Tr) has the form of an instantaneous approximation since 
Eq. (|2.24|) shows that it has no minus-momentum poles besides those of the propagator 
GiR). 

This instantaneous approximation appears in determining the light-front wavefunction. 
Using Eqs. (|TT^ and (|?^ . we have 

oo 

iVij) = \G{R)ViR)Y,[{GiR) - G{R))v{R)Yg{R)\ \^r). (2.25) 

j=0 

From truncating the series in G{R) — G{R) at some j = n — 1 and using a consistent 
approximation to Eq. (|2.22|) . we are led to the n*'^-order approximate solution ip^"^^ 



1^ 



(Gy)"G'|r), 



(2.26) 
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after having used Eq. H2.24() . Thus at any order n in the formal reduction scheme, we 
have iterated the covariant Bethe-Salpeter equation n-times and subsequently made an 
instantaneous approximation via Eq. (|2.24|) . Retaining the minus-momentum dependence 
in V to n^^ order allows for an n^'^-order approximation to the vertex function's poles. 



2.2 Normalization 



The relation Eq. 1)2. 23() between the four-dimensional Bethe-Salpeter wavefunction and the 
reduced light-front wavefunction was presented in [16]. The normalization of the reduced 
wavefunction, however, was not discussed in much detail and we address this below following 
our work in [38]. 

The reducible four-point function defined by G^"^^ = G + GTG has the behavior 



g^^\r) 



-^^3^ ttt: — + nnite, 



R2 _M'^+ ie 

near the bound state pole. Since the four-point function satisfies the equation 

g^^\r) = g^^\r){g-^{r) - v{r))g^^\r), 

the Bethe-Salpeter amplitude must satisfy 

{^r\(g-\R)-V{R))\^r) 
1 = lim —i 



(2.27) 



(2.28) 



R'^-M'^ 



(2.29) 



which is necessarily finite since |^'_r) satisfies the bound state equation: \^ r) = GV\^ r) 
Application of I'Hopital yields the covariant normalization condition [39] 



2iR^ = {^>r\ 



d 
dR„ 



G-\R)-V{R)]\^R 



(2.30) 



The normalization (|2.3U|) takes the form of a diagonal matrix element of a pseudo current. 

In what follows we shall omit total four-momentum labels since they are all identi- 
cally R. As we saw above, the light-front reduction is performed by integrating out the 
minus-momentum dependence with the help of an auxiliary Green's function G{R). The 
normalization condition for the reduced wavefunction is then deduced by using the conver- 
sion Eq. (|2.23|) and the definition of the reduced wavefunction namely 



I*; 



Hence taking the plus component of Eq. ()2.3U|) . we arrive at 



2iR-' 



(7rI 



G\ 1 + W{G-G) 



dR- 



G- 



V 



i + {g-g)w\g 



hr) 



(2.31) 



(2.32) 



The complicated normalization condition is indicative of the effects of higher Fock space 
components. 
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To see this explicitly, we must know the minus-momentum dependence of the interaction. 
We therefore adopt a weakly coupled, one-boson exchange model for V (the so-called ladder 
approximation). Supposing the boson mass is fi and the coupling constant g, we have 

^(^'P) = 7 .,2^^ 2 ^ ■ • (2-33) 

[p — ky — ^'^ + le 

In considering the normalization condition, we will work in (3 -|- 1) dimensions. The bound 
state equation and wavefunctions which result from the leading-order ladder exchanges are 
collected in Appendix El Notice dV/dR^ = 0. For our purposes, it suffices to assume we 
have a light-front wavefunction ^/^(x,k-'") for the lowest Fock state and then investigate the 
form of the constraint given in Eq. (|2.32|) . 

Let us start with the contribution at leading order in G—G to the reduced wavefunction's 
normalization. 

(7ij| g(-—G-')g hR) = l. (2.34) 



-I 



2R 

To perform the integration, we note 

d 



dR- 



dR- 



G-'{k,R) = -2iR+d-\k){l-x), (2.35) 



where we have customarily chosen x = k~^ /R~^ . Evaluation of the integral in equation 1)2.34(1 
is standard and yields 

a simple overlap of the two-body wavefunction. 

To analyze the normalization to first order in G — G, we expand equation (|2.32|) to first 
order 

N^^ + 6N + ... = 1, (2.37) 

where A^'^^ is the integral appearing in (|2.36|1 and the first-order correction arising from 
is 

SN-^bM \a{^a-%G-G)va + Gv{G-a){^G-^)G\ |™). (2.38) 

The presence of G merely subtracts the leading-order result N^^. Considering for the 
moment just the first term in the above equation (and omitting the subtraction G), we 
have 

"^^^ "^^^ -il-x)Y{x,k^\M^)d{k)d{R-kfV{k,p)d{p)d{R-p)j{y,p^\M^). (2.39) 



(27r)4 (27r) 

The minus-momentum integrals above are similar to those considered in deriving the bound- 
state equation to leading order in Appendix El The only difference is the double pole due 
to the extra propagator d{R — k). With x = k~^ /R'^ and y = p^ /R'^, for x > y we avoid 
picking up the residue at the double pole and the result is the same as in Eq. ()2.34l) after 
using the bound-state equation. This term is then subtracted by the G term in equation 
(|2.38() . On the other hand, when x < y we pick up the residue at the double pole. Part 
of the residue is subtracted by the G term; the other half depends on dV{k,p)/dk~ . The 
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Figure 2.2: Graphical depiction of the irreducible five-point function as sum of impulse 
terms and a gauged interaction. By construction, F^ is gauge invariant. 

second term in Eq. (|2.38l) is evaluated identically up to {k p}. Now combining the two 
terms and their relevant 6 functions, we can rewrite the result using the explicit form of the 
light-front time-ordered one-boson-exchange potential Eq. ID. 31 namely 

(2.40) 

Thus although the covariant derivative's action on the potential vanishes, we can manipulate 
the correction to the normalization into the form of a derivative's action on the light-front 
kernel. With this form, we can compare to the familiar nonvalence probability discussed in 
Appendix El [in the frame where R'^ = R~ = M/^/2, with as the eigenvalue of the 

light-front Hamiltonian, denoted p~ in Eq. HE.8|) ]. This correction to the normalization has 
been discussed earlier [3]. 

Here we have seen that the normalization of the light-cone wavefunction includes effects 
from higher Fock states. Since this normalization condition stemmed from a diagonal ma- 
trix element of a pseudo current, we should not be surprised that matrix elements of the 
electromagnetic current, when treated in this reduction scheme, pick up contributions from 
higher Fock states. These contributions to current matrix elements are the subject of the 
next Section. 

2.3 Current 

Here we extend the formalism presented so far to include current matrix elements between 
bound states. We do so in a gauge invariant fashion following [40]. Our notation, however, 
is more in line with the elegant method of gauging equations presented in [41,42]. This 
latter general method extends to bound systems of more than two particles. Our goal here is 
to see how higher Fock contributions appear from the covariant starting point and rule out 
the use of non-wavefunction vertices. Once one accepts the results, one may work directly 
with the three-dimensional light-front kernel and gauge it. The results are the same [43]. 
Consider first the full four-point function G^^^ defined by 

G'(^) = G + GTG. (2.41) 

For later use, it is important to note that the residue of G^^\R) at the bound state pole 
R^ = is — Using the Lippmann-Schwinger equation for T, we can show the 
four-point function satisfies 

G^^'^ = G + GVG^-^\ (2.42) 

To discuss electromagnetic current matrix elements, we will need the three-point func- 
tion for the i^^ constituent, which we write as d^. We define an irreducible three-point 
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function for the i constituent in the obvious way 

df = d-X^^di. (2.43) 

Now we need to relate the one-particle electromagnetic vertex function to the T matrix. 
Let denote the electromagnetic coupling to the constituent particles (since our particles 

are scalars j^^ = d^). Since the electromagnetic three-point function Tf is irreducible, we 
have 

rf = G-^G(^y (2.44) 
and by using the definition of G^^) Eq. (I^ITI) . we have the desired relation 

rf = f + TGf (2.45) 

Notice the right hand side lacks the particle label i. In the first term, the bare coupling 
acts on the ith particle while in the second term the bare coupling does not act on the ith 
particle. For this reason we have dropped the label which will always be clear from context. 

In considering two propagating particles' interaction with a photon, the above definitions 
lead us to the impulse approximation Fq to the current 

To =r^^2' + rfr'r^- (2-46) 

Additionally the photon could couple to interacting particles. Define a gauged interaction 
topologically by attaching a photon to the kernel in all possible places. This leads us to 
the irreducible electromagnetic vertex F'^ defined as (see Figure I2.2jl 

F^ = F^ + (2.47) 

which is gauge invariant by construction. 

Lastly to calculate matrix elements of the current between bound states it is useful to 
define a reducible five-point function (see Figure ESI) 

(-(5)^ = (^(4)rMG(4)_ (2.48) 

Having laid down the necessary facts about electromagnetic vertex functions and gauge 
invariant currents, we can now specialize to their matrix elements between bound states by 
taking appropriate residues of Eq. (|2.48)1 . The form factor is then 

-iiP'f" + P^)F(t) = (^p/|F'^(-A)|^p), (2.49) 

where P'^' = p/^ + A^" and t = A^. 

2.4 Two models 

In this Section, we use the formalism thus far developed for two models. The first model 
has an instantaneous interaction while the second is not instantaneous. Contributions to 
form factors are contrasted and the issue of non-wavefunction contributions is resolved. 
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Figure 2.3: Graphical depiction of the reducible five-point function G^^^ ^. The irreducible 
five-point function is the gauge invariant T^. 




Figure 2.4: Diagrammatic representation of the one-boson exchange potential w{x,y\M'^) 
appearing in Eq. (|2.52|) . 



2.4.1 Wavefunctions 

To be specific, we again work in the ladder approximation for the kernel Eq. 1)2. 33() for which 
the energy pole (with respect to p) of the interaction V{p, k) is 

This interaction is completely non-local in space-time and hence does not have an instan- 
taneous piece. The reduced kernel Eq. (|2.21j) is consequently made up of retarded terms 
(i.e. dependent on the eigenvalue M^) where higher order in G — G means more particles 
propagating at a given instant of light-cone time (see, e.g., [16]). The leading-order equation 
for ip from Eq. ()2.2U|) is the Weinberg equation [44] 

Jo 2(27r)y(l - y) 

with the time-ordered one-boson exchange potential (see Figure ITHl calculated to leading 
order from Eq. 1)2.21(1 

w(x, ylM'^) = -^^e[x(l - x)]e[y(l - y)] Oix - y)D(x, y\M^) - {x ^ y} , (2.52) 
x-y L J 

where 

2 2 2 

D'^(x, y\M^) = M^ . (2.53) 

y X — y I — X 

We can obtain Eq. H2.51() most simply by iterating the Bethe-Salpeter equation once (see 
Section r2.1.3|) and then projecting onto the light cone. See Appendix El for the (3 + 1)- 
dimensional derivation. 

In the limit /i^ ^ m?,M'^, the interaction becomes approximately instantaneous, which 
suggests we separate out an instantaneous piece Vq'. 

V{k,p) = Vo+[V{k,p)-Vol (2.54) 
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where 



Vo = Vo{x,y) 



-g'^e{x - y) 



+ 



(2.55) 



E'^{x — y) — fi"^ E'^iy — x) — /i^ ' 
with E as a constant parameter to be chosen. Of course other choices of Vo are possible. 
We choose the above form of Vo for two reasons. 

First is the form of the instantaneous approximation wavefunction, which we denote by 
4>{x). When we write the potential as Eq. ()2.54|) . we expand Eq. (|2.12() to first order in 



as 



W = Vo + {V -Vo) + {V - Vo)(G -Gj{V- Vo). 
To zeroth order in G — G and V — Vo, the integral equation for (p is 



(j){x) = -Dyvix\M^) 



^ Wo{x,y)4>{y) 
2(27r)y(l-y) 



dy, 



(2.56) 



(2.57) 



where by Eq. (|2.21j) . the reduced instantaneous potential is merely Wo{x,y) = 9[x{l — 
x)]9[y{l — y)]Vo{x, y). In the instantaneous limit, /i^ ^ m?, solutions to Eqs. (|2.5H) and 
(|2.57|) coincide — both wavefunctions approach ~ Dy^{x,Am?'). 

Secondly, we preserve the contact interaction limit by excluding from Vo the factor 
6[x{l — x)]6[y{l — y)]. That is, when {/i^,(7^} oo with j iJ? fixed (or equivalently 
ii^ ^ 0) we have the contact interaction^ Vo g^ j ij? , which rightly knows nothing about 
the momentum fractions x and y. We shall now investigate contributions to the form factors 
for each of the models Eqs. (|2.,SI-i|) and (|2.55|) . 



2.4.2 Form factors 

From Eq. (|2.49|) we can calculate the form factors for each of the models 1)2. 33p and 1)2. 55|) 
in the reduced formalism. Working in perturbation theory, we separate out contributions 
up to first order by using Eq. (|2.2I'{|) to first order in G — G and F^ (|2.47)) in the first Born 
approximation. The matrix element = ('I'p/|F^(— A)|^'p) then appears for a model with 
some kernel V 



(7P'| GiP')[^l + ViP'){GiP')-GiP')y 
X (a'^(-A)(i2 ^ + V{-A)G{-A)^^'{-A)d2 
x(l + {GiP)-G{P))ViP))G{P) |7P) 



Jlo + ^-^t + ^Jf +^J!^)+ 0\V^' 



with the leading-order result 



7^ 



(7P'| G{P')d^{-A)d^^G{P) 



\1P) 



The first-order terms are 

5J^ = (7P'| G{P')d^{-A)d^^(G{P) - G{P)]V{P)G{P) 



\1P) 



5J^ 

5Ji; 



G{P')V{P') G{P') - G{P') a^(- A)d2 'G(P) 



\1P) 



G(P')n^(-A)G(-A)a^(-A) ^2 G(^) 



\1P) 



(2.58) 



(2.59) 



(2.60) 



^This exactly soluble (1 + l)-dimensional light-cone model, in which ipix) oc Dw{x\M'^), has been con- 
sidered earlier in [45]. 
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The labels indicate the intuition behind the reduction scheme (seen in Section [2.1.3p : the 
term 5J^ arises from one iteration of the covariant Bethe-Salpeter equation for the final- 
state vertex followed by an instantaneous approximation (|2.24|) . 5J^ arises in the same way 
from the initial state and 5J^ comes from one iteration of the Lippmann-Schwinger equation 

At this point, we must specify which component of the current vector we are using to 
calculate the form factor. In Section [2.31 we started with the manifestly Lorentz invariant 
decomposition of the current matrix element, cf. Eq. (|2.49|) . Working to first-order in the 
reduction scheme, however, covariance has been lost, e.g., the kernel Eq. (|2.52|) breaks 
rotational invariance [46]. Thus all components of the current are no longer equivalent. In 
order to make connection with the Drell-Yan formula and the Fock space decomposition 
in [26,27], we must choose the plus-component of the current. Agreement with the form 
factor calculated from the minus-component is only achieved when one works to all orders 
in the reduction scheme. Further discussion can be found in [38] and also below in Section 
13:21 

Because the leading-order expression (|2.59p is independent of the kernel V , the result 
will have the same form for both models. Using the effective resolution of unity, the above 
expression converts into 

W? ^^^P'\P^)^P\G^P')^^^-^)d~,MP)\k) (fe+|7p). (2.61) 

Bearing in mind the delta function present in G{R), we have the factor 

{p\d+{-A)\k) = -i{2k+ + A+){2TTf5'^{p -k-A). (2.62) 

Now define x = k^ /P^ and A+ = —(^P~^ as above in Section f2 . 1 . II and denote the reduced 
vertex {k~^\'yji) = 7(x|M^). The leading-order contribution is then 

j^j*{x'\M^)d{k + A) (2x - C)d{k)d{P - k)j{x\M^), (2.63) 

where x' = jE^- Notice Eq. (|2.63|) is quite similar to Eq. (|2.8() . For x > ( evaluation is 
straightforward and leads to 

for the non- instantaneous case. For the instantaneous case Fio{t)i merely replace tp with 
<t>- 

On the other hand, we know from Section 12.1.11 the region x < ^ contains a non- 
wavefunction vertex (see Figure IT^ . However, evaluation of Eq. (|2.63|) with reduced ver- 
tices is quite different than in Section 12.1.11 Rather simply, the term 'y*{x'\M'^) = by 
virtue of Eq. (|2.2()|) because x' < 0. Thus there is no contribution at leading order for 
X < C- 

The first-order terms depend explicitly on the interaction and will hence be considerably 
different for each of the models. We consider each model separately. 
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Figure 2.5: The Z -graph confronting evaluation of the electromagnetic form factor. 






Figure 2.6: Diagrams which contribute to the form factor to first order in G — G for the 
non- instantaneous case (for x > (). 



Non-instantaneous case 



Evaluating contributions at first order for the non-instantaneous interaction Eq. (|2.33|) is 
complicated by the presence of poles in the interaction [cf Eq. (|2.5U|) ]. First we evaluate 
the first Born term 5J^ in Eq. 1)2. 60() . After careful evaluation of the two minus-momentum 
integrals, we have the contribution to 6J^ for x > ( 



6n 



9{x - C) dxdy {2x - () 



IGn'^xx'y^l — y)y 



^^*{y')D{y\x'\M') 



2.g'^e{y-x) 



y-x 



D{y,x\M^)^P{y), (2.65) 



where y' = This contribution corresponds to diagram A in Figure ITItI Additionally 

using x" = x/C, we have for x < C 



6J 



D 



e{C - x) dxdy {2x - C)/C 



167r2y(l - y)y'x"{l - x") 



//u^9^^(y-a;) 



y-x 



D{y,x\M^)i;{y), (2.66) 



which corresponds to diagram D in Figure 0771 

The initial-state iteration term 6J^ in Eq. H2.6U() is complicated by the subtraction of 




Figure 2.7: The remaining diagrams (characterized hy x < Q for the electromagnetic form 
factor at first order in G — G for the non-instantaneous case. 
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the two-particle reducible contribution 

- (7P'| G{P')d+{-A)d^'G{P)V{P)G{P) 



np) 



(7P'| G{P')d + {-A)d^^G{P) 



np) 



(2.67) 



Thus this term merely removes contributions which can be reduced into the initial-state 
wavefunction. Evaluation of the two minus-momentum integrals yields a contribution for 

X > C to 6J^ 



5J 



B 



fir 9,t'"',l'n ' ^\ rW)D(,'yiM^f^^^D(y, .IM^)^. (2.68) 

Ibir^xx [1 — X )y[l — y) y — x 



which corresponds to diagram B in Figure ITBl On the other hand, for x < we have the 
non-wavefunction vertex 7*(x'|M^) for the final state, which of course vanishes. Similarly 
the G{P) term vanishes. Thus there is no contribution to 5J^ for x < C,. 

Finally there is the final-state iteration term 5J'^ in Eq. (|2.6U|) . There are only two types 
of contributions. For x > ("j that which can be reduced in to the final-state wavefunction is 
subtracted by G. The remaining term is: 



5J+ 



(x-C) dxdy'{2x-Q 



167r^x(l — x)x'y'{l — y') 



'-r{y')D{y',x'\M' 



y -X 



D{y,x\M'^)'il){x), (2.69) 



where y = ^ -|- (1 — Oy' ■ This corresponds to diagram C in Figure ESI For x < C, the 
subtraction term vanishes since x' < for which G{P') = 0. The remaining term in 5J'j 
gives a contribution 



5Ji 



e{C-x) dxdy'{2x-C)/C 
167r2(l - x)x"(l - x")y'{l - y') 



-r{y')D^{x"\t) 



g^e{y-x) 
y - x 



D{y,x\M'^)ip{x), 



(2.70) 

which corresponds to diagram E in Figure 0771 To summarize, the non- valence correction 
to the form factor in the non-instantaneous case is 



6Fni 



1 



l-C/2 

and there are no non-wavefunction terms. 



c 



'D 



(2.71) 



Instantaneous case 

The case of an instantaneous interaction is quite different due to the absence of light-front 
energy poles in Vq. Note we are working with Eq. 1)2. 56() to zeroth order in G — G and 
V — Vq. The first term in Eq. (|2.6fl)) we consider is the Born term 5J^ . The pole structure 
leads only to a contribution for x < 

^■^1 = - y 16vr^y(l - y)y'x"{X - x") ^ ^^^^^ "'^^"^^ ^^'^^^ 

and is depicted by the first diagram in Figure 12.81 The interaction along the way to the 
photon vertex is crossed and represents pair production off the quark line (x" > 1). 
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Figure 2.8: Diagrams with crossed interactions necessary to calculate the form factor in the 
region x < C for the instantaneous case. 



The next term 5J^ simplifies considerably due to the absence of light-cone time depen- 
dence in Vo- As above, the final state vertex restricts C, < x < 1. But then we are confronted 



with a factor (A;| G{P) - G{P) 



since x > 0. Thus 5J, 



0. 



The last term we must consider is ^Jj , in which we have an analogous factor for the 



final state [G[P') - G{P')j \k + A). This is zero for x - C > 0, else G{P') = by virtue 
of Eq. 1)2. 14() and Eq. 1)2. 17() . Thus we only have a contribution for x < C which is from 
G{P')\k + A). The expression for this contribution is 



6J+ 



9iC - x)e{y - x)dxdy'{2x - ()/( 
167r2(l - x)x"(l - x")y'{l - y') 



<P*{y')Vo{y\x')D^{x"\t)<l,{x) 



(2.73) 



and is depicted on the right in Figure 12.81 The interaction again is crossed (x' < 0) and 
represents pair production. With Eqs. (|2.72|) and 1)2. 73|) . we have both bare-coupling pieces 
of the full Born series for the photon vertex (further terms in the series, which result from 
higher-order terms in the expansion of W ^ add interaction blocks to each diagram on the 
quark-antiquark pair's path to annihilation). In this way we recover the full four-point 
Green's function from summing the Born series [29,30]. Notice also the above form of 
is what one would obtain from extending the definition of 7 as a non-wave function 
vertex [29]. For the case of an instantaneous interaction, the light-front Bethe-Salpeter 
formalism automatically incorporates crossing. 

To summarize: we have found the non-valence contribution to the instantaneous model's 
form factor, namely 

5Fi = {5Jt + 54)/{l-Cm. (2.74) 

This expression involves crossed interactions or equivalently terms which could be called 
non-wavefunction contributions. 



Comparison 

Now we compare the form factors for the cases of instantaneous and non-instantaneous 
interactions. In the instantaneous limit /i^ ^> m?,M'^, however, we understand the be- 
havior of both wavefunctions. The wavefunctions Eq. ()2.51|) and Eq. ()2.57|) both become 
narrowly peaked about x = 1/2 in the large limit, cf the behavior of Dw(x,4m2). Since 
we are investigating non-valence contributions to form factors, we choose additionally to 
solve for the instantaneous wavefunction to first order in V — Vo (see Eq. I2.56|) to put the 
wavefunctions on equal footing: i.e. (/)(x) ^ ip{x), because the difference V — Vo is presumed 
small. This has the efficacious consequence of producing identical leading-order terms, cf 
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Table 2.1: Numerical solution of the bound-state equation Eq. (|2.51l) for various values of 
IJ?. The coupling constant a = 0.100. 





4 - M'^/m^ 


0.100 


6.57 X 10-1 


0.316 


2.72 X 10-^ 


1.00 


6.34 X 10"^ 


3.16 


8.79 X 10"^ 


10.0 


9.57 X lO""* 


31.6 


4.81 X 10"^ 



Eq. 1)2. 64() . and eliminates the issue of normalization. Furthermore, the optimal choice of 
the instantaneous interaction is not under investigation here. So we shall simplify the issue 
by choosing £^ = in Eq. (|2.55|) . 

For a few values of /i^, we solve for the wavefunction using Eq. (|2.5H) . In Table we 
list the values of fj? used as well as the corresponding eigenvalue (all for the coupling 
constant a = 0.100, where a = g'^ /Airm'^). We then calculate the form factors in the 
instantaneous and non-instantaneous cases. We arbitrarily choose C = 0.707. In (1 -|- 1) 
dimensions, this fixes A'^ = — (^■^M^/(l — ^). The ratio of these form factors — defined using 
Eqs. ((TMl), (Eini), and (EHH), namely 

{Flo + 6Fi)/{Flo + 6Fm) (2.75) 

— is plotted versus 

logio(/iV"i^) in Figure 1^ The Fi gure indicates the form factors are the 
same in the large /i'^ limit. Since the leading-order contributions are identical, the Figure 
additionally shows the ratio of non-valence contributions to the form factor, namely 

6Fj/6Fni. (2.76) 

This ratio too tends to one as fj? becomes large, of course not as rapidly. Thus for finite /i^, 
we must add the correction V — Vq in Eq. (|2.56j) which results in higher Fock contributions 
present in the non-instantaneous case. 

Lastly we remark that all of this is analytically clear: having eliminated different /i^ 
dependence hidden in ip and (j), expressions for form factors in the instantaneous and non- 
instantaneous cases are identical to leading order in l//i^, i.e. 6J^,6J^,6J^ ~ 1//U^ while 
the remaining non- valence contributions 6J^,5J^ match up with SJ^^ and respectively, 
to l//i2. 

Of course, we have demonstrated the replacement of non-wavefunction vertices only to 
first order in perturbation theory. Schematically, we have required 

\V -Vo\<.\V{G-G)V\. (2.77) 

This is merely the condition that corrections to Vo from finite //^ are larger than second- 
order corrections from the reduction scheme in Eq. 1)2. 56() . Quantitatively this condition 
translates to m/4:TTfi <C a. If this condition is met, the leading corrections are all of the 
form V — Vo and the instantaneous model becomes the ladder model (and non-wavefunction 
vertices disappear when calculating the form factor). When the condition is not met, but 
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Figure 2.9: Comparison of form factors in instantaneous and non-instantaneous cases. The 
ratio Eq. (|2.75|) of form factors (LO + NV) is plotted versus logiQ{^'^ / m^) at fixed C = 0.707. 
Additionally the ratio Eq. 1)2. 76(1 of non-valence contributions to the form factor (NV) is 
plotted. For fj? ^ m?,M'^, the non- instantaneous contributions becomes approximately 
instantaneous. 



holds to \y{G — Gy\^V , say, we must work up to second order in the reduction scheme and 
use the second Born approximation. We do not pursue this lengthy endeavor here, since for 
the non-instantaneous case 7(x|Af^) oc 9[x{\ — x)] necessarily excludes non-wavefunction 
contributions in time-ordered perturbation theory. In fact at any order in time-ordered 
perturbation theory we expect to find a complete expression for the form factor in terms of 
Fock component overlaps devoid of non-wavefunction vertices, crossed interactions, etc. We 
shall make this correspondence explicit when considering the application of the light-front 
BSE to GPDs in Chapter EJ 



2.5 Summary 

Above we investigate current matrix elements in the light-front Bethe-Salpeter formalism. 
First we present the issue of non-wavefunction vertices by taking up the common assump- 
tions of light-front constituent quark models. By calculating the form factor in frames 
where (" 7^ 0, Lorentz invariance mandates contributions from the vertex function's poles. 
Quark models which neglect these residues and postulate a form for the wavefunction are 
assuming not only a pole-free vertex, but a vertex independent of light-front energy. These 
assumptions are very restrictive. 

This leads us formally to investigate instantaneous and non-instantaneous contributions 
to wavefunctions and form factors necessitating the reduction formalism Eq. (|2.23|) . We 
provide an intuitive interpretation for the light-front reduction, namely it is a procedure for 
approximating the poles of the Bethe-Salpeter vertex function. This procedure consists of 
covariant iterations followed by an instantaneous approximation, cf Eq. (|2.26)) . where the 
auxiliary Green's function G enables the instantaneous approximation, see Eq. (|2.24() . In 
order to calculate form factors in the reduction formalism, we construct the gauge invariant 
current Eq. 1)2. 47() . As results in the reduction scheme are only covariant when summed to 
all orders, our expressions violate Lorentz symmetry. To extract the form factor we use the 
plus-component of the current in order to make contact with the Drell-Yan formula and the 
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Fock space representation. 

Using the ladder model H2.33(l and an instantaneous approximation H2.55() we compare 
the calculation of wavefunctions and form factors in the light-front reduction scheme. Cal- 
culation of form factors is dissimilar for the two cases. In the ladder model, which is 
not instantaneous, non-wavefunction vertices are excluded in time-ordered perturbation 
theory. For the instantaneous model, however, contribution from crossed interactions is 
required. Moreover, these instantaneous contributions derived are identical in form to non- 
wavefunction vertices used in constituent quark models. As a crucial check on our results, 
we take the limit fi^ ^ m? , for which the ladder model becomes approximately instan- 
taneous. In this limit, calculation of the form factors for the two models is identical. 

The net result is an explicit proof, for (H-l) dimensions, that non-wavefunction vertices 
are replaced by contributions form higher Fock states if the interactions between particles are 
not instantaneous. We call this a replacement theorem [47]. The analysis relies on general 
features of the pole structure of the Bethe-Salpeter equation in light-front time-ordered 
perturbation theory. Therefore it is trivial to extend the theorem to (3 -|- 1) dimensions. 
When one properly considers the entire coupled tower of Dyson-Schwinger equations, the 
light-front reduction produces light-front versions in light-cone time-ordered perturbation 
theory [34,35]. We will make this connection in our approximation of neglecting the self 
energy when we apply the light-front BSE and its reduction scheme to GPDs in Chapter |31 
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Chapter 3 

Generalized Parton Distributions 



The light-front BSE development and resolution of the non-wavefunction controversy dis- 
cussed above have natural applications. In Chapter |21 we worked through matrix elements 
of the electromagnetic current in a frame where A"*" 7^ 0; we can now make the connection 
to generalized parton distributions (GPDs). These distributions, which in some sense are 
the natural interpolating functions between form factors and quark distribution functions, 
turn up in a variety of hard exclusive processes, e.g., deeply virtual Compton scattering 
(DVCS) and hard electro-production of mesons [48-52]. There are a number of insightful 
review articles on this subject from a variety of perspectives [12,53-56]. The scattering 
amplitude for these processes factorizes into a convolution of a hard part (calculable from 
perturbative QCD) and a soft part which the GPDs encode. We will use DVCS to motivate 
the introduction of GPDs. Since light-cone correlations are probed in these hard processes, 
the soft physics has a simple interpretation and expression in terms of light-front wavefunc- 
tions of the initial and final states [26,27]. After analyzing DVCS, we will see this explicitly 
for the Wick-Cutkosky model. 

Below, we take the light-front BSE model of Chapter [21 into (3-1-1) dimensions and make 
clear the connection to higher Fock components by explicitly constructing the three-body 
wavefunction from our expressions for form factors. The immediate application of our de- 
velopment for form factors (in a frame where the plus-component of the momentum transfer 
is non-zero) is to compute GPDs. Connection is made to the Fock space representation and 
continuity of the distributions is put under scrutiny. We show how continuity is ensured by 
light-cone time-ordered perturbation theory for the Wick-Cutkosky model. The light-front 
Bethe-Salpeter formalism can also be employed to obtain time-like form factors. This latter 
application is interesting since no Fock space expansion in terms of bound states is possible. 
In the interest of organization, a summary of our investigation in [38] for timelike processes, 
appears in Appendix IfI 

The Chapter is organized as follows. First we review the basics of GPDs in Section EI 
including the field theoretic properties they must obey. The introduction of GPDs is framed 
in the context of DVCS. This discussion of DVCS is the natural extension of our review of 
deep-inelastic scattering in Section [TTTl Next in Section we use the light-front reduction 
scheme to calculate generalized parton distributions in the Wick-Cutkosky model. We do so 
by using the integrand of the electromagnetic form factor calculated in an arbitrary frame 
[these expressions in (3-1-1) dimensions are collected in AppendixEl where the leading-order 
bound-state equation for the wavefunction also appears] . We discuss the continuity of these 
distributions in terms of relations between Fock components at vanishing plus momentum. 
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Figure 3.1: Diagrammatic depiction of the DVCS process. An electron of momentum k is 
incident on a proton of momentum P. The leading-order electromagnetic process is virtual 
photon exchange between the electron and proton. As a result, the electron's momentum is 
altered to k' . The reaction produces an intact proton of momentum P' and a real photon, 
g'2 = 0. 

Connection is also made to the overlap representation of generalized parton distributions. 
Construction of the higher Fock components directly from light-front perturbation theory 
is reviewed in Appendix El Further applications of the reduction formalism to timelike 
processes are contained in Appendix El We conclude our study of the light-front BSE with 
a brief summary (Section 13.2.3)1 . Finally we discuss limitations in modeling GPDs from a 
few different perspectives. 

3.1 Definitions and properties 

In this Section we define GPDs in terms of a light-like correlation of hadronic matrix 
elements. In order to motivate this definition, we consider DVCS which is the theoretically 
cleanest process into which the GPDs enter. The discussion of DVCS below builds on 
that of DIS presented in Section 11.11 In contrast to DIS, DVCS is an exclusive process 
where the final states are detected. Despite this difference, DVCS and DIS are quite similar 
theoretically. 

3.1.1 DVCS 

In DVCS off the proton (or other targets), one is interested lepton scattering at high mo- 
mentum transfer. The basic reaction is (see Figure l3T|) 

eik) + p{P) e{k') + p{P') + 7(g')- (3.1) 

Thus the situation of the reaction is similar to DIS in setup, however, afterward one is 
interested in measuring not only the final electron, but also recovering an intact proton 
(altered to some momentum P') in addition to a real photon, q'^ = 0. Experimentally 
it is easier to deduce the existence of the final state photon by missing mass techniques 
rather than to detect it directly. Furthermore the existence of a plethora of other channels 
requires high luminosity beams and long run times in order to filter out enough DVCS 
events exclusively. 

As in DIS, we label the momentum transfered to the proton as q = k — k' . We will again 
work to leading order in QED, for which the momentum is transfered via a single virtual 
photon exchange. In the deeply virtual limit, the photon's virtuality becomes arbitrarily 
large, q'^ — oo. The amplitude for this process can be written in terms of the off-diagonal 
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Compton amplitude T^'^(g, P, A), where A = P' — P is the momentum transferee! to the 
final state proton. The off-diagonal Compton amplitude is 

r^'^(g,P,A) = I d^ye''i-y{P'\T{j^^{0)r{y)}\P), (3.2) 

and differs from the DIS case only by the final state proton being off the momentum-space 
diagonal. The overall amplitude A for the DVCS process has the form 

^DVCS = e;{q')L,{k,q)T^'''{q,P,A), (3.3) 

where e*(q'') is the polarization vector for the final state photon, and L^{k,q) describes 
the QED virtual photon exchange off the electron line. Notice q' = q — A hj momentum 
conservation. We include the form of the amplitude Eq. H3.3|l in order to contrast with 
the case of DIS, where the amplitude squared was under investigation. Moreover we are 
interested in the amplitude as a whole, not just the imaginary part. 

Our setup so far has paralleled the case of DIS. We shall thus choose our variables in a 
similar fashion. This choice corresponds to the asymmetric frame employed by Radyushkin. 
Nearing the end of this section we will also provide expressions in the symmetric frame of Ji. 
In the rest of this work we will make use of both frames and their correspondingly different 
sets of variables. In the asymmetric frame, we have P"*" = 0, as in Section 11.11 and the 
decomposition of the initial state's momentum P^ is identical to Eq. H1.4() . Let us define 
the Bjorken variable for DVCS as 

Further let us choose g • n = 0, so that —q^ = gives the photon's virtuality. This fixes 
the photon momentum to be 

Now for the t-channel momentum transfer A'^, we have A^ = t. Furthermore A is restricted 
from the constraint q'^ = from the final photon. Thus we have the decomposition 

A/^ = -Cp'^ + A^ + ^^^n'^, (3.6) 

which holds only as q^ — oo. Notice the plus component of the final state photon 
momentum is q'^ = C-P^- 

As we know from DIS, in the deeply virtual limit the Compton amplitude is dominated 
by the leading light-cone singularity in the product of currents. Taking this contribution 
generates the so-called handbag dominance due to asymptotic freedom, see Figure l3^ The 
only differences compared to Section Fl. II are that the final states have altered momenta, the 
final photon is real, and we are not just looking at the cuts of the diagrams. There is an 
additional requirement so that the amplitude maintains a clear separation of scales, namely 
the stipulation that \t/q'^\ <^ 1. 

Using the handbag dominance, we arrive at the leading contribution to the off-diagonal 
Compton amplitude in the deeply virtual limit 



^ {k + qf + ie^ ^ (fc-g + A)2 + ie' 



M(fc,P,A)|, 
(3.7) 
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Figure 3.2: Graphical depiction of the handbag mechanism in DVCS. There are two con- 
tributions to the off-diagonal Compton amplitude: the handbag shown on the left, and the 
crossed bag shown on the right. Figure adapted from [12] 



where we have written down contributions from both handbags. The off-diagonal free 
quark-proton Green's function appearing in Eq. 1)3. 7j) is 



P) 



(3.8) 



which again differs from the DIS case only by the momentum A transfered to the final state. 

We now use the collinear approximation to obtain the leading contribution from this 
Green's function in the deeply virtual limit. The leading contribution in Eq. ()3.7|) can be 
isolated by taking k, (k ■ n)p^. Further we shall define X = k ■ n, i.e., X = k'^ /P^. In 
the collinear approximation, the off-forward Compton amplitude now appears as^ 



T^^''{q,P,A) = ^{g^^''-pV 



dX 



1 



+ 



1 



X-Q + ie X + ieJ 
where we have introduced the function M{X, t) defined by 

(fk 



)M(X,C,t), (3.9) 



M^p{X,C,t) 



(27r)' 



-^5{X -k-n)M^p{k,P,A). 



(3.10) 



As in the DIS, the integral over y that is buried in Map{k, P, A) can be evaluated trivially 
and we thus arrive at 



M^f,{X,C,t) 



dX 
2^ 



(P'|V^(0)V'a(An)|P) 



(3.11) 



Following Section Fl.H we now separate Mq^(X, ^, t) into irreducible structure functions. 
We will do this only schematically. First we write out the general dependence on Dirac 
structures 



1 



1 



M^^iX, C, t) = -i>^p Fix, C, t) + -Kb FsiX, C, t) + -fi^^ M'F^{X, C, t) 



(3.12) 



^As in DIS, the are perturbative QCD corrections to the formula Eq. 13.911 that resuh in the factorization 
scale dependence of both the hard scattering kernel and the soft matrix element, here M{x, (,t, fj,^), such 
that T'^'^ {q, P, A) is scale independent. Additionally similar to DIS, we can interpret fj, as the resolution 
scale for the partonic process. For a complete discussion of 0(as) corrections to DVCS and beyond, see, 
e.g. [12,56]. 
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Next we observe that by taking various traces we can write 

/ ^ e*""" (P'|V^(0)7'^V'(An)|P) = 2 [F{X, (, t) + F^{X, (, t) A^^ + U'^F^iX, C, t) n^] . 

(3.13) 

Finally we identify F{X, t) as the leading structure function for DVCS.'^ The contributions 
from the other structure functions F^{X^ (, t) and F4{X, t) are of higher twist (suppressed 
by |A_l|/(5 and M^/Q^, respectively). The GPDs are contained in the function F(X, ^,t). 



3.1.2 GPDs 

In the last Section, we reduced the amplitude for DVCS down to one leading contribution 
in the deeply virtual limit. This leading contribution is a convolution of a hard scattering 
part [seen in Eq. (|3.9j) ]. and a soft hadronic matrix element F{X,C,,t). The latter contains 
the GPDs. In this Section, we will inspect the GPDs more closely and we shall also switch 
frames from Radyushkin's asymmetric frame to Ji's symmetric one. The switch to the 
symmetric frame illuminates further properties of the GPDs. 

Utilizing Eq. (|3.13j) . we can write the leading- twist structure function as 

F{X,C,t) = \ I ^e^^^(P'|V^(0)riV(An)|P). (3.14) 

From this relation we see immediately that the PDFs are contained in the limit A'^ — > 0, 
namely F{X, 0, 0) = fi{X). For the case of the proton (or any other spin-i target), we can 
decompose the matrix element in terms of two independent GPDs H{X, t) and E(X, C, t), 
viz. 

^H{X,C,t) + ''"'y.f'' EiXX,t) 



F{X,C,t) = u{P') 



2M 



u{P). (3.15) 



Notice by virtue of the Gordon identity (see Appendix^), a possible third structure propor- 
tional to n • A is not independent. Furthermore, the reduction relation to PDFs involves 
only H{XXit), i-e., fi{X) = H{X,0,0) The independent function E{X,(,t) is simply not 
visible in DIS because its presence in any amplitude requires momentum transfer. 

The decomposition of the function F{X,(,t) into GPDs H{XX,t) and E{X,(,t) in 
Eq. 1)3. 15(1 is reminiscent of the Dirac and Pauli electromagnetic form factors. Moreover, 
in DVCS the proton remains intact which suggests an even deeper connection with the 
proton's elastic form factors. This connection is afforded by the "sum rule" that produces 
the local limit of the bilocal field operators in Eq. p.l4|) . Notice the effect of integration 
over X in Eq. ()3.14|) produces the local limit 



e^^^^(0)/i^(An) =5(A)V'(0)^V(0). (3.16) 

2tt 

Thus we find the sum rules 

J dXH{XX,t)=Fi{t) (3.17) 
dX E{XX,t) = F2{t), (3.18) 



/ 



■^As in DIS, beyond leading order in QCD perturbation theory, the bilocal operator in Eq. IjH.l^^ must 
renormalized. This re: 
the twist decomposition. 



be renormalized. This renormalization leads to fji^ dependence of the structure functions but does not affect 
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where Fi{t) is the Dh'ac form factor and F2{t) is the Pauli form factor. The ( dependence 
drops out^ upon integration over all X. This is a manifestation of Lorentz invariance 
because in calculating form factors, one can choose a frame in which the momentum transfer 
A+ = — C-P^ = 0. For DVCS, however, there is no remaining freedom to use such a frame 
since the virtual photon's momentum was chosen so that = 0. The requirement C 7^ 
implies that DVCS is sensitive to physics otherwise absent from electromagnetic form factors 
due to Lorentz invariance. 

In fact, there are more stringent constraints imposed by Lorentz invariance. The n^^ 
moment of each GPD with respect to X must be a polynomial in <^ of order n. This 
constraint is referred to as the "polynomiality" property of GPDs. We will devote ample 
space to the issue of polynomiality and the special case of sum rules in Chapters |1] and |S1 

There is a further property that the GPDs must satisfy: they are bounded by the 
positivity of the norm on Hilbert space. This condition implies a bound, referred to generally 
as a positivity bound or very imprecisely as just the positivity property of GPDs. Looking 
at the field operators in Eq. (|3.14|) . we see that the definition of F{X, t) is similar to that 
of a mixed density. As such it is constrained by the PDFs. Using the Cauchy- Schwartz 
inequality, one finds the bound 




9iX-C)F{X,C,t) <y/i(^Xj/i(^-|) (3.19) 

that we shall refer to as positivity.^ The positivity bound is a property that shall haunt us 
throughout Chapters and [3 and thus we only briefly mention it here. 

As a final point in our introductory discussion of GPDs, we switch reference frames 
and variables to define GPDs that are more symmetric with respect to the initial and final 
states. To do so, we define the average momentum between the states as = ^(P + P')^ 
so that the initial momentum P = P — A/2 and the final momentum P' = P+ A/2. Now we 
choose a frame of reference in which P-*- = 0. This frame is related to the asymmetric frame 
by a boost and so strictly speaking A_|_ is different from above. It is then more convenient 
to use a modified Bjorken variable. Recall ^ = — A+/P^ is a measure of the longitudinal 
momentum transfer relative to the initial proton's momentum. In the symmetric frame, a 
useful quantity to consider is the amount of longitudinal momentum transfer relative to the 
average proton momentum. Thus we define the variable ^ by 

A • n A+ , , 

C = — = — = — . (3.20) 
2P-n 2P+ 



' The local limit also produces an operator of zero anomalous dimension and hence the sum rules Eq. ijlHTjl 
and 13.1811 are independent of the renormalization scale fi. Physically the form factors do not depend upon 
the resolution scale chosen because it is the (transverse) distribution of total charge that is probed. 

*The positivity bound is stable under QCD evolution, i.e. with Eq. 13.191 modified to include fi^ de- 
pendence in all three functions. This can be demonstrated by studying the evolution kernels. Alternately 
one can use the light-cone wavefunction representation where the Fock components themselves acquire scale 
dependence (because gluon radiation mixes Fock components). Since the positivity bound is manifest in the 
wavefunction representation, writing down dependence in the light-cone wavefunctions in F{X,(^,t, n^), 
one immediately sees 



e{X-0\F{XX,t,iJ.')\ < Jf,{x,t^^)f^[^^,tI^ 
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The conversion between C and ^ is 




(3.21) 



The properly symmetric definition for the GPDs appears as 




1 /" ^ i\x 
2. 27r 



{P'\ip{-Xn/2)fiiP{Xn/2)\P) 



(3.22) 



which features the variable x instead of X 



The two are related by 



X 



1 + e 



(3.23) 



because x is measured relative to the average plus momentum. The operator combination 
in Eq. (|3.22|) is Hermitian and so -F(x,^,i) is real. Under the transformation A'^ — > — A^, 
we thus have F{x, — ^, t) = F{x, ^, t) which is commonly referred to as .^-symmetry. This 
symmetry naturally extends to the GPD moments which are hence required to be even 
polynomials in ^. 

As was the case for PDFs, a natural physical picture for the GPDs emerges when we 
use light-cone quantized fields. This is because the bifocal operator entering into the GPD 
definition is an equal light-cone time operator and not an equal time operator. Despite this 
similarity to PDFs, there is not a general probabilistic interpretation for the GPDs because 
the requisite matrix element has more the form of a mixed density operator. Decomposition 
of the GPD in terms of light-cone quark and anti-quark modes is presented in Appendix lO 
In the special case ^ = 0, one can manipulate the GPDs into true probability distributions 
of quarks with momentum fraction x in impact parameter space [57, 58]. The case for 
^ 7^ does not have a probabilistic interpretation, however, some physical insight can be 
gained [59]. 

Having motivated the definitions of GPDs from analyzing DVCS and having presented 
a brief description of their properties, it makes sense to consider a toy model example of 
GPDs. We do so in the following Section. 



In this Section, we cast our results for form factors in the Wick-Cutkosky model [38] in 
the language of GPDs and the light-cone Fock space expansion. The (3 + l)-dimensional 
expressions for form factors are presented in Appendix ^ These results are the higher 
dimensional generalizations of the (1 + l)-dimensional expressions we derived in Section 
12.41 Additionally one can obtain these results directly from time-ordered perturbation 
theory using two-body projection operators as explicated in Appendix IeI 

In contrast to the proton case (to which we return in Chapter ISJ, the GPD^ F{x,Cj't) 
for our meson model is defined by a non-diagonal matrix element of bilocal field operators 



■''Notice we have readily abused notation. Here we use the variable x to denote the active quark's plus 
momentum fraction in the asymmetric frame. Above in Section 13.11 we used X but now return to the 
Wick-Cutkosky model of Chapter |21 where we employed x throughout. 



3.2 GPDs and the light-front BSE 




(3.24) 
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where (/'(x) denotes the quark field operator and d'^ = — is the electromagnetic 
coupHng to scalar particles. Here we are considering only the twist-two contribution^ and 
have taken accordingly the plus-component of the vector current in Eq. (|3.24|) . Comparing to 
the current matrix element in AppendixEl the definition of the GPD leads immediately 
to the sum rule 

-F{x,Ct)=Fit), (3.25) 



l-C/2 

where F(t) is the electromagnetic form factor for our model. Hence one can calculate 
these distributions from the integrand of the form factor because the non-local light-cone 
operator merely fixes the plus momentum of the active quark. In this way, the light-cone 
correlation defined in Eq. (|3.24)) has a natural description in terms of light-front time- 
ordered perturbation theory, e.g., for x > ( the relevant graphs contributing to the GPD 
are in Figure 12. 6( and those for x < ( are in Figure 12.71 Conversion of the contributions 
to the form factor into GPDs is straightforward using Eq. p.25p : we merely remove —2iP~^ 
and J dx from Eqs. HD.13ID.18|) . These cumbersome expressions are not worth generating 
here. 



3.2.1 Continuity 

In order for the deeply virtual Compton scattering amplitude to factorize into hard and soft 
pieces (at leading twist), the GPDs F(x,C,i) must be continuous ai x = C,, see Eq. (|3.9|) . 
There is no proof of this from QCD, but the power-law scaling seen experimentally allows 
us to infer the GPDs continuity. Maintaining continuity at the crossover is more pressing 
because experiments which measure the beam-spin asymmetry are limited to the crossover 
[60]. Let us scrutinize continuity in our results for the Wick-Cutkosky model. 

The leading-order expressions we have derived for the Wick-Cutkosky model are con- 
tinuous. This is easy to see since the contribution for x < ^ is identically zero. The 
valence wavefunction contribution for x > C is a convolution of wavefunctions one of which 
is ■0*(x', . . .) which is probed at the end point since x' = — > 0. From the bound-state 
equation Eq. (|D.5|) . we see the two-body wavefunction vanishes quadratically at the end 
points. Taking into account the overall weight x'~^, the valence piece vanishes linearly at 
the crossover. At leading order then, continuity is maintained at the crossover, while the 
derivative is discontinuous. Working only in the valence sector, valence quark models will 
never be of any use to beam-spin asymmetry measurements since the value at the crossover 
requires one wavefunction to be at an end point. ^ 

Let us now check the next-to-leading order contributions to the GPD for continuity. 
First we shall deal with the term stemming from iterating the Bethe-Salpeter equation for 
the initial state HD.16|) (see diagram B of Figure 12.6(1 . Since there is no corresponding 
Z-graph generated from iterating the initial state, we expect this contribution to vanish. 
Looking at the expression, we see again the factor ip*{x' , . . .)/x' which vanishes linearly as 
X ^ Q. Moreover there are the interaction terms: first D(y, p"*"; x, k-'"|M^) which is finite 
as X ^ C, and then 

I^^V^k'^lM^) (3-26) 

® Additionally by ignoring self-energy corrections, there is no dependence in the model. 
In the three-body bound state problem (e.g. the nucleon), the valence GPD will vanish only if the 
three-body interaction is non-singular at the end points (which is physically reasonable and perturbatively 
true) [61]. 
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which vanishes at the crossover. Thus not only does the initial-state iteration term vanish 
at the crossover, its derivative does so as well. 

Now we investigate the Born terms Eq. (|D.14|) (see diagram A of Figure 12. 6|) and 
Eq. HD.15|) (see diagram D of Figure 12.7(1 at the crossover. Approaching from above, 
we use Eq. ()D.14|1 . and we have the finite contribution at the crossover 



9^0iy-C)/iy-C) 



(k^ + A^)2 + m2 



Z)(y,p^;C,k^|M2)V;(y,p^). (3.27) 



On the other hand, approaching the crossover from below in Eq. HD.15|) we have to deal 
with singularities as x" = x/C ^ 1. Writing out the propagator for the quark-antiquark 
pair heading off to annihilation, we see 

This linear vanishing cancels the weight (1 — x")~^. Taking the limit x — > then produces 
equation H3.27|) and thus the Born terms are continuous. 

Lastly we must see how the final-state iteration terms match up at the crossover. Using 
Eq. (|D.17|) to approach (" from above (see diagram C of Figure 12. 6|) . we have the contribution 



(167r3)2(l - cy(l - y' 



^'/^ r^(2/,P^;C,k^|M2)V(C,k^). (3.29) 



(k-L + A-L)2 + m2' 



Approaching from below (see diagram E of Figure 12.7(1 . we utilize equation ()3.28() in 
taking the limit of ((D.18() . The result is (|3.29() and hence we have demonstrated continuity 
to first order, i.e. 

F{C, C, t) = F(C, C, tf°''' + F{C, C, tf^^' (3.30) 
no matter how we approach x = 

3.2.2 Fock space representation 

We now write the derived GPDs in terms of Fock component overlaps. In the diagonal 
overlap region x > ( this will be a mere rewriting of our results, while there is a subtlety 
for the non-diagonal overlaps. To handle the zeroth-order term in Eq. ((D.13() . we define the 
two-body Fock component as 

i^2{xi,'ki,X2,'k2) = 1. — V(a:i,k4i), (3.31) 

^/XlX2 

noting that the relative transverse momentum can be defined as k^j = X2k]'- — xik.2 ■ In 
terms of Eq. ((3.31() . the GPD appears as 

F(x,C,t)^° = ^-^j^ I [dx]2[d\.% Kx-x,W2{<M^)^^Mx^Ml (3-32) 
V S J -I o , x',x^ 



j=i,2 y 
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where the primed variables are given by 



. ' , , for i / j and <^ ^ , , (3.33) 

and the integration measure is given by 

N N 



dx]N = Wdxi5(l-^Xi\ (3.34) 



i=l 1=1 

N N 

t'^'^^]^ = [2(2^)3]^-i n < E ^0 • (3-35) 

Notice the sum over transverse momenta in the delta function is zero since our initial meson 
has P"*" = 0. The sum over j in Eq. (|3.32|) produces the overall factor of two for our case 
of equally massive (equally charged) constituents. 

To cast the next-to-leading order expressions for x > ( in terms of diagonal Fock space 
overlaps, we must write out the three-body Fock component. Looking at the diagrams in 
Figure 1221 it is constructed from the two-body wavefunction in the following manner 



: 2(27r)3 f . Myj,Pf) 

^/xiX2X3 J y/ym 



r2\ 



e{yi - xi)x3S{y2 - X3)6{pt -kt)D{yi,pt;xi,kt\M' 

+ e{y2 - X3)xi6{yi - xi)5ipi - ki)D{y2,pt,X3,H\^^) 



(3.36) 



where i runs from one to three and the label j, which stems from the integration measure, 
runs from one to two. We discuss how to obtain this three-body wavefunction directly from 
time-ordered perturbation theory in Appendix El Using ip^ in Eq. (|3.36|) . the terms in the 
GPD at first order in the weak coupling can then be written compactly as 

F(x,C,t)^^° = ^^^f;^^ [[dxUdk^js Six-x,)r^{x',M^)^^M^^,'i^i)■ (3-37) 



j=l,3 \j^3^3 

One can verify that the diagrams in Figure 1221 are generated by Eq. (|3.37)) . Additionally 
there is a fourth diagram generated by Eq. 1)3. 37() which does not appear in the figure. 
This missing diagram is characterized by the spectator quark's one-loop self interaction 
and is absent since we have ignored f{k'^) and the scale dependence of light-cone wave 
functions.'^ The absence of this diagram does not affect continuity at the crossover. The 
missing diagram vanishes at x = C since the final-state wavefunction is tl^* {x' , . . .) , which 
is probed at the end point. Lastly we note the above Fock component overlaps satisfy the 
positivity constraint for a composite scalar composed of scalar constituents [63]. Moreover, 
with the inclusion of the self-energy loop, the form factor is gauge invariant [43] and hence 

^Recently the scale dependence of light-cone Fock components in QCD has been investigated [62]. In the 
Wick-Cutkosky model, the scale dependence is of course different and leads to different evolution compared 
to perturbative QCD. 
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the GPDs are properly gauged, in the sense that the bilocal Ught-cone operator is inserted 
in all possible places. 

Now we must come to terms with the non-diagonal overlap region, x < At first 
order, the diagrams of Figure ITTl correspond to four-to-two Fock component overlaps. We 
have been cavalier about time ordering, however. The expressions Eqs. HD.15() and HD.18() 
do not correspond to time-ordered graphs. Both terms contain a product of time-ordered 
propagators: one for the two quarks leading to the final-state vertex and another for the 
quark-antiquark pair heading off to annihilation. But for an interpretation in terms of a 
four-body wavefunction, all four particles must propagate at the same time. This is a subtle 
issue as a graph containing the product of two independently time-ordered pieces (where 
one leads to a bound-state vertex) corresponds to a sum of infinitely many time-ordered 
graphs. It is easiest to write out the relevant pieces in terms of the propagators' poles. The 
quark, anti-quark heading to annihilation have propagators d{k) and d{k + A) and poles we 
label k~ and k~ , respectively. The remaining propagators of interest d(p') and d(P' — p') 
have poles p~ and p^ , respectively. We can then manipulate as follows 

1 1 _ 1 ( ^ ^ 

kc ka Pfj Pa Pfy Pa ~\~ kc ka Pa kc — ka 

In this form, we have produced the correct energy denominator for the instant of light-front 
time where four particles are propagating. Multiplying this denominator by the three-body 
wavefunction yields the four-body wavefunction (up to constants). This is the part of the 
four-body wavefunction relevant for GPDs (there are additional pieces for two-quark, two- 
boson states, see Appendix 0. In the resulting sum Eq. ()3.38|) . the first term will produce 
the two-body wavefunction for the final state and we have a genuine four-to-two overlap. 
We do not write this out explicitly. 

The second term in Eq. 1)3. 38() . however, contains again the propagator for the pair 
heading to annihilation. Using the light-front Bethe-Salpeter equation for the vertex (which 
contains infinitely many times) we can introduce a factor of the time-ordered interaction. 
The resulting product of independent time orderings can again be manipulated as in Eq. 
1)3. 38() . The result produces another overall four-body denominator which contributes to the 
four-body Fock component of the initial state. Since we iterated the interaction, however, 
this new contribution is no longer at leading order and can be neglected. Thus the second 
term in (|3.38|) does not contribute at this order. 

Having manipulated the GPDs into non-diagonal overlaps for x < C, we must wonder if 
continuity at the crossover is still maintained. In the limit x — > the light-front energy of the 
struck quark goes to infinity. Consequently k~ , which contains this on shell energy, is infinite 
and dominates the four-body energy denominator. This is identical to the reasoning in 
Eqs. (|D.15|) and (|D.18|) where instead of the four-body denominator, we have the propagator 
D^^{x" + x"A.-^\t), which is dominated by k~ at the crossover. Either way, we arrive 
at the expressions found above for the crossover Eqs. (|3.27|1 and (|3.29|) . 

Having cast our expressions for generalized parton distributions in terms of the Fock 
components, we can enlarge our understanding of the sum rule and continuity at the 
crossover. Both must deal with the relation between higher Fock components. The way 
^-dependence disappears from Eq. (|3.25|) mandates a relation between the diagonal and 
non-diagonal Fock component overlaps that make up the GPD. Such a relation between 
Fock components must follow from the field-theoretic equations of motion. Continuity itself 
is a special case of the relation between Fock components, specifically at the end points. 
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Above we have seen our expressions are continuous (and non-vanishing) at the crossover 
and exphcitly that the three- and four-body components match at the end point (where 
X — = 0). This weak binding model for behavior at the crossover is a simple example of 
the relations between Fock components at the end points^. More general relations must be 
permitted from the equations of motion to guarantee Lorentz covariance (e.g. in the struc- 
ture of the Mellin moments of the GPDs, of which the sum rule is a special case). Here, of 
course, Lorentz symmetry is violated. Infinitely many light-cone time-ordered graphs are 
needed in the reduced kernel to reproduce the covariant one-boson exchange ()2.33|1 . Thus 
exactly satisfying polynomiality requires not only infinitely many exchanges in the kernel 
but contributions from infinitely many Fock components. It should be possible, however, 
to show how the sum rule and polynomiality are improved order by order. 

3.2.3 Summary 

Above we have investigated the application of the light-front reduction of current matrix 
elements to GPDs. In Appendix El we reviewed the derivation of the form factor at next- 
to-leading order in the (3 -|- l)-dimensional ladder model. These expressions were then 
converted into the GPD for the model. Continuity of these distributions at the crossover 
(where the plus momentum of the struck quark is equal to the plus component of the 
momentum transfer) was explicitly demonstrated, cf Eq. (|3.3fl|) . Connection was made to 
the overlap representation of GPDs by constructing the three-body wavefunction to leading 
order in perturbation theory. As a check on our results, we also reviewed the construction of 
higher Fock states from the valence sector in old-fashioned time-ordered perturbation theory 
(Appendix El • The derived overlaps Eqs. (|3.32|) and 1)3. 37(1 satisfy the relevant positivity 
constraint. The non-vanishing of the GPDs at the crossover could then be tied to higher 
Fock components, specifically at vanishing plus momentum, and are hence essential for any 
phenomenological modeling of these distributions. This rewriting allowed us to understand 
how continuity arises perturbatively from the small-x behavior of Fock state wavefunctions. 
In perturbation theory, the diagonal valence overlap vanishes at the crossover, while the 
higher Fock component overlaps do not. In general the n-to-n overlap matches up with the 
(n + l)-to-(n — 1) overlap at the crossover due to the dominance of the rebounding quark's 
infinite energy. In essence, no matter how complicated the higher Fock component is, the 
events in the last instant of light-cone time must simplify. Either the active quark does 
not interact in this time, for which the contribution vanishes, or else there is an interaction 
which then reduces to a contribution of a Fock component with two fewer particles. This 
is true to all orders in perturbation theory and holds for bound states of any number of 
particles due to the nature of the time-ordered kernel and light-cone perturbation theory. 

Unfortunately issues involving Lorentz invariance (such as the sum rule for the electro- 
magnetic form factor and the polynomiality constraints) are left untouched. To maintain 
covariance one would need infinitely many time-ordered exchanges in the kernel as well as 
infinitely many Fock components. It should be possible, however, to understand pertur- 
batively how the C-dependence disappears from Eq. (|3.25|) . This requires further relations 
between Fock components and these should be afforded by the field-theoretic equations of 
motion. An additional application of the reduction formalism for currents is to calculate the 
time-like form factor. As seen in Appendix El and in [38], we calculate the ladder model's 
generalized distribution amplitude, which is related to the time-like form factor via a sum 

'The known Fock component end-point relations in gauge theory derived in [64] do not shed hght on the 
continuity of GPDs. 
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rule similar to Eq. ()3.25|) . This is in contrast to the non-existent Fock space expansion for 
these types of processes. Thus the reduction formalism succeeds in generating contributions 
where the light-cone Fock component wavefunctions cannot be employed. 

With the formalism explored here, one could use phenomenological Lagrangian based 
models to explore both generalized parton distributions and generalized distributions am- 
plitudes within the light-front framework. Such an investigation is interesting not only for 
testing phenomenological models, but also for anticipating problems for approximate non- 
perturbative solutions for the light-cone Fock states. Nonetheless more model studies are 
warranted before truly realistic calculations can be pursued. In particular, in light of our 
discussion continuity of GPDs may be transparent in light-cone time-ordered perturbation 
theory, but one must go beyond the toy model considered above and include spin degrees 
of freedom. 

3.3 Limitations 

In order to assess future data on GPDs one needs phenomenological parameterizations for 
the new distributions. The goal of the rest of this work is construct one. For ordinary 
quark distributions, one has a variety of forms invented by hand, as it were. For GPDs 
however, the constraints do not permit such cavalier parameterizations. Additionally the 
parametrization should be simple enough to be readily accessible to experimentalists. Let 
us first review a few types of models that could be utilized for GPDs. This will lead us to 
double distributions (DDs) which are taken up in subsequent Chapters. 

Above we have seen the difficulties encountered when applying light-cone wavefunctions 
to calculate GPDs. The kinematical region — ^ < x < ^ necessarily involves higher Fock 
components of the bound state. In light-front phenomenology these are rarely considered. 
Although we derived these components perturbatively in the Wick-Cutkosky model, one 
ultimately must sacrifice Lorentz covariance to make the problem numerically tractable, 
even at weak coupling. Violation of Lorentz symmetry is undesirable for GPD models. 
The essence of the distributions is the dependence on the skewness variable ^. In terms 
of the form factor, we can understand GPDs via the calculation of the form factor in an 
arbitrary frame labeled by ^ (we will appreciate this better in Chapter Viewed from 
this perspective, the sum rule's lack of dependence on ^ is due to the freedom to choose 
arbitrary frames in which to calculate the Lorentz invariant form factor. A virtue of GPDs 
is the ability to measure correlations which otherwise vanish due to Lorentz invariance in 
single current matrix elements. The polynomiality restrictions of GPD moments are further 
manifestations of Lorentz invariance and are additionally spoiled by light-front calculations. 

These observations about light-cone wavefunction calculations of GPDs are generally 
true of quark models as well. The reason quark models fail to be useful for GPDs stems 
from their non-relativistic nature coupled with the restriction to a fixed number of particles. 
If one were to utilize a Poincare covariant framework to calculate GPDs in terms of a fixed 
number of constituents, the kinematic region — ^ < a; < ^ will be a complication. In this 
framework, one must ultimately deal with constructing dynamical operators consistent with 
Poincare covariance in order to express the GPDs. In order to ensure Lorentz invariance, one 
must return to the general two photon amplitude and construct according to the Poincare 
symmetries. This procedure is likely to be inconsistent with asymptotic freedom as the 
dynamics required for Poincare invariance of the two photon amplitude have nothing to do 
with perturbative QCD in the asymptotic limit. 
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Figure 3.3: Proper gauging of the DSE for the three-point function with the non-local light- 
cone operator (denoted by a cross). The Ball-Chiu Ansatz cannot be used to resolve the 
active quark . 

In non-relativistic scenarios, one also has the problem of support. Ordinary PDFs cal- 
culated in non-relativistic quark models extend beyond x G [0, 1] and this trouble naturally 
carries over to GPDs. Moreover, one is lead non-relativistically to little or no contribution 
in the region — ^ < x < ^, and the distributions will likely vanish at the crossover x = ^. 
Removing these weaknesses of model GPDs by using evolution is unsatisfactory because one 
then relies solely upon perturbative QCD for the non-perturbative physics absent from the 
original model. Besides, this is considerable work when alternatively the gross qualitative 
features of PDFs can be obtained more simply from suitable evolution of the distribution 
fi{x) = 6{x — 1/3). By analogy one can forgo the entire the issue of ^ dependence, use the 
Ansatz: H{x,^,t) = fi{x)Fi(t), and hope for the best. Ultimately this parametrization 
will not be quantitative enough for interpreting data. 

Beyond quark models, an approach that we considered in some detail was the calculation 
of GPDs from Dyson-Schwinger equation (DSE) models [65, 66]. The Dyson- Schwinger 
approach features Lorentz covariance and attempts to model self-consistently the solutions 
to field theories. Both of these qualities are desirable from the point of view of GPDs because 
relativity is borne in and there is no fixed particle restriction in field theory. Unfortunately, 
current Dyson-Schwinger model approaches, which fall into two classes, are incompatible 
with the properties of GPDs as we now explain. 

The first class of DSE models involves solving the covariant DSEs in a truncation scheme, 
such as the rainbow-ladder truncation. The solution of these equations requires treatment 
in Euclidean space and the connection to light-cone dominated amplitudes is hence gen- 
erally lost. The singularity structure of propagators and vertices must be known in the 
whole complex plane to enable calculation of PDFs or GPDs. Clearly this cannot be done 
with a numerical solution in a limited domain. Imagining that GPDs could be obtained 
in the rainbow-ladder truncation scheme, one has the additional complication of the pos- 
itivity bounds. It is unclear whether non-perturbative truncation schemes are consistent 
with positivity. Part of this uncertainty stems from a lack of knowledge of the analytic 
continuation, in particular the way to assess contributions from the infinite light-cone sin- 
gularities of the bound-state vertex. On the other hand, an interesting application of these 
DSE models would be to calculate the moments of PDFs and GPDs since the calculation of 
matrix elements of twist-two operators can be carried out in Euclidean space. Surprisingly 
this territory is largely unexplored. 

The second class of DSE models relies upon using Ansdtze for the various vertex func- 
tions and propagators. Additionally field theoretic constraints arising from Ward-Takahashi 
identities are imposed, although the solutions to these constraints are not uniquely specified 
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by the model elements. This class of models has the same drawbacks for calculating GPDs 
as the first class. One again has difficulty connecting the Euclidean space formulation of 
these models with light-cone correlation functions. Moreover, entire functions are often 
used for model elements and these in turn defy a description of PDFs and GPDs because 
the integral over the relative light-cone energy does not converge in the complex plane. It 
is clear for this class of models that the positivity bounds must be generally violated. This 
is because an ad hoc bound state vertex cannot be gauged with the non-local light-cone 
operator. Correctly inserting the non-local operator at the bounds state vertex, which en- 
sures the positivity bound is respected, requires knowledge of the underlying field theoretic 
dynamics. 

Lastly there is an additional problem present in this second class of DSE models due to 
the incorrect treatment of higher point functions. The well-known Ball-Chiu form for the 
electromagnetic vertex [67] is used in these models to describe an electromagnetic three- 
point function that is consistent with the Ward-Takahashi identity. The Ball-Chiu vertex, 
however, cannot be used for PDFs or GPDs because one needs to gauge the three-point 
function with the non-local light-cone operator and this requires knowledge of the four- 
point function [29,30,68], see Figure 1213 As is, the Ball-Chiu Ansatz cannot fix the plus 
momentum of the active quark because it models the sum of the ungauged diagrams in 
the Figure. The individual contributions cannot be resolved and hence cannot be properly 
gauged. 

Clearly the consistent modeling of GPDs is a difficult problem. Most standard hadronic 
models cannot be used in their present forms to calculate GPDs. Furthermore investigating 
how to rectify such models with the properties required of GPDs is lengthy and not par- 
ticularly illuminating, when one considers that the models have no transparent relation to 
QCD. Thus it makes sense to consider only simple models whereby one can guarantee the 
properties of GPDs will be satisfied from the onset, and the resulting parametrization will 
be efficaciously tractable for comparison with experiment. To devise such models, we turn 
our attention to double distributions. 
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Chapter 4 

Double Distributions 



Modeling GPDs directly from their definition leads to light-cone quantization and light- 
cone wavefunctions. Another approach is to use the formalism of double distributions 
(DDs) [69, 70] and this is the discourse taken up in the rest of this work. The formalism 
of DDs elegantly explains the polynomiality conditions required of GPDs and thus gives 
one the ability to construct models consistent with known properties — although insight into 
model construction has often been limited to factorization Ansdtze. Very few DDs have ac- 
tually been calculated, although ad hoc DDs are almost exclusively used for phenomenology. 
Recently two-body light-front wavefunction models of the pion have been used to obtain 
GPDs [71] based on DDs. Without modifying the quark distribution, the resulting GPDs 
violate the positivity constraints [63]. This inconsistency is attributed to missing contribu- 
tions from non-wavefunction vertex diagrams the contribution of which are unknown when 
one uses non-covariant vertices. In general these diagrams are a substitute for higher Fock 
components, see [72]. Thus the method used by [71] must at least use covariant wavefunc- 
tions to be sensible. The covariant models we used [73], however, allowed us to test the 
uniqueness of this ostensible construction. In this Chapter, we present the two models which 
were used to test the construction of DD functions. The two models are both simple models 
for scalar bound states, one with scalar constituents and the other with spin-^ constituents. 

Indeed we found that appealing to Lorentz invariance is enough to determine only one 
component of the double distribution in the two-component formalism (even for C-odd 
distributions, where the Polyakov- Weiss D-term [74] is absent). Moreover, the component 
determined from the reduction relations in the scheme presented by [71] is completely am- 
biguous. Below we show that in the scalar constituent model, missing the second component 
leads to incorrect GPDs. The same is true for a spin-^ constituent model of the pion, where 
additionally the positivity bound is violated. The correct DDs unique to each model can be 
calculated from non-diagonal matrix elements of twist-two operators and is demonstrated. 
On the other hand, exploiting the ambiguity inherent in defining one component DDs (which 
is akin to gauge freedom [75]) one can generate infinitely many different GPD models which 
share the same form factor and quark distribution as well as satisfy polynomiality (and 
likely positivity). We also calculate DDs for a model with a propagator having complex 
conjugate energy poles. Such a parametrization is an efficacious way to model confinement. 
The calculation produces no trouble in Euclidean space, while the difficulties in Minkowski 
space are dealt with in Appendix O 

The organization of the Chapter is as follows. After defining DDs in Section 14.11 we 
explicitly derive the GPD for the scalar triangle diagram with point-like vertices in Section 
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Km This is the GPD for the toy BSE model in Section THm Next we show the DD 
for this model extracted from the form factor in the Drell-Yan frame via [71] leads to 
incorrect GPDs. Having encountered this problem, we calculate the missing component 
of the DD for the scalar triangle diagram in Section 14.2.21 Then in Section 14.2.31 we 
present a calculation of the DD for a similar model where the propagator is extended to 
have complex conjugate poles. The calculation is done in Euclidean space here and in 
Minkowski space in Appendix [Ul In Section [4.31 we move on to the model for the spin-i 
case. Next in Section [4.3.21 we regularize the current and then extract this model's GPD. 
Although not manifest, this model satisfies polynomiality, which is demonstrated in Section 
14.3.31 Using [71] as a guide, we construct a DD for this model in Section 14.3.41 Similar 
to Section [4. 2. H this one-component DD too gives rise to a different GPD than the light- 
front projection. Additionally positivity is not satisfied by this one-component DD (Section 
14.3.5j) . We calculate the complete two-component DD from matrix elements of twist-two 
operators in Section [4.3.61 Lastly we conclude with a brief summary (Section 14.41) . 



4.1 Definitions 

In this Section we define DDs for scalar bound states via their moments by following the 
two-component formalism of [74,75]. Focusing on the ambiguities inherent in defining one- 
component DDs from two component objects, we will understand why the DD constructed 
according to [71] leads to an incorrect GPD. Moreover, we shall calculate the correct DD 
for the scalar triangle diagram from matrix elements of twist-two operators. We remark 
in passing that the two components of the DD (below F and G, or F and D-term in the 
standard formalism) can be viewed as projections of a single function of two variables [76]. 

In defining GDPs and DDs there are two convenient choices of variables to use. Above 
we have worked in Radyushkin's asymmetrical frame and asymmetrical variables which 
are ideal for perturbation theory and are a natural generalization of ordinary parton dis- 
tributions. The non-diagonal matrix elements of twist-two operators are, however, more 

conveniently expressed in variables symmetric with respect to initial and final states.^ To 

> — > < — 

this end we define the average momentum P^' = {P + P'y/2. Let D^' = ^{d^ - (9^). Then 
for a pion of spin-i constituents we have 

{P + A/2\^i;{0)-f{^'iD^'' ■ ■ ■ W''"} xP{0)\P - A/2 ) = 

where the action of ^ on Lorentz indices produces only the symmetric traceless part and 
T-invariance restricts k in the first sum to be even and odd in the second. For a pion 

of scalar constituents, replace with 2iD^ . As it stands there is considerable freedom 
in the above decomposition, e.g., one could rewrite the above with kBn^k-iit) / {n — k + 1) 
as a contribution to Ajik{t)- Carrying this out for all k, puts the bulk in the first term 

^Good discussion of the conversion from symmetrical and asymmetrical variables and distributions can 
be found in [77]. Additionally advantages and disadvantages of both are presented. 
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and renders the second term proportional to the symmetric traceless part of only (n + 1) 
A's — moments of the Polyakov- Weiss D-term [74]. This is the usually encountered form 
of the DD with D-term. Calculationally, however, we find our results directly in the form 
of Eq. Kn . 

The DDs are defined as generating functions for the moments^ 

Ank{t)= d(3 da P''-''a''F{/3,a;t) (4.2) 

J-l J-l+\l3\ 

1 ri-m 

•n—k„ k/ 



Bnk{t)= d(3 da P''-^a''G{P,a;t). (4.3) 

As a consequence of the restriction on k in the sums, the function -F(/3, a; t) is even in a 
while G{(3, a; t) is odd. Also for n-even, there is no contribution from the I?-term to the 
function G{P, a; t). 

These DD functions then appear in the decomposition of matrix elements of light-like 
separated operators. Summing the moments to produce these operators, we have 

( P + A/2 I 'tp{-z-/2) / ip{z-/2) I P- A/2 ) 

= 2P-z dp da e-^/^^^-^+^^^-^Z^ F(/3, a; t) 

J-i J-i+m 

/■I fl-|/3| - 

-A-Z dp da e-^^^-^+*"^-^/2 G{P, a; t), (4.4) 

J-l J-1+\I3\ 



where = 0. 

Denoting ^ = — A^/2P"'", the GPD in symmetric variables reads 

Hix,U) = ^ I dz-e^^^'^'- {P + A/2\iji-z- /2)j+i^{z- /2)\P - A/2). (4.5) 

Inserting Eq. (|4.4j) into this definition yields 

-1 rl 

-1+ 

By integrating over x, we uncover two sum rules: the sum rule for the form factor 



/■I /-l-l/SI r T 

H{x,C,t)= dp da5{i- P-ia) F{p,a;t)+iG{P,a]t) . (4.6) 

J-l J-i+\l3\ 



j dp j daF{P,a;t) = F{t) (4.7) 
and what we call the G-sum rule 

^ dp I da G{P,a;t) = 0, (4.8) 



which is trivial since G is an odd function of a. Eq. 1)4. 8|) has non-trivial consequences 
however, e.g., it shows the method employed by [71] leads only to the F DD in Eq. (|4.4j) . 
This function integrates to the form factor via Eq. 1)4. 7|) and in the forward limit {^, t — > 0} 
reduces to the quark distribution (when integrated over a). Thus F{P,a;t) should be 

^In general the twist-two form factors and hence the double distributions have renormalization scale 
dependence that arises from perturbative QCD corrections. 
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p (b) P' 

Figure 4.1: Covariant triangle diagram for the meson electromagnetic form factor. 



properly termed the forward- visible DD, which encompasses more than just neglecting the 
D-term. From Eq. 1)4. we see that G{(3,a;t) does affect higher moments of the GPD. 
This will be the source of the discrepancy between the method employed by [71] and the 
correct calculation of the DDs from Eq. 1)4. 1|) . We now calculate model DDs and expose 
this discrepancy. 

4.2 Simple models 
4.2.1 Scalar model 

For the meson model with scalar constituents, we choose a point-like Bethe-Salpeter vertex 
T{k, P) = 1, where the coupling constant is assumed to be absorbed into the overall normal- 
ization. This is the model considered in Section [1.2.21 Furthermore, we choose derivative 
coupling of the photon to charged scalar particles. 

The meson electromagnetic form factor for this model can be calculated from the Feyn- 
man triangle diagram. In order to derive the GPD, however, we need to choose the kine- 
matics specified in Figure 14.11 with k as the momentum of the struck quark. Using the 
stated meson vertex and taking the plus-component of the current, we have 

Fit) = I d'k=^^^^^±^{[k'-m' + ze] [{k + Af - + ^e] [{P - kf - + ^e]y\ 

(4.9) 

where the momentum transfer is t = and the skewness is defined relative to the ini- 
tial state A"*" = — C-P"*" < 0. Physically C plays the role of Bjorken variable for DVCS. 
Additionally we work in the asymmetric frame where P"*" = 0. 

To turn Eq. (|4.9j) into an expression for the GPD H{x, t), we insert 6{k'^ /P^ — x) to 
fix the momentum of the struck quark and keep C 7^ 0- This forces 

m = lm.,(.t)^.. (4.10) 

Lastly we integrate over k^ to project onto the light cone. Doing the contour integration 
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to extract H{x,C,,t) in Eq. H4.9() . we are confronted with the poles 

' u- — u- _ H 

k-=p- + {k-P)-^-^^ , (4.11) 
,fc- = -A- + (fc + A)„- -f 

where the on-shell energies are = ^ 2p+^ the abbreviation x' = is used. Thus 
the non-vanishing contribution to the integral is 

2m e[x{l - x)] [e{x - C) Res(A:^-) - ^(C - x) Res(A;-)] , (4.12) 

which leads to 

(1 - (:/2)H{x, C, t) = e{x - C)Hi{x, C, t) + 9{C - x)H2{x, (, t). (4.13) 

Using k'-*- = k"*- + (1 — x') A-*- for the relative transverse momentum of the final state, the 
functional forms are 

Hi{x,C,t) = {2x-0\N\^ j dk^Dw(x,k^|M2)Dw(x',k'^|M2)/x(l -x)x' (4.14) 
H2{xX,t) = i2x-C)\N\^ J dk^Dw{x,k^\M^)Dw{x" ,k"^\t)/Cx"{l - x"){l - x), (4.15) 

where x" = x/( and k"-*" = k"'"-|-x" A"*" are the relative momenta of the photon. Additionally, 
we use the replacement 

DMx,k^\M'r' = M'- (4.16) 

x(l — X) 

which is the propagator of the Weinberg equation [44]. 

Comments about the GPD H{x, C, t) in Eq. (1^131) are in order. The model is covariant 
and thus the sum rule and polynomiality conditions are met (see Section 14.3.31 below for 
clarification). We have checked this explicitly and suitable discussion can be found in [78,79]. 
Consideration of this model without derivative coupling was first done from the perspective 
of DDs, see e.g. the toy model of [69]. This DD model was revisited recently with derivative 
coupling at the photon vertex in the Appendix of [71] and the same DD also appears in [78]. 
To derive the DD for this simple model, we first appeal to Lorentz invariance as suggested 
in [71], recalling along the way the relevant properties of DDs in asymmetric variables. 

First consider the form factor. In the C = frame, Eqs. (|4.14M1?T3|) reduce to the 



Drell-Yan formula [36,80] via the definition in Eq. I|4.1 

-^^—-^DMx,k^\M^)Dy,{x,k^ + (1 - x)A^|m2), (4.17) 

where t = — A^. The form factor is Lorentz invariant and has a decomposition in terms of 
the Lorentz invariant DD F{x,y;t), namely 

F{t) = [ [ F{x,y;t)dydx. (4.18) 
Jo Jo 
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Figure 4.2: Comparison of covariant GPDs for the scalar triangle diagram. The GPDs 
Eq. (|4.1H|) (denoted LC) and Eq. (|4.2:-{|) (DD-based) are plotted as a function of x for fixed 
C, = 0.7 and t = —rin? for the mass = m < 2m. We also plot the difference between 
the two curves (6). The area under the curves is identically F{—m'^) for LC and DD-based 
GPDs, and hence zero for their difference 5. 



The DD satisfies the following relations [70]: support property 

F{x,y;t)^9{l-x-y), (4.19) 
reduction to the quark distribution at zero momentum transfer 

pl—x 

q{x) = / F{x,y;0)dy (4.20) 

and is Miinchen symmetric [81] 

F{x,y;t) = F{x,l-x-y;t). (4.21) 

Using Eq. (|4.17)1 . we can write F[t) in the form (|4.18)) with 

- x\N\'^6{l - x - y) 

^' ^' ~ w?-x{l-x)M^-y{l-x-y)t' ^ ^ ^ 

which satisfies Eqs. H4.19fH^^T]) . The ingenuity of DDs comes about when we construct the 
GPD via ^ ^ 

HixX,t)= [ [ F{z,y-t)5{x-z-Cy)dydz. (4.23) 

JO JO 

In this form the sum rule and polynomiality conditions follow trivially. 

In Figure 14.21 we plot the GPD Eq. 1)4. 13|) as weh as the GPD derived from DD via 
Eq. (|4.23)) . Surprisingly the two are different despite the fact both models are covariant and 
posses the same form factor and quark distribution. Additionally we plot their difference 
(5) as a function of x. 

4.2.2 Scalar model, revisited 

Based on the definitions in Eq. 1)4. 1() . we need more than the form factor to determine the 
GPDs. The procedure above thus missed one component of the DDs. To correctly derive 
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both F and G DDs for the scalar triangle diagram of Section 14.2.11 we must now consider 
the action of the operator 2iD^^^iD^^'^ ■ ■ ■ iD^"> . This produces a factor 

— {2k + A)^''{2k + A)''^ ■■■{2k + A)^"> (4.24) 
2"- 

in the integrand of Eq. (|4.17() . which we now take in the symmetric frame. After the 
integration over k is performed, we are left only with 

rl fl-l3 _ _ 

I d(3 da D{f3,a;t){2(3P -aA)^f'{f3P -aA/2)f'^ ■■■{PP -aA/2Y"\ (4.25) 

Jo J-i+p 

where we have used the replacement 

= m--fi{X-mi^-t\{^-fif-o?\l^ ^'-''^ 

Using the binomial expansion, we can identify -F(/3, a; i) and G(/3, a; t) via Eqs. (|4.2() 
and ()4.3() . namely'^ 

F(/3, a;t) = fi Q{\ - fi - a)D{P, a; t)\N\'^ (4.27) 
G(/3, a;t) =ae{l-(3- a)D{p, a; t)\Nf. (4.28) 

To compare with the results of Section r4.2.1l we must revert to asymmetric variables which 
is accomplished by the transformation /3 — > x, y — > (a — /3+l)/2. The denominator common 
to both terms becomes 

D{x, y; t) = — — ^ -, (4.29) 

— x(l — x)M^ — y{l — x — y)t 

and we have 

F{x, y;t)=x9{l-x- y)D{x, y; t)\N\^ (4.30) 
G{x, y; t) = {2y + x- 1) e{l-x- y)D{x, y; t)\N\\ (4.31) 

Notice the function G{x, y; t) is Miinchen antisymmetric, i.e. G{x, y; t) = ~G{x, l — x — y;t), 
which is required because G{P, a; t) is odd with respect to a. 

To construct the GPD H{x,C:t) we must also convert Eq. 1)4. 6(1 to asymmetric variables. 

H{x, C, t) = C dz C ^ dy5{x-z- Cy) \f{x, y; t) + ^Mx, y; t)] (4.32) 
Jo Jo L ^ — t, J 

We can now plot the GPD calculated from Eq. (|4.32|) using the two DDs in Eqs. 1)4. 30() and 

(|4.31j) . The result agrees with Eq. (|4.13|1 depicted in Figure E21 Moreover, the contribution 

from is identically the difference 5 plotted in the figure. In the DD formalism, the 

reduction relations do not determine the GPD. We will illustrate this further with a spin-i 

model for the pion. Before doing so, however, we will take a detour to pursue an interesting 

extension of the above scalar model. 



■^Notice in this case we can see that F and G are projections of a single function D{/3, a; t) = /(/?, a) in 
the notation of [76]. Determination of f{[3,a) for the spin-i model of Section f4.3l however, cannot be done 
by inspection. 
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4.2.3 Complex conjugate poles 

In this Section, we extend the calculation of the model DD above to include complex 
conjugate singularities in the constituents' propagators. Complex conjugate singularities 
present in solutions to Dyson-Schwinger equations have been studied in the connection 
with the violation of Osterwalder-Schrader reflection positivity and confinement [82-89]. 
Recent work [90] in solving the Bethe-Salpeter equation with a quark propagator consisting 
of pairs of complex conjugate singularities shows that the width for meson decay (into free 
quarks) generated from one pole is exactly canceled by the contribution from its complex 
conjugate. Additionally recent studies have modeled Euclidean space lattice data with 
propagators that have time-like complex conjugate singularities [91-93]. In [94] we pursued 
the calculation of space-like amplitudes in Minkowski space using a simple model propagator 
consisting of a pair of complex conjugate poles. Specifically, we were interested in the 
calculation of light-cone dominated amplitudes for this type of model which necessitates 
a treatment in Minkowski space. This treatment is contained in Appendix [H] because 
there is considerable subtlety. Nonetheless, covariant calculation in Euclidean space is 
straightforward and provides another chance to obtain DDs properly. 
In Euclidean space, the model propagator is chosen to be 

SM*) = Ejj^^. (4.33) 

Here and below we use the shorthand e = it to denote the pair of complex conjugate 
poles. Notice e is finite and can be chosen to be positive. Unlike in Minkowski space where 
the measure is imaginary, contributions to Euclidean space amplitudes are real and one 
has no difficulty in calculating form factors and distribution functions using Ssik) in the 
relevant diagrams. The simplicity of the model at hand will allow us to calculate its double 
distribution analytically and thereby determine the quark distribution and electromagnetic 
form factor, since these functions are related to the double distribution by the so-called 
reduction relations. The difficult task of calculating of these quantities in Minkowski space 
by projecting onto the light-cone is contained in Appendix O 

GPDs are not Lorentz invariant objects, however, they stem from a projection of a 
Lorentz invariant double distribution function [69] . These functions are particularly attrac- 
tive from the perspective of model building [71], though one must be careful that the starting 
point is indeed covariant [63], otherwise desirable distribution properties and straightfor- 
ward physical interpretation may be sacrificed. The model under consideration is fully 
covariant, and thus the DD representation is an ideal testing ground for our model propa- 
gator. Hence we proceed to calculate the model's Euclidean space DD, recalling along the 
way the relevant properties of DDs. In this Section, by contrast to Eq. (|4.H) . we decom- 
pose the matrix elements of twist-two operators in the asymmetric frame, with asymmetric 
variables^ _^ ^ 

Let L"^ = d^ — d^. For this scalar model, we define the twist-two operator of spin-(n-|- 1) 
as above 

Q/^w-A'n ^ (j){{))iD^^'iD^'^ . . . iB^">(/.(0), (4.34) 

where the action of ^'"^ on Lorentz indices produces only the symmetric traceless part. 

We work in Radyushkin's asymmetric frame with P as the momentum of the initial 
state, P + A that of the final and t = A^. The initial and final states are on-shell: = 
(P + A)^ = M^. The non-diagonal matrix element of Qi^i^i—t^n ^^^^ |-,g decomposed into 
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P + A 



Figure 4.3: Diagram used to calculate non-diagonal matrix elements of twist-two operators 
(denoted by a cross). 



Lorentz invariant moment functions Ankit) and Bnk{t). Asymmetrically this decomposition 
reads 

( P -)- /\\Q^^^^^■■■^^"\p ^ = 

(2P + A)^'' ^' ,., ^nfc(t) (2P + A)''' ■■■{2P + AY'"-" (-A)''"-'=+i • • • (-A)^"^ 

- A^^y , ^' , Bnkit) (2P + A)^i •••(2P + A)'^"-'=(-A)^"-^+i •••(-A)^">, (4.35) 
Klfn — ky. 

Hermiticity forces the matrix elements of Qi^i^'^—i^n to be invariant under the transformation 

P ^P ^ A 
A^-A 

Consequently the values of k are restricted to be even in the first sum and odd in the second. 
As it stands there is again considerable freedom in writing Eq. ()4.35p . e.g., one could rewrite 
the above with kBn^k-\(i) j (n — A; -1- 1) as a contribution to A„fc(t). Carrying this out for 
all /c, puts the bulk in the first term and renders the second term proportional only to the 
symmetric traceless part of (n-|- 1) A's — moments of the Polyakov- Weiss D-term [74]. This 
is the usually encountered form of the DD with D-term. Alternatively one can also express 
the moments as projections of a single Lorentz invariant function [76]. Calculationally, 
however, Eq. (|4.35j) becomes the most practical to work with [73]. 

The F and G DDs can be defined as generators of the coefficient functions 

Ank{t) = C dx C " dy x"-^(x + 2y - lfF{x, y; t) (4.36) 
Jo Jo 

Bnk{t)= f dx f dyx^-\x + 2y-lfG{x,y;t). (4.37) 
Jo Jo 

As a consequence of the restriction on k in the sums, the function F{x^ y; t) is Miinchen 
symmetric [81] , i.e. F{x,y;t) = F{x,l — x — y;t), while G{x, y; t) is Miinchen antisymmetric. 
Also for n-even, there is no contribution from the D-teim to the function G{x, y; t). 
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These functions then appear in the decomposition of matrix elements of hght-Hke sepa- 
rated operators 

( P + A I 0(0) iz-D^{z-) \P) = 

{2P-Z + A-Z) C dx C \y e-'''^-'+'y^-'F{x,y;t) 
Jo Jo 

-A-z [ dx f dy e-"-^-^+*^^-^G(x, y; t), (4.38) 
Jo Jo 

where z'^ = 0. 

Denoting = — > 0, the GPD in asymmetric variables reads 

2^i2-Q ^ ^ + ^ ' m^D+^in I P )■ (4.39) 

Physically C, plays the role of Bjorken variable for deeply virtual Compton scattering. In- 
serting Eq. 1)4. 38() into this definition yields 

H{x,C,t)= [ dz [ dy6{x-z-Cy)\F{z,y;t) + -^G{z,y;t)\. (4.40) 
Jo Jo ^ — C J 

which is what we found circuitously above. By integrating Eq. 1)4. 40(1 over x, we uncover 
two familiar sum rules for the DDs: the sum rule for the form factor 

[ dx [ dy F{x,y-t) = F{t) (4.41) 
JO Jo 

and the G-sum rule ^ ^ 

[ dx I dyG{x,y;t)=0, (4.42) 
Jo Jo 

which follows after a change of variables since G is Miinchen antisymmetric. Eq. (|4.42|) is 
important and mandated by current conservation. As we saw above, the G DD function is 
all too frequently overlooked and treated as identically zero. Lastly, the model quark distri- 
bution function q{x) can be found from the DD at zero momentum transfer, see Eq. (|4.39() . 

qix) = / dyF{x,y;0). (4.43) 



Jo 

We can use the decomposition in Eq. (|4.35|) to calculate our simple model's DD. Pa- 
rameterizing the momenta as in Figure the non-diagonal matrix element of O^") reads 

( P -)- /\\Q^^^^^■■■^^"\p ) = 

2N ^ r (2fc + A)i/-(2fc + Ar • • • (2fc + A)/-"> 

7r2 ^ y Ik"^ + m? - ie] \{k + AY + - ie'] \{k - PY + m? - ie"] 

e,e',e"=it 

The normalization constant N is chosen by the condition -F(O) = 1. Let us denote the 
propagators simply by 2t = {k — PY +m? — ie" , 55 = {k+ AY +itY — is' and C = k"^ -\-rnY — ie. 
We introduce two Feynman parameters {x, y} to render the denominator specifically in the 
form [xSt -|- y*B + {1 — x — y)^]~^. One then translates k^ to render the integral (hyper-) 
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spherically symmetric via the definition = l^^ + xP^ — yA^. The resulting integral over 
/ can be evaluated directly (remember we are in Euclidean space). 

Binomially expanding the result of the integral, we can make contact with Eq. ()4.35|) and 
subsequently determine the F and G double distributions by inspection from Eqs. 1)4. 36j) 
and (|4.37|) . Defining the auxiliary functions 



Do{x, y- t)=m^ - x{l - x)M'^ - y{l - x - y)t (4.45) 

and 

ri(rr irA — AT \^ 

Z)„(x,y;i)2+e2(i_2z)2' 



z=0,x,y,x+y 

the DDs can be written simply as 

F{x,y;t) = xD{x,y;t) (4.47) 
G{x, y; t) = {x + 2y- l)D{x, y; t). (4.48) 

We treat factors of 9{\ — x — y) as implicit above. Accordingly F is Miinchen symmetric and 
G is antisymmetric. Notice although e is finite, corresponding results using the standard 
perturbative propagator can always be recovered in the limit e — > 0. For example, the 
correct F and G DDs are recovered in the limit e ^ [73]. 

The model GPD can be derived by utilizing Eq. (|4.4Uj) . although the integral must be 
performed numerically. The quark distribution can be found via the reduction relation 
Eq. namely 



g(x) = iV ^ 



z=0,x 



x{l -x)Do{x, 0;0) X ^^e{l-2z) 



(4.49) 



Lastly the form factor can be found from the sum rule Eq. 1)4. 4ip 

In Figure 1^31 we plot the quark distribution and electromagnetic form factor for various 
values of e in GeV^. We have arbitrarily chosen the other model parameters as M = 0.14 
GeV and m = 0.33 GeV. Additionally in Figure 14.51 the GPD is plotted: first at fixed 
C and t for various values of e and then at fixed t and e for various values of C- Curves 
corresponding to e = are the standard results for a propagator with one real pole. 

Above we considered calculation of amplitudes for space-like processes using a scalar 
propagator with one pair of complex conjugate poles. Moreover the analysis can be extended 
easily to the case where vertex functions have complex conjugate singularities. Such models 
cannot be directly employed in Minkowski space, they must be analytically continued from 
Euclidean space. 

In Appendix [GJ the problems of using a propagator with complex conjugate poles in 
Minkowski space are discussed at the level of the quark distribution function. If the model is 
defined in Minkowski space, one will generally violate the support and positivity properties 
of the quark distribution. Above we showed the model is perfectly well defined in Euclidean 
space by calculating non-diagonal matrix elements of twist-two operators. This leads us to 
the model's double distribution which we used to calculate parton and generalized parton 
distributions as well as the electromagnetic form factor. 
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Figure 4.4: On the left, the quark distribution Eq. H4.49jl is plotted as a function of x for 
a few values of e in GeV^. On the right, the form factor calculated from Eqs. ()4.41|) and 
(|4.47|) is plotted as a function of —t for a few values of e in GeV'^ . The model parameters 
are arbitrarily chosen as: M = 0.14 GeV and m = 0.33 GeV. 




Figure 4.5: Plots of the GPD calculated from Eqs. (|07l) . and (|On)) . On the left, 
the GPD is plotted as a function of x for a few values of e (in Gev^) at fixed C = 0.5 and 
t = -1.0 GeV2. On the right, the GPD appears at fixed e = 0.17 GeV^ and t = -2.0 GeV^ 
and is plotted as a function of x for a few values of The model parameters are arbitrarily 
chosen as: M = 0.14 GeV and m = 0.33 GeV. 
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4.3 Double Distributions for the pion 



Having calculated model DDs for scalar constituent models, we shall now add spin degrees 
of freedom to the constituents. In this Section, we explore a model for the pion both 
on the light cone and from the covariant description in terms of DDs. To calculate the 
DD functions, we again highlight the ambiguity and limitations inherent in the method 
employed by [71] and then properly determine the DDs. 

4.3.1 Bethe-Salpeter amplitude and pion wavefunction 

For the spin-^ model, we choose the trivially qq symmetric Bethe-Salpeter vertex 

T{k,P) = -igj'^. (4.50) 

Here we have assumed only 7^ coupling at the quark-pion vertex with coupling constant g, 
whereas four Dirac structures exist [95]. This simple coupling is suggested by an effective 
interaction Lagrangian of the form (see, e.g., [96,97]) 

£j = -ig TT . qj^rq, (4.51) 

where the coupling constant g = m//7r, with m the constituent mass and /j^ the pion decay 
constant. Notice the (ladder approximation) kernel is independent of light-cone time. Thus 
this model (as well as the scalar triangle model in Section 14.2.1(1 are special cases of the 
instantaneous formalism described by [30] in the impulse approximation. The relation of 
the vertex to the Bethe-Salpeter wavefunction is given by 

^{k, P) = \. {-ig)l \ ' (4.52) 

The valence wavefunction can be found by projecting the Bethe-Salpeter wavefunction 
onto the light-cone = 0, see e.g. [98]. Using the normalization convention of [3], we have 



V'(rc,k4i; A, A') 



2P+ 



27r 



(4.53) 

where x is the fraction of the pion's plus momentum carried by the quark (x = k'^/P'^), 
and the relative transverse momentum is k; 
found from Eq. (|4.53() to be 



rcl 



xP . The valence wavefunction is 



V^(x,k^;A,A') = ^\f^^ D^{x,k^\M^) k^x Sx.y - Xm 5x,^x' 
x[l — X) 



(4.54) 



where we have employed the notation kx = k^ + iXk'^. As a result of the contour integra- 
tion, we have a factor of 9[x{l — x)] implicitly in Eq. ()4.54|) . Additionally the wavefunction 
is symmetric under interchange of x and 1 — x. As is known, introduction of orbital an- 
gular momentum into this wavefunction leads to divergent quark distributions and form 
factors which will be handled below. Since this model is non-renormalizable, the choice of 
regularization scheme influences the dynamics. 
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4.3.2 Form factor and generalized parton distribution 

The pion electromagnetic form factor for this model can be calculated from the Feynman 
triangle diagram. In order to extract the GPD, however, we need to choose the kinematics 
specified in Figure 14.11 with k as the momentum of the struck quark. As it stands, using 
the wavefunction Eq. ()4.54|1 . the triangle diagram diverges. Following the approach of [99], 
we covariantly smear the point-like photon vertex in Figure 14.11 in a way reminiscent of 
Pauli-Villars regularization 

Eq. (|4.55j) is a simple way to model non-qq components of the wavefunction.^ Alternatively 
one could smear the qq — ir vertex in a covariant manner [31,32, 100]. This smearing should 
additionally respect the qq symmetry of the vertex. We do not pursue this option here since 
positivity constraints (see Section I4.3.5() are generally violated. On the other hand, one 
could use Pauli-Villars subtractions to regulate the theory, however, positivity would also 
be put into question. Because our concern is with model comparisons not phenomenology, 
we shall choose K = m merely for simplicity. Although not obvious from inspection, results 
for h.^ m exhibit the same features investigated below. Most noteworthy, positivity remains 
satisfied when m. 

Considering matrix elements of the current operator between pion states, the model 
Eq. 1)4. 55() conserves current. This can be demonstrated most easily by calculating A • J 
in the Breit frame. Additionally since the model is fully covariant, we can extract the 
electromagnetic form factor from any component of the current. In particular, potential 
end-point singularities present in matrix elements of J~ have been removed by the photon 
vertex smearing Eq. 1)4. 55|) [99]. Using the plus-component of the current, we have the 
expression 

F{t) = ^ I ' 



+ m)-f^{jt - p + m)7^(^ + A + "T')7^ 



1-| J (2^)4 [A:2-m2 + k]^[(A: + A)2-m2 + ie]^[(P-A;)2-m2 + ie]' 

(4.56) 

where the momentum transfer is t = A2 and the skewness C is defined relative to the initial 
state: A+ = -(P^ < 0. 

To calculate the GPD, we follow the procedure described in Section 14.2.11 The result 
can be written as 

(1 - C/2)H{x, C, t) = e{x - OHesix, C, t) + 0(C - x) [{Hi^stix, C, t) + /7nval(x, C, t)] . (4.57) 

where is the piece determined by the effective two-body wavefunction, i/inst is the 
contribution from instantaneous propagation of the spectator quark, and the remaining 
contributions we term non-valence (although strictly speaking the instantaneous piece is 
also of the non- valence variety). It is a peculiarity of this model that explicit instantaneous 
terms are not present for x > (. 

^For all its positivity preserving virtues, the vertex smearing has unphysical drawbacks. Firstly the 
Ward-Takahashi identity is messed up. Secondly the large momentum behavior of the induced constituent 
quark form factor is inconsistent with asymptotic freedom. We will come back to these issues in Chapter]^ 
for the proton. 
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Figure 4.6: On the left, the GPD Eq. H4.57() is plotted at fixed t = —m? for a few values of 
C^. On the right, the same GPD is plotted for fixed C = 0.5 for a few values of t. The model 
parameters are arbitrarily chosen as: M = 0.14 MeV and m = 0.33 MeV. 



The functional forms are 



■^1^ "^efil^;', k'-^ ; A, M)'4)es{x, k-^ ; A, A') 



A / dk- 



4CDwHx"M\t) 
(1 -x") 

2{k^ -k'^ +m^)Dw{xM\M^)Dw\x",k"^\t) 



x{l - x)x"{l - x")x'{l - x'){l 
2CD^Jv{x" ,k"^\t) + Dw{xM\M^) 



C) 



where we have defined the effective wavefunction 



Afi{x,k ;A,A') 



— x) 



k-x dx.y - Am 5a,-A' 



Dy^'ix,k^\Mi 



(4.58) 
(4.59) 

(4.60) 
(4.61) 



and made the abbreviation A = g'^\C\'^ Nc'm'^ /2{2tt)'^ . It is sensible to think of Eq. ()4.61|) as 
an effective wavefunction since x — > 1 — x symmetry has been lost. Moreover the wavefunc- 
tion vanishes at x = 1 and is non-vanishing at x = 0. This is a desirable addition to the 
dynamics stemming from the regularization. Notice the ladder kernel (|4.51|) is momentum 
independent and hence does not vanish at x = 0, 1. This is the dynamical reason why the 
un-regularized wavefunction Eq. (|4.54|) does not vanish at the end points. Continuity of the 
GPD at X = C follows directly from Eqs. (I4.58H4.6()I) . 

In Figure lOl we plot the GPD for the parameters: M = 0.14 MeV and m = 0.33 MeV. 
On the left, the graph shows the GPD for a few values of ^ as a function of x for fixed t, 
while on the right we have fixed ^ and t varying. 



4.3.3 Sum rule and polynomiality 

Since not manifest, one should check the covariance of the model Eq. (|4.57() . With the 
covariant starting point Eq. ()4.56|1 . we anticipate polynomiality will be satisfied which pro- 
vides a useful check on our expressions Eqs. H4.58fH?Bm) . First we define the moments of 
the GPD with respect to asymmetric variables 

PniCt) = I x"//(x,C,t)dx. (4.62) 
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Figure 4.7: Polynomiality conditions checked for the GPD Eq. (|4.57j) . The moments 
Pn{C,t = 0) from Eq. (|4.62j) are plotted as a function of Q for n = 0, 1,2,3. Additionally 
X's denote the simple (n + l)-point polynomial fit to the moment Pn- 



Polynomiality requires the moments P„ to be of the form 

n 

Pn{C,t)=Y,(^nj{t)e- (4.63) 
j=0 

The zeroth moment is merely the sum rule for the form factor, hence aoo(0 = F'it)- For 
simplicity, we check a few of the lowest moments for the polynomiality condition Eq. (|4.63|) 
at t = 0. In Figure we plot the moments P„(C)0) for n = 0,1,2,3 which appear as 
smooth functions. Additionally we plot simple (n+l)-point polynomial fits to the moments, 
which line up nicely with the integrals 1)4. 62() . 

4.3.4 Double distribution 

To construct the DD, we shall first proceed incorrectly by appealing to the Lorentz invari- 
ance of the form factor as in Section f4.2.11 This will at least lead to one component of the 
DD satisfying the reduction relations and can be compared to Eq. ()4.57|) . 

Using Eq. (|4.58j) in the A"*" = (Drell-Yan) frame, we can write F{t) in the form (|4.18|) 
with 



F{x,y;t) = [3m^-M^x{l-x)+y{l-x-y)t 



2ttA e{l-x- y) y{l-x- y) 



r 1 ^ ' 

(1 — x) -m? — M^x(l — x) — y{l — x — y)t 

(4.64) 

Aside from factors arising from spin, this DD is basically the same as that considered [71] 
which one can realize by utilizing = — M^/4 + m? . Not surprisingly, then, this DD 
satisfies Eqs. H4.19til?^T|) . For reference we give the quark distribution function 

(1 - xf 3m^ - Mlx{l - x) 



= Tr2rT^ V13 ' ("^-65) 

6 \ra^ — M^x\ \ — x)Y 

which could be calculated directly from -(/'eff in Ec[- (|4.61|) . 

In FigureEiHl we plot the GPD Eq. (|137)) as weh as the GPD derived from DD Eq. 
via Eq. (|4.23)) . As in Section [4.2. 11 the two are different despite the fact both models are 
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Figure 4.8: Comparison of covariant GPDs for the spinor triangle diagram. The GPDs 
Eq. (|4.57|) (denoted LC) and Eq. (|4.(i4|) (DD-based) are plotted as a function of x for fixed 
C, = 0.9 and t = —Am? for the mass = 0.15m. We also plot the difference between the 
two curves {5). The area under the curves is identically F{—Am?) for LC and DD-based 
GPDs, and hence zero for their difference 5. 
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Figure 4.9: Comparison of GPDs: GPD calculated from the DD Eq. (|4.64|) via Eq. (|4.23j) 
(DD-based) compared with the light-cone projection of the form factor Ea. ()4.57|) (Light- 
cone) for fixed Q = 0.4 at t = 0. Here we plot the ratio i?(x, Q) appearing in Eq. (|4.67j) as a 
function oi x > C,. Positivity constrains this ratio to be less than one. 



covariant and posses the same form factor and quark distribution. Additionally we plot 
their difference function of x. 

4.3.5 Positivity constraints 

Here we demonstrate another difference between the GPD in 1)4. 57|) and the one stemming 
from the one-component DD Eq. (|4.(i4j) . To do so, we look at the positivity constraints. 
Originally these constraints appeared in [70, 101] and were derived from the positivity of the 
density matrix by restricting the final-state parton to have positive plus-momentum [and 
ignoring the contribution from E{x,(,t) for the spin-^ case]. Although the matrix elements 
involved for GPDs are off diagonal, they are still restricted by positivity and their diagonal 
elements. Correcting the constraints for the presence of the ^'-distribution was first done 
in [26] . By considering the positivity of the norm on Hilbert space, stricter constraints for 
the spin-| distributions H and E have recently appeared as well as constraints for the full 
set of twist-two GPDs [102]. 

For the scalar pion case there is of course no contribution from the non-existent E- 
distribution and hence the original bounds are actually correct (modulo factors due to the 
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difference of a spin- 2 proton versus a spin-0 pion). Given the matrix element definition of 
the GPD consistent with equation ()4.1UI) . namely 



H{xX,t) 



1 f dz- 




(7r(P')lV5(0)7+^(^")K(^)) 



(4.66) 



l-C/27 47r 



the spin-0 positivity constraint (for x > C,) reads 



i?(x,C)^ (l-C/2) 



i/(x,C,0) 



< 1 



(4.67) 



where q{x) is the model distribution function in Eq. (|4.65|) . Of course the above result 
holds for finite —t, however since the function decreases with —t, Eq. (|4.67|) is the tightest 
constraint. Notice for 7^ 0, the limit t = is in an unphysical region. If we treat this limit 
as formal, however, and analytically continue our expressions, we can use Eq. ()4.67|) . Such 
continuation is consistent with the light-cone Fock space representation of GPDs [26,27]. 

Given the constraint Eq. (|4.67j) . we can test whether GPDs calculated from the light- 
cone projection ()4.58|) and DD ()4.64|1 satisfy positivity. In Figure we plot R{x,C) for 
each GPD as a function of x for the fixed value of C = 0.4. There is noticeably different 
behavior in the figure: positivity is violated by the DD-based model. As above (Section 
I4.2.2() . we must carefully derive contributions from the other component G{x,y;t). 

4.3.6 Derivation of the correct DDs 

To derive both F and G DDs for the spin-^ pion model, we must consider the action of 

the operator ^^^^iD^^ ■ ■ -iD^"^ between non-diagonal pion states. Inserted into Eq. (|4.56j) 
which is now taken in the symmetrical frame, this operator produces 



complicates evaluating Eq. (|4.68|) by requiring contributions from diagrams reduced by one 
propagator. Since we have smeared the photon via (|4.55|) . the reduced diagrams are finite. 
Let us denote the propagators simply by 21= (fc — P+A/2)^ — m^+ze, *8 = {k+A)'^ — m'^+ie 
and C = /c^ — + ie. To correctly evaluate Eq. (|4.68|) . we must write the trace as 



and evaluate each term separately canceling propagators in the denominator of Eq. (|4.68|) . 
These cancellations enable the DD to be read off by inspection. 

These integrals can easily be evaluated using Feynman parameters. For example, let us 
consider the non-reduced contribution from Eq. (|4.70j) . The denominator of Ea. l4.68l aDPears 
as 21 and so we introduce two Feynman parameters {x, y} to render the denominator 

specifically in the form [a;2l + y^B -|- (1 — x — y)(t]~^. One then translates A;^ to render the 




Tr{ . . .j'^ = 4 [P'' (| - ^ (53 + £)) + ^ (M^ - I - 21 + ^ - e:)) + A:^(m2 - I - 21) 1 (4.70) 
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integral (hyper-) spherically symmetric via the definition = l^^ + f3P^ — (a + l)A^/2. 
Here (3 = x and a = x + 2y — 1. Using a Wick rotation to evaluate the resulting integral 
over I, we can cast the contribution to Eq. H4.68() from non-reduced terms in the form 



ttA 



dp 



1-/9 r _ ^{^ 

da[{l - Pf - a^]D{P, a; tf 2P{PM^ + (1 - p)t/2) - Aa{M^ - t/2) 



1+/3 



E 

fc=0 



k\{n-k)\ 



A\ A*"} 



(4.71) 



with D{P, a; t) given by Eq. (|4.26|) . Given the form of Eq. (|4.71|) . we can identify {/?, a} as 
DD variables and hence read off contributions to F and G DDs. 

6F{P, a; t) = ttA [(1 - pf - a^] D{P, a; tf {pM^ + (1 - P)t/2) (4.72) 
6G{p, a; t) =7:Aa [(1 - pf - a^] D{p, a; tf {Ml - t/2) . (4.73) 

Notice that these contributions respect the properties of DDs, namely 5F is even in a while 
6G is odd. This need not be the case, however, for each intermediate step of the calculation, 
e.g. contributions from OS-reduced terms and C-reduced terms are individually neither even 
nor odd in a while their sum is even and difference is odd. 

Ignoring for the moment contributions from 2l-reduced terms, we arrive at the DDs 



F{P, a;t)=T: A D{P, a; tf (1 - p)m'^ - Pa^Ml + {I - p)[{l - pf - a^]t/A 



G{P, a;t) = -TT Aa D{P, a; tf 



m 



Ml{l -P-a^) + [(1 - pf - o?\tl\ 



(4.74) 
(4.75) 



The contribution from 2l-reduced terms has the form of a D-term in that it is proportional 
to 5(/3). Using Feynman parameters for the denominator iB^ ^ and suitable changes of 
variables, we arrive at the contribution to Eq. (|4.68)1 



1 , 7r^a"+i(l -a2) 
da-, 



Li K - (1 - a2)t/4]2 
from which we can identify the D-term 

a(l — o? 



A\ Ml 

2"/ 



A\ Mn} 
'2) 



D{a;t) = ttA- 



m? - (1 - a2)t/4]2 



(4.76) 



(4.77) 



Although strictly speaking, the D-term is a contribution to the G-DD, we shall treat it sep- 
arately for ease. Furthermore, the complete D-term arises from the /3-integral of G{P,a;t), 
see Appendix IHI and [75]. 

Switching now to asymmetric variables, we have 



F{x,y;t) 
G{x,y;t) 
D{y;t) 



TT A D{x,y;tf 

X (1 - x)m^ - x{x + 2y - ifU^ + (1 - x)y{l - x - y)t 
-■K A{x + 2y- 1) D{x, y; tf 



m 



y{l - X - y)t - Ml{l 



{x + 2y- 1) 



7r^y(l - y){2y - 1) [m^ - y{l - y)t\ 



(4.78) 

(4.79) 
(4.80) 
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Figure 4.10: Contributions to the covariant GPD from DDs. The hght-cone GPD Eq. 1)4. 57() 
and the DD-based Eq. I|4.81j) are identical, denoted (LC) and plotted as a function of x for 
fixed C = 0.9 and t = —Am? for the mass = 0.15m. We also plot the individual 
contributions from F in Eq. (|4.78j) and from 5^^+ 2(2-^) ^ ™ Eas. l4.7^ and l4.8ni denoted 
(G + D-term). 



where the function D{x,y;t) is given by Eq. (|4.29() . Accordingly F{x,y;t) is Miinchen 
symmetric and G{x,y;t) is antisymmetric, while D{y;t) is antisymmetric about y = 1/2. 
Notice F{x,y;t) in Eq. 1)4. 78() is not that of Eq. ()4.64|) . Given the ambiguity inherent in 
defining F versus G DDs [cf. Eq. ()4.8)) ]. there is no reason to believe the Fs would be the 
same. In principle, we could construct a quasi-gauge transformation [75] to render the Fs 
the same. This would enable identification of the missing G function unique to Section 
14.3.41 We shall not pursue this tangential point, ^ but simply point out that the procedure 
in [71] leads only to one component of the DD and precisely which component [in the 
decomposition of Eq. 1)4. 1|) ] is unknown and hence useless. 

The function F{x, y; t) satisfies the reduction relations: it reduces to the quark distri- 
bution via Eq. 1)4.20(1 and integrates to the form factor Eq. (|4.21|) — the latter can only be 
checked numerically. Lastly then it remains to see whether the DD-based GPD lines up with 
true GPD calculated in Section 14.3.21 To construct the GPD we use the form of Eq. ()4.32j) 
modified to handle the extra D-term in Eq. ((4. 80(1 separately 



Hix,C,t)= I dz I dy 5{x-z-Cy) F{x,y;t) + j^G{x,y;t) 




''"^lD(x"-t) 



2(2 - C) 

(4.81) 

In Figure KWl we plot the GPD Eq. (|137|) and the DD-based Eq. (|OT|) . They are 
identical. We also plot the individual F and G + D contributions to the GPD. Even in the 
absence of the D-term, the contribution from G cannot be neglected in ascertaining the DD. 
Again calculation of the GPD from the full two-component DD agrees with the light-cone 
GPD definition. The argument of Lorentz invariance used to calculate an i^-type DD that 
satisfies the reduction relations does not determine the GPD. 

For a given model, there exists a unique two-component DD in the decomposition of 
Eq. (|4.1|1 . This DD produces a GPD which is then consistent with the light-cone projection 
and satisfies all the necessary reduction relations. On the other hand, if one is interested in 
constructing different GPD models satisfying the reduction relations, polynomiality, etc., 



Notice the contribution to the GPD from our D-term resembles that of iifinst in Eq. l|4:59t but is not 
identical. Both terms originate from a reduction of the spectator's propagator. In the case of the D-term, the 
spectator's propagator is completely removed by the 2l-reduction. For //inst, there is residual a;-dependence 
stemming from the light-cone instantaneous propagator 7^/2P"''(l — x). 
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one could exploit the ambiguity of one-component DDs. Consider the following. The 
contribution to the GPD from F in Figure B.IUI is markedly different from Eq. H4.57() and 
the analogous contribution from Eq. ()4.64|1 . These three GPDs, however, share the same 
form factor and quark distribution, since their respective G functions do not contribute to 
the reduction relations. Thus one can always use DD-gauge freedom [75] to transform a 
given F into a new function and throw away contributions from G. The result is an infinite 
set of different GPDs with identical form factors and quark distributions. As pointed out 
in Section [4.3.51 one must be careful to maintain positivity although it is likely that there 
still is an infinite set of GPDs which would. 

4.4 Summary 

Above we consider two covariant models for pious: one with scalar constituents and the 
other with spin-^. The spin-^ model requires regularization and we choose the method 
in [99] in order to maintain positivity. For each case we derive the GPD from its matrix 
element definition which forces us to consider the triangle diagram for the form factor with 
the plus momentum of the struck quark kept fixed in a general A"*" ^ frame. 

We also construct the DDs for each model. The approach of [71] leads only to one 
component of the DD (the "forward- visible" piece) which is itself ambiguous, thus result- 
ing GPDs are incorrect and need not satisfy positivity. This fact remains true even for 
C-odd distributions. This means the reduction relations alone are not enough from which 
to calculate a model's DD (and hence GPD). To obtain both components of the DD unam- 
biguously, we calculate the matrix elements of twist-two operators. The resulting DD-based 
GPDs then agree with those calculated on the light cone. The "gauge freedom" inherent 
in defining F vs. G DDs could be exploited, however, for phenomenological studies where 
one is interested in a mathematical fit to data, rather than a calculation of a given model's 
DD. Having investigated the calculation of DDs for pious, we can now turn to the case of 
the proton. 
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Chapter 5 

Double Distributions for the 
Proton 

In this Chapter, we extend the simple scenarios so far pursued to build a DD model for the 
proton. We treat the proton as a bound state of a residual quark and two quarks strongly 
coupled in both the scalar and axial-vector diquark channels. The resulting light-cone 
wavefunction of the proton has appropriate spin structure: containing correlations where 
the residual quark carries the spin of the proton as well as correlations where the quark and 
diquark are in a relative p-wave. We use this model to derive DDs for the proton. Inclusion 
of the spin structure into double distributions is crucial if one wishes to make contact with 
the spin sum rule for hadrons. While this (in essence two-body) model for the proton is 
crude, model parameters can be tuned so that the electromagnetic form factors are fit at 
small momentum transfer. Quark distributions can be matched and resulting GPDs are 
suitable for phenomenology. 

This Chapter has the following organization. First we review our conventions for DDs 
for spin-^ particles and the relation of these DDs to GPDs in Section [5.11 In order to deal 
with Z?-terms efficiently, we opt to calculate three DD functions. The ambiguity inherent 
in defining DDs for spin-^ systems is detailed in Appendix ^ Next in Section 15.21 we 
present the quark-diquark model under consideration. DDs are calculated in this model 
in both the scalar and axial-vector diquark channels. Relevant identities are gathered in 
Appendix ^ while the details of the derivation appear in Appendix EI Intuition about 
the model and its construction is provided in Appendix [jj where the effective light-cone 
wavefunction is extracted from projecting onto the light-cone. Section [5 . 31 presents potential 
phenomenological uses for the model. The model is tuned to reproduce the Dirac and Pauli 
form factors of the proton for small momentum transfer. Resulting quark distributions and 
GPDs are sensible at a low scale but do not contain the relevant non-perturbative small-x 
behavior relevant for GPD phenomenology. One could then follow [71] and augment the 
DDs with empirically parametrized quark distributions and re-tune the model's parameters 
or use perturbative evolution. Consequently these simple model GPDs then satisfy all 
known constraints, including positivity. 

5.1 Definitions 

To begin, we set forth our conventions for DDs and their relation to GPDs. Moments of 
DDs appear naturally in the decomposition of twist-two operators' matrix elements that 
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are non-diagonal in momentum space; moreover, they provide an elegant explanation of the 
polynomiality property of GPDs. To make any progress in calculating DDs, we must use 
the parton model simplification for the gauge covariant derivative: D'^ — > d^. 

The non-diagonal proton matrix elements of twist-two operators can be decomposed in 
a fully Lorentz covariant fashion in terms of various twist- two form factors A^k^t), B^kit) 
and Cnk{t), namely^ 

{P', A'|V^(0)7^^iB^i • • • iD^"^ij{0)\P, A) 



2 



- ^A^j^-^iP) E ^^c„.«)P''. . . . p--. (--) ...(-- 

(5.1) 

where the action of ^ • • • ^ on Lorentz indices produces the symmetric traceless part of the 
tensor and P is defined to be the average momentum between the initial and final states 
= ^{P' + P)^ and A is the momentum transfer, A'^ = {P' — P)'^. T-invariance restricts 
k in the first two sums to be even and odd in the last sum. There are three Dirac structures 
in the above decomposition since in general the twist- two currents are not conserved thus 
allowing for a structure proportional to A^'. The ambiguity of DDs for spin-| systems is 
addressed in Appendix ^ Such difficulties in constructing DDs have been addressed in 
the literature [73-76], and we find the construction in Eq. (|5.1j) the easiest to work with in 
actual calculations. 

The above decomposition can be used to define three double distributions as generating 
functions for the twist-two form factors 

'1 



Ank{t)= [ d(3 [ daP''-''a''F{l3,a;t) (5.2) 
J-i J-i+\P\ 



-1+1/31 

Bnk{t)= d(3 da(3''~^a^K{(3,a;t) (5.3) 

J-l J-1+\I3\ 

Cnk{t)= dp da(5''-^a^G{P,a;t). (5.4) 

As a consequence of the restriction on k in the sums, the functions F{I3, a; t) and K{I3, a; t) 
are even in a, while G{f3, a; t) is odd. The F{f3, a; t) and K{I3, a; t) DDs are similar in 
form to the functions originally employed in [69]. The difference is due to the third DD, 
G{f3,a;t), which incorporates the D-term [74] among other things. 

Summing up the moments in Eq. 1)5. If) , these DD functions then appear in the decom- 

^Beyond the parton model, this decomposition still holds. The twist-two form factors, however, ac- 
quire scale dependence due to the renormalization of the twist-two operators. This scale dependence must 
ultimately be dealt with for phenomenological applications. 
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position of the light-like separated quark bilinear operator 



(P',AV^(-z/2)^, {z/2)\P,\) 



dp 



da e-*/3-P-+^"^-/2 



1+1/31 



icT^^ z A A • z 

a; t) + -^^K,{P, a; t) - -^G,{(3, a; t) 



AM 



ux{P), (5.5) 



where we have appended a flavor subscript q in the relevant places and 2^ is a lightlike 
vector, z'^ = 0. 

Now we define the light-cone correlation function by Fourier transforming with respect 
to the light-cone separation z^ 



1 



dz-e'^""' 



'{P',X'\^g{-z-/2)-f+iP,{z-/2) \P,X). (5.6) 



Above the variable ^, or skewness parameter, is defined by ^ = — A+/2P+. As is customarily 
done, we assume without loss of generality that ^ > 0. The correlation function in Eq. (|5.6|) 
can be written in terms of the two independent GPDs Hq{x,^,t) and Eq{x,$^,t) as 



2P 



-uAP') 



^+Hqix,tt) + ^^Eqix,tt) 



nxiP). (5.7) 



Unlike the DDs, these GPDs are quantities that enter directly into the amplitude for deeply 
virtual Compton scattering (DVCS), for example. Inserting the DD decomposition Eq. (|5.5|) 
into the correlator in Eq. 1)5. 7() . we can express the GPDs as projections of the DDs 

/■I n-m 

Hgix,tt)= d/3 da6{x-/3-Ca)[Fg{P,a;t)+CGgiP,a;t)] (5.8) 



Egix,^,t) 



dp 



-1 



da 5{x- P- ia) [Kg{p, a; t) + ^Gg{P, a; t)] (5.9) 



-n 



from which we can view the ^-dependence of GPDs as arising from different slices of the 
underlying Lorentz invariant DDs. Due to the symmetry of the DDs with respect to a, the 
GPDs Hq{x,(,,t) and Eq{x,(,,t) are both even functions of the skewness parameter ^. 

The GPD Hq{x,^,t) has an important reduction property. Taking the diagonal limit of 
the light-cone correlator Eq. (|5.7jl . we have 



fq{x) = Hq{x,^,Q) 



l~\x\ 
■l+\x\ 



da Fq{x, a; t) 



(5.10) 



In DVCS, the relevant current operators produce the charge and flavor structure 
since there are two photons and thus the charge squared weighted GPDs enter in relevant 
physical amplitudes. To discover the relation of GPDs to electromagnetic form factors it is 
advantageous to consider instead the single photon structure Cq and define 



H{x,i,t) = Y,eqHq{x,i,t) 

1 

E{x,^,t) = Y,e,Eq{x,^,t). 



(5.11) 
(5.12) 
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Since Gq{P,a;t) is an odd function of a, we have J^^dp J^i]^]^p^daGq{P,a;t) = and 
consequently the sum rules 



/ dxH{x,^,t) = y^ eg [ dp [ da Fg{p,a;t) = Fi{t) (5.13) 
J-i ^ J-i J-1+\I3\ 

rl rl .1-1/31 

/ dxE{x,^,t) = } eg dp da Kg{p,a;t) = F2{t) (5.14) 



which relate the zeroth moments to Dirac and Pauli form factors. Higher x-moments of 
Hg(x,£,,t) and Eg(x,(,,t) are even polynomials in ^ as can be seen directly from Eqs. ()5.8() 
and (|5T9|) . 

The GPDs satisfy further constraints: the positivity bounds [26, 70, 101-105] . These 
bounds are particularly important for comparing with experiment. Model GPDs which 
reduce to the experimental quark distribution in Eq. (|5.10j) but violate the positivity bounds 
are not worth considering because one knows from the start that rate estimates predicted 
by the model GPDs are automatically wrong. Usually violation of the positivity bounds is a 
signal that the model is inconsistent with the underlying field theory. Perhaps surprisingly 
violation occurs frequently in more-or-less all standard hadronic models. Of interest to us 
are the basic bounds for both the spin-flip and non-flip amplitudes 



e(x-0|K'^'(^'e,t)| <^/,(^)/,(f^), (5.15) 

which we use as a stipulation in constructing our model. 

Lastly we need to address the negative range of the DD variable p. Experimentally and 
diagrammatically P is strictly positive and crossing symmetry can be used to extend the 
range to positive and negative P, or in our case reduce the range to strictly positive. To 
this end, we define two functions for each DD 

F^{P,a;t) = Fg{P,a;t)±Fg{P,a;t) (5.16) 
K^{P,a;t) = Kg{P,a;t)±Kg{P,a;t) (5.17) 
G^iP,a;t) = Gg{P,a;t)±Gg{p,a;t), (5.18) 

where the antiquark contributions are defined by crossing 

Fg{P,a;t) = -Fg{-p,a-t) 
Kg{P,a;t) = -Kg{-P,a;t) 
Gg{P,a;t) = -Gg{-p,a;t). 

Thus the plus DDs [F^{P,a;t), K^{P,a;t), and Gg{P,a;t)] are odd functions of P and 
the minus DDs are even functions. In partonic language, the minus DDs correspond to a 
difference in quark and antiquark DDs (the valence configuration) while the plus DDs are 
a sum of quark and anitquark DDs. 

By virtue of the above definitions Eqs. (|5.16j) - l|5.18|) . we can remove the explicit negative- 
P parts from DDs and consequently the GPDs. We do so by defining 

H^{x,i,t) = Hg{x,tt)±Hg{x,^,t) (5.19) 
E^ix,^,t) = Egix,^,t)±Eg{x,tt), (5.20) 
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where the antiquark contributions are defined by the crossing analogous to the DDs above. 
In this form, we can rewrite the reduction relations in a more familiar way 



f^{x) = H^{x,0,0) = f \aF^{x,a;t), (5.21) 

J~l+X 



where fg (x) is the valence quark distribution, namely fg (x) = fq(x) — fq{x). The sum 
rules are thus 

/ dxH~{x,i,t) = y2^q dp daFg{P,a;t)=Fi{t), (5.22) 

Jo g Jo j-l+f3 

I dxE-{x,^,t) = y2(iq dp daKg{P,a;t)=F2it). (5.23) 

^0 q Jo J-1+/3 

In our simple valence model for the proton, we only address quark configurations. The 
double distribution variable /? as well as the momentum fraction x are positive below. The 
above positivity bounds in Eq. (|5.15j) also hold for valence and plus amplitudes. 

5.2 Model double distribution 

To calculate DDs for the proton, as a first step we use only a simple model consisting of 
two quarks strongly coupled in the scalar and axial-vector diquark channels along with 
a residual quark. This model can be considered as loosely based on relativistic quark 
models [106] or on the Nambu-Jona Lasinio model of the proton in the static approximation, 
see, e.g., [107,108]. We keep full Lorentz covariance in order to preserve the polynomiality 
of the moments Eq. (|5.1jl . Without covariance, we would not even be able to deduce the 
DDs. 

In order to further simplify this calculation, we mandate only tractable one loop con- 
tributions. To this end, we treat the scalar and axial-vector diquark T-matrices as free 
particle propagators, namely of the forms 

fe^ — m^jj + le 

^'^^^^ = 12 i-— - (5-25) 

fc^ — + le 

respectively. For simplicity we have neglected the term proportional to k^k'^ in the massive 
vector propagator. Evaluation of the contributions generated by this term are considerably 
complicated and will be the subject of a future publication. The proton Bethe-Salpeter 
vertex for our model is thus 

r(A;, P) = ^ X^''^D{P -k)^ rW (fc, P) + ^ x2 D'^'^iP -k)^ [k, P), (5.26) 
where the diquark vertices are direct products of spin and isospin factors 

X^'^ = ^iil5C)^^{iT2) (5.27) 

xS = ^(i7MC)®i=(zr,r2), (5.28) 
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Figure 5.1: Impulse approximation to the twist-two matrix elements. Here the twist-two 
operators with momentum insertion are denoted by a cross. The diquark spectator is 
depicted by a double line and the initial- (final-) state proton has momentum P {P')- 



and we do not append propagators for the first and second quarks. For simplicity we choose 
the quark-diquark vertex functions to be point-like, namely 

r(^) = 101 (5.29) 
rg=757.®ri. (5.30) 

This choice corresponds to modeling only a subset of the possible structures for the proton 
wavefunction, see [109, 110] for a complete discussion. The vertex function also contains 
an overall color anti-symmetrization which we suppress throughout. The conjugate vertex 
T{k, P) = Cr{—k, —P)'^C^. One could modify the point-like vertex with a form factor as 
is commonly done, however, this generally violates the positivity bounds. 

In this model, the axial diquark channel does not contribute to the proton's electromag- 
netic form factors. Thus we can only determine the parameters m and msn by fits to the 
Dirac and Pauli form factors. The parameter niAD could be tuned by fitting the quark distri- 
butions at some scale, however, we shall pursue a simpler course and set mso = "mAD = n^D- 
Alternately ruAD could be tuned from neutron form factor data. 

To derive the F{P,a;t), K{P,a;t) and G{P,a;t) DDs, we consider the action of the 
twist-two operator ^mmi-'-Mh between non-diagonal proton states. Working in the impulse 
approximation (see Figure we have the contribution in the scalar diquark channel, for 
example, 

^4^ ux'{P'){^ + A + m)T>^''^-''-{it + m)nA(P) 

[ik2-m2 + ie]«[(/fc + A)2-m2 + ie]«[(fc-P + A/2)2-m|) + ie]' ' 

where we have chosen an analytic regularization (so that a > 1) respecting the positivity 
bounds. 2 Since the NJL model is non-renormalizable, the choice of scheme is incorporated 



^Alternate schemes using Pauli- Villars subtractions often regulate such models and are also attractive 
from the perspective of DDs since Lorentz covariance is maintained. These subtractions, however, generally 
violate the bounds in Eq. (15.151 . For example, in the NJL model of the pion with two subtractions [79], the 
positivity bounds, which were ignored by the authors, are violated for small values of —t. For the case of 
a quark-diquark model regularized via Pauli- Villars subtractions, the violations are more severe and persist 
for all values of —t due to the mismatch of end-point and crossover behavior. This commonly encountered 
problem is discussed in [63]. For these reasons, we have used the analytic regularization above, which we 
also employed previously for the pion [73] in Chapter 2] 
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into the dynamics and hence the choice of regularization should maintain desired properties. 
For phenomenological estimates of GPDs, the positivity bounds are of crucial importance, 
and our regularization choice respects these bounds, see Appendix |j] for details. One can 
view the regularization a > 1 as mimicking the non-local, non-perturbative structure of the 
twist-two, quark-antiquark vertex. This choice of smearing maintains current conservation 
but does not respect the Ward-Takahashi identity. Thus the normalization of amplitudes in 
Eq. (|5.26j) is only approximately preserved. To remedy this feature, we add an a-dependent 
prefactor to the axial diquark contribution. This factor is then adjusted so that A^^ = 1 
and consequently the u-quark distribution is correctly normalized, = 2. An additional 
drawback of the regularization scheme is that the induced quark form factors do not become 
point-like for large momenta. While this is inconsistent with asymptotic freedom, the model 
is meant only for use in the low momentum transfer region because the diquark substructure 
cannot be resolved. 

The symmetric, traceless tensor T in Eq. ()5.31() is 

pMMi...Mn ^ ^{/.(^ ^ ^/2)^i ■■■ik + A/2)^"> (5.32) 

The involved numerator structure complicates calculation of the symmetric traceless part 
of the tensor, however, this can be calculated directly without recourse to explicitly writing 
out such tensors of rank n. See AppendixEI In order to compactly write out the DDs, we 
define the auxiliary functions 

DoiP, a; t) = pmj) + (1 - /3)m^ - - (3)M^ + [(1 - 0f - a^]t/4, (5.33) 

which is the typical energy denominator in both channels and 

D(/3, a; t) = 3N^^^^ [(1 - pf - a'] DM a; (5.34) 

which is the typical prefactor for all DDs in this model. 
Calculation of DDs for the scalar diquark yields 



(m + /3M)2+[(l-/3)2-«2]| + zZ^^ 

6g^uD{P,a;t) \ 2M{1 - (3){m + (5M) ]■ (5.35) 

4Ma{m + (3M) 



In the axial-vector diquark channel, we find 



Kjj''\(3,a-t) I =lz{a){26,,d + Ku)Dif3,a;t) 



X 



'{2m + (3Mf - 3m2 - [(1 + -a']{ + 

2Mp[2m - (1 - (3)M] | , (5.36) 

-4Ma{2m + /3M) 



where Z{a) is the regularization dependent factor akin to wavefunction renormalization. As 
commented above, the value of Z(a) is chosen to preserve = 1. 
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Figure 5.2: Comparison of fits for GEit) and GM{t) with empirical parameterizations. The 
ratios GE{t)/GD{t) and GM{t) / fipGoit) are plotted vs. -t in GeV^. The curves "Arr:l" 
and "Arr:2" correspond to the parameterizations of Gsit) given in Tables I and II of [113], 
respectively. The curves a = 2 and a = 3 are our model fits for different regularizations. 



5.3 Phenomenological applications 

Our philosophy is to tune the parameters m and mo so that proton electromagnetic form 
factor data at low momentum transfer are reproduced. This is particularly simple, since 
the axial diquark contribution cancels out of the the proton form factor. We choose to 
fit to the charge radius and magnetic moment and the experimental values we use are 
< r"^ >= (0.870 fm)^ and fip = 2.79fi]\j. For the electromagnetic form factors of the 
proton, there is high precision data from JLAB [111,112] and a recent global analysis and 
parametrization of [113]. The Sachs form factors are known experimentally to about 2% 
accuracy in the small momentum transfer regime and are given by 

GEit)=F^it) + ^F2it) (5.37) 
GM{t)=Fi{t) + F2{t). (5.38) 

We treat the regulator parameter a rather unconventionally; we use it as a means to 
explore different covariant forms for the wavefunction. For a = 2, magnetic moment and 
charge radius are reproduced < 0.1% for m = 0.437 GeV and mo = 0.726 GeV. The Sachs 
electric and magnetic form factors match up < 5% for —t < 0.2 GeV^. For a = 3, magnetic 
moment and charge radius are reproduced < 0.1% for m = 0.565 GeV and mo = 0.825 GeV. 
For these parameters, the Sachs electric and magnetic form factors also match up < 5% 
for -t < 0.2 GeV2. In the fi gures, we compare the phenomenological form factor fits to the 
two parameterizations of [113]. The fits for GEif) and GMit) are plotted in Fig. 15.21 As is 
standard, we plot ratios of electric and magnetic form factors to the empirical dipole form 
factor, namely 

Gn{t)={l-^J\ (5.39) 

where the dipole mass squared is M|, = 0.71 GeV^. 

We can also determine the u and d quark distributions in our model. Since we do not 
have antiquarks, the plus and minus distributions are identical f^{x) = fq{x). In Fief. 15.31 
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Figure 5.3: Quark distributions for the proton model. The u- and d-quark distributions are 
plotted as a function of x for the a = 2 and a = 3 fits to the form factors. 



we plot the u and d quark distributions as a function of x for the a = 2 and a = 3 fits. The 
distributions are properly normalized so that = 2 and = 1. As commented above 
this normalization requires a relative a-dependent factor for the axial-diquark contributions. 
Again this is required because the regularization scheme we choose does not preserve the 
Ward-Takahashi identity. Without the extra factor, the violation is ~ 5-10%. Notice the 
distributions do not vanish at the end-point x = 0. This is typical of NJL type calculations. 
The kernel is independent of momentum and hence the wavefunction should be non-zero 
at both end-points. The fact that the distributions vanish at x = 1 is due to our choice of 
regularization. Physically it is thus reasonable to think of the regularization as mimicking 
contributions from higher Fock states. Moreover, the model's scale is a-dependent, cf Figure 
15.31 Additionally in Figs. 153] and EIHl we have plotted the H^{x,^,t) and E^{x,^,t) 
GPDs. The figures show the GPDs at fixed -t = 0.1 GeV^ for a few values of ^. The GPDs 
plotted are for the a = 2 fit; the distributions are qualitatively similar for the a = 3 model 
parameters. 

The distributions plotted in Figs. I5.3H5.5I are presumably at some low hadronic scale 
intrinsic to the model. One way to determine the scales of the a = 2 and a = 3 models is 
to use the evolution equations to evolve empirical parameterizations down to a scale where 
the first few moments of our model distributions agree. This procedure is not unique: many 
models can reproduce the empirical quark distributions at higher scales. Also the use of 
perturbative evolution is questionable at low scales. While the evolution kernels for GPDs 
are known at next-to-leading order [114-118], perturbative evolution cannot generate the 
non-perturbative small-x physics which our model lacks and the small-x physics is crucial 
for relating to DVCS data. The the leading-twist DVCS amplitude is a convolution of GPDs 
with a hard scattering kernel that emphasizes regions where the final-state wavefunction 
is evaluated near the end-point. In fact, the imaginary part of the amplitude is directly 
proportional to an overlap of light-cone wavefunctions, where the final state is at zero x. 

An attractive alternative to using perturbative evolution to define the scale of the model 
is to implant a realistic quark distribution at some scale via factorization of the DD [71]. 
With a proper choice of scale, perturbative evolution can then be used from jj, up to the 
experimental scale. This seems necessary for the parametrization to be most useful for 
experimental comparison. Ultimately the experiment must guide phenomenology in con- 
structing parameterizations that are general enough to account for the observed correlations 
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among GPD variables. While the DD certainly does not factorize, the Ansatz in [71] avoids 
factorizing the variable dependence of the GPDs. Further phenomenological consideration 
of this quark-diquark model is currently underway and will be presented elsewhere. 
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Chapter 6 

Conclusions 



Above we have presented a comprehensive but more-or-less pedagogical treatment of mod- 
ehng GPDs. In the course of attempting to model these distributions, we came across 
subtle issues about the light-front BSE and its reduction scheme. Bent on applications to 
GPDs, we reviewed our work concerning current in the reduction formalism in Chapter El 
There we proved a replacement theorem for apparent non-wavefunction type contributions. 
The resolution of this issue was crucial to enable proper calculation of GPDs in light-front 
perturbation theory. 

In Chapter |31 we detailed our calculation of GPDs in the light-front Wick-Cutkosky 
model. This calculation provided a non-trivial example of the light-cone Fock space rep- 
resentation of DVCS, in which form the GPDs' positivity bounds are manifest. The non- 
vanishing of GPDs at the crossover between kinematic regions was tied to higher Fock 
components. Furthermore we saw how the features of light-cone time-ordered perturbation 
theory entail the continuity of GPDs. Continuity of GPDs ensures the QCD factorization 
theorem for DVCS. The proton case is likely to proceed along the lines outlined for the 
Wick-Cutkosky model, however, incorporation of spin degrees of freedom is difficult due to 
divergences. 

Next in Chapter [IJ we considered an alternate approach to GPDs based on covariant 
DDs. This approach allowed us considerable insight as to the Lorentz covariance properties 
required of GPDs. Unfortunately in this formalism the positivity bounds are not transparent 
but can be used as a guide in model building or calculation of DDs. Using simple models, 
we not only compared covariant and light-front calculations of GPDs, we also corrected the 
inconsistencies of the method proposed in [71] for the calculation of DDs. Understanding 
the DD formalism is pertinent because ad hoc forms for DDs are almost exclusively used to 
model GPDs. We have provided complete calculations of DDs for several models of scalar 
bound states. Such calculations fill a void in the literature. 

Finally in ChapterEl we used a simple two-body model to calculate DDs for the proton. 
This is the first such study. The model consists of two quarks strongly coupled in the scalar 
and axial diquark channels along with a residual quark. The simplicity of such a model 
allows for analytic computation of the DD functions which contain appropriate spin struc- 
ture. The simplicity also allows for us to make contact with the light-cone wavefunctions 
and the equivalent overlap representation of GPDs. This toy model study thus enables a 
comparison between the physical intuition of the light-cone Fock space representation of 
GPDs and the manifestly Lorentz invariant decomposition of DD functions. 

We were careful to choose a regularization scheme that satisfies both Lorentz invariance 
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as well as the positivity bounds required of GPDs. Our model, although toy-like in nature, 
can be used to match the electromagnetic form factors of the proton and can either be 
augmented with realistic valence quark distributions or be evolved to higher scales. 

The scope of our investigation is two fold. One direction is to calculate as many DDs 
in simple scenarios as possible. This will give modelers a better sense of the form and 
behavior of DDs and may assist with empirical parameterizations of GPDs. In another 
nearly orthogonal direction, one can improve upon the proton model used here in order to 
see how various features of a realistic proton wavefunction manifest themselves in processes 
like DVCS. Further extension of the model presented in Chapter |3 and its application to 
future data are under investigation. 
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Appendix A 

Conventional Notations and 
Notational Conventions 



For a general Lorcntz vector a'^ = {a^,a^,a^,a^), we define the light-cone variables as 
the linear combinations 

a± = -^(aO±a^). (A.l) 
Now consider the scalar product of two Lorentz 4- vectors, and b^. 



(A.2) 



This implicitly defines the metric tensor g^j^y, namely gj^ = g |_ = 1, gn = g22 = — 1 and 

all unspecified entries vanish. Notice in our conventions gi^^/ = gi^" . For convenience we 
shall also refer to transverse two- vectors, namely a-*- = (a^,a^). For an on-shell particle of 
mass m, the dispersion relation k"^ = m? takes on a form devoid of the square root 



k- = 



2k+ ' 



(A.3) 



For Dirac gamma matrices, when explicit forms are needed we use the Bjorken and Drell 
convention: 



and 



r 



where i = 1, 2, 3 and the Pauli matrices cr* are given by 
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(A.4) 

(A.5) 

(A.6) 
(A.7) 



92 



and 

One can similarly form plus and minus gamma matrices in the analogous way 

7^ = ^(7° ±7')- (A.9) 
Also useful when considering fermions are the projection operators A± that are defined by 

A± = ^7^7^- (A.IO) 



Notice (7"*")^ = (7 )^ = and the projector is simple to express in terms of 7^^, namely 

V2 



A± = :^7°7^. 
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Appendix B 

Instant and Front Form Dynamics 



In this Appendix we review the features of relativistic dynamics and the Poincare algebra. 
Above in Section 11.2.21 we used a toy model to illustrate the characteristics of both the 
instant and front form of dynamics. Here we give a more theoretical description of the 
differences between the two forms and thereby define what it means to specify a form of 
dynamics. 

A relativistic system is characterized by the ten operators of the Poincare group. Full 
relativistic invariance stems from the algebra formed from the Poincare generators. The 
generators of the Poincare group consist of the four momenta P^, and the six operators 
embedded in the antisymmetric tensor M^'^. The latter are the angular momenta, specified 
by 

V = ^e'^'^M^^; (B.l) 

and the Lorentz boosts, specified by 

= (B.2) 

where the indices i, j, and k run from 1 to 3. 

These ten generators satisfy the Poincare algebra, namely 

i[pA',p^] = 0, (B.3) 
i[M^"',pP] = g'^Ppf' - gf'Pp'', (B.4) 
ilMP" , MP^"] = g''^M''P-g''PM^'' -g^^MPP+g^PMP". (B.5) 

This algebra enforces that the PP are independent and transform as a four-vector under 
the Lorentz transformation, while M^'^ must transform as an antisymmetric tensor. To 
find a representation of the Poincare algebra, one constructs the generators from dynamical 
variables: positions, momenta, spins, etc. 

To make our discussion concrete, we restrict our attention to the case of a single free 
scalar particle. The dynamical variables at our disposal are the position and momentum 
four-vectors xP and pp that obey the quantum mechanical commutator 

i[xP,pP]= gP". (B.6) 

This commutation relation will furnish a representation of the Poincare algebra provided 
that we construct the generators as 

PP = pP (B.7) 
MP" = xPp'' - x^pP. (B.8) 
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We also must impose the relativistic constraint = w? to have Lorentz covariance. 

In order to specify the form of dynamics of the system, one needs to create a fohation of 
space-time into space and time. This is done by choosing the initial surface S on which the 
general time variable r is zero. Choices for S are restricted by the condition that all possible 
world lines may intersect S exactly once. Clearly this is true also in classical mechanics 
where Galilean time t is completely uncoupled from three-dimensional Euclidean space and 
thus the instant t = specifies a unique point on the space-time trajectory of all particles. 
As we know from relativity, the world lines of particles are restricted to lie on or within 
the fore and aft light cones. Because relativity restricts the class of possible world lines, it 
enlarges the possible initial surfaces S and hence broadens our ability to define "time" in 
relativistic dynamics. 

Once we have specified the initial surface S, we have also to construct the normal N^^. 
The significance of the normal is that • P serves as the Hamiltonian in the sense that 
it evolves the system away from the initial surface on parallel slices in r. Thus in close 
analogy with non-relativistic quantum mechanics, the Hamiltonian is dynamical. 

Now having specified the dynamics of the system, we imagine the initial conditions 
(wavef unctions, etc) are defined on the initial surface S. The Poincare generators now 
fall into two categories: kinematical operators and dynamical operators. The kinematical 
operators have the property of mapping S — > S. They are called kinematical because the 
behavior of states under these operators does not require dynamical input since r remains 
zero. On the other hand, operators that do not act transitively on S are dynamical, i.e., 
the effects of these operators on states require knowledge of the system away from S. 

To give a familiar example of a form of dynamics, we consider first the instant form. 
The initial surface is specified by the conventional instant of time S : = 0. The conjugate 
to time is the energy 

= N ■p= A/p2 + m2. (B.9) 

To find out which operators are kinematical and dynamical under this foliation, we use the 
representation in Eq. (|B.7|) evaluated on the surface S. The Hamiltonian, which for the 
instant form is = P^, is automatically dynamical in any system. Its appearance in the 
other Poincare generators indicates that they too are dynamical. Thus we find six kinematic 
operators 

P'=p' 

= x'P^ - X^P' 

and additionally three dynamical ones 

= (B.12) 

This situation of instant form dynamics could be anticipated directly from special relativity. 
We know that spatial translations and rotations leave the surface = unchanged. Hence 
the generators of these transformations must be kinematical. On the other hand, boosts in 
any direction mix the coordinate of that direction with x^, hence dynamics are required in 
order to get the system off the surface S. 

Another example of a form of dynamics is the front form. In this form, the initial surface 
is specified by a plane tangential to the light-cone S : = 0, see Figure Ib!T1 The normal 
to this surface is given by the lightlike vector 

iV'^ = -^(1,0,0,-1) (B.13) 
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(B.IO) 
(B.ll) 



Figure B.l: The hypersurface S : = that defines the front form of dynamics. The 
surface S is the plane tangential to the light-cone, = 0. Figure adapted from [129] 

and hence the Hamiltonian is given by 

N-P = p-=p-. (B.14) 
The kinematic operators in the front form are the transverse and longitudinal momenta 

p^=p^, P+=p+, (B.15) 

an angular momentum 

M^"^ =x^P^ -x^P^, (B.16) 

and the light-front boosts 

= -x-^P+, M+- = x~P+. (B.17) 

There are thus only two additional dynamical operators given by 

M--^ =x~P-^ -x'P~. (B.18) 

The physics of the partition into kinematic and dynamic classes of generators is clear. 
The two dynamical operators involve rotations about x^ and x^ axes. These rotations 
change E and hence require dynamics. There are three independent translations that leave 
the surface x~^ = alone and hence the three kinematical operators P^ and P~^. Rotations 
within the surface S that leave x"*" = unchanged are generated by the angular momentum 
M^^. There are two special combinations of rotations and boosts (from the instant form 
mindset) that make up the kinematical light-front boosts M~^-^. Finally there is the third 
kinematical boost M ^ = —K^ that is identical (up to sign) to the instant form boost in 
the x^ direction. Under such a boost of rapidity r] the coordinates get mapped as follows 



/xO\ 




fx^ 


cosh r] + x^ sinh r]\ 


x' 




x' 


X2 












\x^ 


sinh r] + x^ cosh r] J 



(B.19) 
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Thus the hght-cone time x'^ — ^ e^x'^ is merely rescaled under this boost. The dynamics are 
diagonal with respect to K'^ boosts. Notice, consequently, the initial surface S is indeed 
invariant under the action of M ^ since = on S. 

Of course our discussion of the forms of dynamics above was limited to the case of a free 
scalar particle. Inclusion of dynamics leads to interaction terms in the Poincare generators. 
Finding such terms that are still consistent with the Poincare algebra is the real struggle 
of constructing a relativistic dynamical system [1]. Indeed Poincare invariance alone is 
not enough to guarantee relativistic covariance. Despite this, a few special cases that do 
maintain full relativistic covariance have been arduously worked out, see e.g. [2]. This issue 
at play here is how to impose relativistic covariance on quantum mechanical, fixed particle 
number systems with interaction potentials. At heart these two are incompatible: action 
at a distance and particle conservation have stepped aside in favor relativistic field theory. 
Nonetheless, the discussion in this Appendix motivates the use of front form dynamics and 
the light-cone quantization of field theories. 
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Appendix C 

Light-Cone Quantization 



In Section 11.11 we derived an expression for the DIS structure function /i (x) . Here we 
elucidate the physics of this distribution function using Hght-cone quantization in the hght- 
cone gauge. For reference the leading-twist function encountered in DIS is defined as 

fi{x) = \j ^e^^--(^l?(0)7+^(An)|P). (C.l) 

Notice that the correlator is not an equal-time Green's function since the quark field is 
annihilated at = A/\/2A and created at = 0. During this time, fluctuations of the 
residual proton state occur that deny the matrix element a probabilistic interpretation. On 
the other hand, both quark fields are evaluated at x"*" = which tells us the function 
is an equal light-cone time correlator. To see the density interpretation of and hence 

justify the name PDF, we must accordingly use light-cone quantization. 

On the light cone, spinor fields have components that are dynamically constrained. To 
separate out the dependent fields we write 

V' = V+ + V'-, (C.2) 

where V'+ is the good component of the spinor field, and ijj- is the bad component that is 
dynamically determined from the QCD equations of motion. Explicitly the good and bad 
components are defined from the projectors (given in Appendix^, namely 

V'i = AiV'- (C.3) 
The dynamically independent good fields -0+ are then quantized at equal light-cone time 



= -1={A+U6{x- - y~)<5(x^ - y^), (C.4) 

x+=y+ 

with all other possible anticommutators (among the good fields) vanishing. This anticom- 
mutation relation is naturally solved by a Fourier decomposition in terms of plane-wave 
modes 



^+(--,x^) = / ^?^^(^+) E h(A:Mc^)^+A(A;)e-(-'=^-"-''") 



dk+dk^ 

A=± 

+ 4(A:+,k^)r;+A(A:)e^(""'=''-"^-'^^)l . (C.5) 
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The above plane- wave decomposition is merely a projection of the four-dimensional solutions 
to the Dirac equation onto the hypersurface = 0. The spinors in the above decomposition 
are good projections of the light-cone spinors in Appendix (jj namely u-^-x(k) = A^ux{k) 
and similarly for v^x{k). 

The quark mode operators b and 6^ must satisfy the anticommutation relation 

[b{{k+,k^),by{k'+,k'^)} = 2k+6x,x'S{k+ - k'+)6{k^ - k'^), (C.6) 

and analogous relations hold for the antiquark creation and annihilation operators S and 
d. Anticommutators between all remaining pairs of quark and antiquark mode operators 
vanish. 

Using the relations in Appendix we can convert the bilocal operator in Eq. H1.19() to 
a more useful form in terms of good fields 

V^(0)7+V'(An) = ^/2^;^(0)V'+(An). (C.7) 

Thus inserting the mode expansion Eq. HC.5() for the fields above we arrive at 

r ..^J_^ f dk^i{P\b{{xP+M)bx{xP+M)\P), 0<x<l 

) Attx ^ J {2Tr)^ \{P\d{{xP+,k^)dx{xP+,k^)\P), -1< x < ' 

The above expressions show that fi{x) is determined from quark and antiquark density 
operators and hence has a probabilistic interpretation as the name PDF suggests. One can 
go further and expand the proton state in the light-cone Fock basis to generate contributions 
to the PDF in the form of an infinite sum of squares of light-cone Fock state wavefunctions. 

We can attempt the same procedure for the higher- twist structure function /^{x) defined 
implicitly through Eq. H1.18() 



h{x) = -J -e^-(P|V;(0)7-V'(An)|P). (C.9) 

The quark bilinear operator, however, only has bad components of the spinor field 

^{0)-f-^{Xn) = V2ipl (O)V'- (An). (C.IO) 

These components are dynamically determined from the equations of motion for the '0+ 
field, namely 

i|^V-(An) = i7+(-zp^)V^+(An). (C.ll) 

Clearly there is no partonic density interpretation for the f4,{x) structure function. The 
partial derivative term is kinematically of higher twist. The gluonic part of the covariant 
derivative leads to dynamical effects. Intermediate state gluons generally lead to more 
propagators compared to the kinematic terms and usually more suppression. 

A final point with respect to intermediate gluons must be made. One can imagine 
including longitudinal gluon scattering into the handbag diagram depicted in Figure 11.31 
Inclusion of these extra propagators does not necessarily lead to power-law suppression since 
large momentum flow can still be accommodated exclusively in the handle of the bag. These 
multiple gluon exchanges can thus be of the same order as the leading twist contribution. 
Summed to all orders these longitudinal gluon scatterings produce the gauge link 

V^(0)V'(An) ^ ^(0)P (^e-'3lx''-^dx''^ V'(An) (C.12) 
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Figure C.l: Diagrammatic depiction of the kinematic regions contributing to the GPDs. 
Figure adapted from [12]. 



which is necessary for the combination of bilocal fields to be gauge invariant in non-Abelian 
field theory. Above we have tacitly assumed light-cone gauge, n-A = 0, in which longitudinal 
gluon scattering is a gauge artifact. There is more subtlety in this choice than we have 
presented here [120,121]. 

In summary we have seen the leading- twist PDF fi{x) describes the density of partons 
in the infinite momentum frame and in light-cone gauge. One can go further and ask for 
the decomposition of GPDs in terms of the infinite momentum frame partons in light-cone 
gauge. Starting from the matrix element definition 

Fix,^,t) = \j ^e*^-(P'|V^(-An/2)7+V'(An/2)|P), (C.13) 

we again insert the mode expansion Eq. (|C.5() for the good fields and arrive at 

1 ^ r dk^ 

{p'\b\[ix - + A^Mix + op^,k^]\p), e < X < 1 

{P'\dx[i-x + OP^,-k^-A^]b-x[{x + OP^,^^]\P), -C<x<C • (C.14) 



X 



{P'\dl[{-x-OP M + A^]dx[{-x + OP\'i^^]\P), -1<^< 



-X 



Contributions to F{x,(^,t) are shown as partonic cartoons in Figure I(TT1 For x > one has 
the amplitude for removing a quark of a certain momentum fraction from the initial state 
and then inserting it with altered momentum into the final state. The same interpretation 
is true for anitquarks in the region x < — ^. When — < x < one has the amplitude 
for removing a quark-antiquark pair from the initial state. Notice as ^ — > 0, one recovers 
a density of partons with momentum fraction x. Similar to the case of PDFs, one can 
expand the proton states in the light-cone Fock basis. Doing so, one derives the overlap 
representation of GPDs [26,27]. 
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Appendix D 



Wavefunctions and Form Factors in 
(3 + 1) Dimensions 



In this Appendix we collect results relevant for Cliapter|51for wavefunctions and form factors 
in the (3 + l)-dimensional ladder model Eq. (|2.33|) . 

D . 1 Wavefunct ions 

Using the Bethe-Salpeter equation H2.1() and the definition of the hght-cone wavefunction 
\ipji) Eq. (|2.31j) we have the light-cone bound-state equation 



\i>R)=g{R)w{R)\i>R), 
where w{R) is the reduced auxihary kernel 

w{R) = g-^{R) G{R)W{R)G{R) g~^{R), 

with W{R) defined in Eq. (ITn|) . 

To leading order in G — G, one calculates w{R) for the ladder model to be: 

V{x,k^;y,p^\M^) = -{ xR+,k^\ w{R) \ yi?+,p^) 

- \e{x-y)D{x,k^;y,p^\M^) 



(D.l) 



(D.2) 



x-yL 



|(;r,k^) ^ (y,p^)}]e[x{l - x)]e[y{l - y)], (D.3) 




-I- 




Figure D.l: Diagrammatic representation of the one-boson exchange potential V appearing 
inEq. (jDJl). 
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Figure D.2: Light-front time-ordered diagrams contributing to the second-order kernel. 
Notice that diagrams A-D can be summed into the square of the first-order interaction and 
are consequently removed by the G subtraction term. The remaining diagrams E and F are 
the first of infinitely many required to maintain covariance. 



where we have defined 



I 2 o 



x-y 



1 



X 



(D.4) 



and taken R-*- = 0. Graphically this one-boson exchange potential Eq. (|D.3|) is depicted in 
Figure Ell 

For reference, the bound-state equation ()D.H) appears as 



V'(x,k^) = Dwix,k^\M^) j 



dydp^ 



2(2^)32/(1-2/) 



V{x,k^;y,p^\M^)iP{y,p' 



(D.5) 



A few comments about these results are in order. Truncation of the series in G — G vio- 
lates Lorentz symmetry because G is intrinsically related to projection onto the light cone. 
Although the covariant one boson exchange generates at leading order the two light-cone 
time-ordered one boson exchanges depicted in Figure irm there are infinitely many terms 
in the covariant ladder which are missing. To make this clear, we write out the second-order 
kernel in the reduction scheme. 

At second order, the reduced auxiliary kernel reads 



u;(2)(i?) = g~\R) G{R)V{R) \g{R) - G{R)\ V{R)G{R) g~\R) 



(D.6) 



For the moment, we ignore the subtraction term, insert the effective resolution of unity 
and evaluate the expression using Cauchy's residue theorem. The horrendous details are 
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contained in [72]. Contributions to w^"^^ are depicted as light-cone time-ordered diagrams 
in Figure ID.2I Notice that the diagrams A-D can be summed into the square of the first- 
order interaction. In this sense they are reducible contributions and are subtracted off by 
the G{R) term in Eq. HD.6|) . The remaining two diagrams (E and F) are characterized by 
four quanta propagating at a given instant of time. These contributions are irreducible 
and represent the first of a string of infinitely many interactions omitted from the covariant 
ladder. To second order the bound state equation reads 

«.,k-) = Av(x.k-|M^) / ,p y) 2(2,F t(i- .) ^'■•y.v^^.Wm-.'i'l 

(D.7) 

here the interaction V{x, k-*-; y, p-*-; z, q-*") is given by 

V{x, k^;y, p^;z,q^)= 2{2TTfy{y - l)5{y - z)5{p^ - q^)F(x, k^; q^jM^) 

- (xi?+,k^| u'(2)(ii) I zi2+,q^), (D.8) 

and the explicit functional form of the matrix element above is given in [72] and is repre- 
sented by diagrams E and F in the Figure. For brevity we have omitted writing the 
dependence. 



D.2 Form factors 



To calculate form factors, we use the electromagnetic vertex F'^ constructed in [38] up to 
the first Born approximation (notice that the ladder model's gauged interaction = 0) 



where denotes the electromagnetic coupling to scalars. 

Now using Eq. (12:231) to first order in G-G, the matrix element = (*p/|r^(-A)|^'p) 
then appears 

-1 



X ( a^'(-A)d2 ^ + y(-A)G(-A)a'^(-A)d 



X 1 + (G(P) - G(P))y(P) G(P) 



\1P) 



with the leading-order result 

Jilo = M \G{P')d^{-A)d^'G{P) |7P;. 
The first-order terms are 

64 = M\G{P')d'^{-A)d,'[G{P)-G{P))v{P)G{P) 
54 = (7p,| G{P')V{P')(G{P')-G{P')Xd^'{-A)d:^^G{P) 



IP) 
ip) 



54 = (7P'| G{P'){V{-A)G{-A)d^^{-A)WG{P) 



\1P) 



(D.IO) 



(D.ll) 



(D.12) 
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p P' 

Figure D.3: The leading-order diagram for the electromagnetic form factor. 



As outlined in [38], Eqs. (|D.11|) and ()D.12|) can be evaluated by residues being careful 
to remove two-particle reducible contributions by utilizing (|D.5|) . Here we state the results 
of these calculations in (3 + 1) dimensions. We denote as the momentum transfer and 
define A"*" = — C-P"*" (see Figure ID. 3|) . The leading-order result appears 

Jto = -^^P^ j 2(ty3^|l''-l)x- ('" - or (x-,k-^)^(x,k^), (D.13) 

where x' = and k'^ = k-*- + (1 — x') A-*- denotes the momentum of the final state. Using 
= + P'YF{t), Eq. (lDT3l) reduces to the Drell-Yan formula [36,80] for C = 0. 
The first of the leading order corrections is the Born term 5J^ . For x > Q we have 

, , +2iP+ f e(x - C)dxd'k^dydp^(2x - 

oJJ — 



7 (^>C) (I67r3)2 J xx'y{l - y)y' 

xV*(y',p'^)i^(y',p'^;x',k'^|M2)^!fc^Z)(y,p^;x,k^|M2)V;(y,p^), (D.14) 

y -X 

where y' = and p'"*" = p"*" + (1 — y')A.-^. This contribution corresponds to diagram A 
in Figure IT^ On the other hand, for x < C, we have 

, _ +2iP+ f e{C-x)dxdk^dydp^{2x-C)/C 

u 



7 (^<C) (I67r3)2 y y(l-y)ylx"{l-x") 

X V.*(y',p'^)Z)w(:r",k"^|t)^^^^^^Z)(y,p^;x,k^|M2)V;(y,p^), (D.15) 

y -X 

where x" = x/C, and k"-*" = k"*" + x" denotes the photon's relative momentum. This 
expression is diagram D in Figure ITTl 

The next leading-order term is the initial-state iteration 5J^ . The only contribution is 
for X > C) namely 

, _ +2iP+ f e{x - C)dxdk^dydp^{2x - C) 

1/ 



(167r3)2 J xx'{l — x')y{l — y) 

y -X 



xV^*(x',k'^)Z)(y',p'^;x',k'^|M2)^!^i^^Z)(y,p^;x,k^|M2)V'(y,p^), (D.16) 
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which corresponds to diagram B of Figure IT^ 

Lastly there is the final-state iteration term SJ~I^- For x > C we have 

, _ +2iP+ r 6{x-C)dxdk^dy'dp'^{2x-C) 
^'^f ix>0 - (I67r3)2 J x(l - x)x>y'{l - y') 

X ^*(y',p'^)Z)(y',p'^;x',k'^|M2)^!^i^^Z)(y,p^;x,k^|M2)^(x,k^), (D.17) 

y -X 

where implicitly y = C + (l~C)y' p"*" = p'"*" — (1 — y') A"*". This corresponds to diagram 
C in Figm'eEni While for x < the expression 



, _ +2iP+ r e{C - x)dxdk^dy'dp'^{2x - C)/C 
^-^f {x<0 - (I67r3)2 J {l-x)x"{l-x")y'{l-y') 



X r{y',P^)Dw{x",k"^\t /^^^ ''^ D{y,p^;x,k^\M^Mx,k^), (D.18) 



y -X 



corresponds to diagram E of FigureO Eqs. (|D.13HD.18() are then the complete expressions 
for the form factor up to first order ignoring the self-energy loop. 
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Appendix E 



Light-Front Time-Ordered 
Perturbation Theory 

The results of Chapters [21 and Ol can similarly be achieved directly from "old-fashioned" 
time-ordered perturbation theory in a form which utilizes projecting onto the two-body 
subspace of the full Fock space. For a nice, complete discussion of this formalism for 
the light-cone ladder model, see [122]. In this Appendix, we show how to derive higher 
Fock space components in this formalism, thereby demonstrating the generation of higher 
components from the lowest sector we found indirectly for GPDs and form factors in section 

We write the light-cone Hamiltonian as a sum of a free piece and an interacting piece 
which carries an explicit power of the weak coupling g. In an obvious notation this is 

p-=p-+gPf. (E.l) 

The free term P~ is diagonal in the Fock state basis, while the interaction generally mixes 
components of different particle number (in the scalar model we consider above, the inter- 
action is completely off-diagonal since there are no instantaneous terms). Let us suppose 
that in the full Fock basis, we have an eigenstate of the Hamiltonian, i.e. 

[P- +gPr)\i^)=p-\i^), (E.2) 

where the eigenvalue is labeled by p~ . Since the coupling is presumed small, the mixing 
of Fock components with a large number of particles will be small. Thus one imagines our 
bound state will be dominated by the two-body Fock component. 

To make this observation formal, we define projection operators on the Fock space V 
and Q in the usual sense. The operator V projects out only the two-particle subspace of 
the full Fock space and hence Q projects out the compliment. Let us define the action of 
these operators on our eigenstate 

r\i^) = IV'2) (E.3) 

Q\^) = \^q)- (E.4) 

As is well known, combination of Eqs. (|E.3|) and ()E.4|) leads to the following equation 
for the two-body Fock component 

P;s\^2)=p-\^p2), (E.5) 
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where the effective two-body Hamiltonian is 



^cff = Pvv + ^off (E.6) 

= P-p-p + PvQ~ ^=~PqV' (E-'^) 

and we have defined the fohowing notation for any operators A and B, P^b ^ AP~B. 
The effective two-body interaction V^s defined in equation HE.6|) is dependent upon the 
energy eigenvalue p~ since we have suppressed the degrees of freedom of the Q subspace. 
The relation between the Q-space probability (i.e. the non-valence contribution) and the 
effective interaction appears 

(^qI^s) = --^{HVeM- (E.8) 

In a weak-binding limit, we can series expand the effective interaction in powers of the 
coupling and thereby re-derive the light-front potential. Given that every boson emitted 
must be absorbed in the two-quark sector, we can have only an even number of interactions 
and hence ^ 

Ves = rgPrQ L^f2(-^\J] QgPfV. (E.9) 

P~ - Po \P -Po J 

So, for example, at leading order we have all possible ways to propagate from the two-body 
sector and back with only two interactions in between. The diagrams in Figure corre- 
spond to the two possibilities distinguished by the action of _ ^ between interactions. 

At the next order, we have all possible ways to propagate from two bodies to two bodies 
with four interactions in between, etc. 

To generate higher Fock components from the two-body sector, we necessarily must look 
at the Q-space state arrived at from Eqs. HE.3|1 and HE.4|1 

IV's) = PQp\i^2). (E.io) 

To generate an n-body Fock component from this state, we merely act with an n-body 
projection operator which we shall denote Q„. Similar to the above, we expand in powers 
of the coupling to find 

IV^n) = Qn-r^f: (-^] Q9Pr\i^2). (E.ll) 
P n=0 \P ^° J 

For example, the leading-order three-body state is obtained by attaching a boson to a quark 
line in the only two possible ways (and adding the light-front energy denominator at the 
end). With these three-body states, we can consider all possible three-to-three overlaps that 
would contribute to the form factor. These are depicted in Figure 1?^ (with the exception 
of a quark self- interaction) . The four-body sector is richer since there are two-boson, two- 
quark states as well as four-quark states. The two-to-four overlaps required for GPDs must 
have four quarks. At leading order, we generate the diagrams encountered above in Figure 
12.71 Not surprisingly directly applying time-ordered perturbation theory from a light-front 
Hamiltonian agrees with our derivation above from covariant perturbation theory in the 
Bethe-Salpeter formalism projected onto the light cone. 
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Appendix F 

Application to Timelike Processes 



In this Appendix, we pursue timelike processes in light-front dynamics. These processes do 
not have a Fock space wavefunction representation. Below we study time-hke form factors 
in particular to demonstrate the versatility of the light-front reduction of the BSE. This 
will allow us to make the connection to the generalized distribution amplitudes (GDAs), 
which we explicitly calculate for the ladder model. Analogous to GPDs, GDAs encode the 
soft physics of two-meson production and can thus be thought of as crossed versions of the 
GPDs. The GDAs enter in convolutions for various two- meson production amplitudes, see 
e.g. [124-127] . These distribution functions as well as time-like form factors are a theoretical 
challenge for light-front dynamics, since there is no direct decomposition in terms of meson 
Fock components alone. Furthermore, we shall see the leading-order expressions are non- 
valence contributions (which necessarily excludes a description in terms of constituent quark 
models). 

The time-like form factor F{s) for our model meson is defined by (see Figure If!T|) 

i^p |o) = -i{p - p'YF{s), (F.i) 

where s = (p + p')"^ is the center of mass energy squared. Now define P'^ = p^ + p'^ and 
C = p~^ /P~^ . We can work out the kinematics of this reaction in a frame where P"*" = and 



P 




Figure F.I: The triangle diagram for the time-like, pion form factor. 
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see that 



2P+' 

-_ (i-C> 

P 2P+ ' 

V>l = s{l-OC-M\ (F.2) 

where M is the meson mass. 

Similar to GPDs, the GDA for our model has a definition in terms of a non-diagonal 
matrix element of bilocal field operators^ 

^(^,C,^) = / ^e-^^--(vI/p ^p.\ qiO)i d+qix-) |0). (F.3) 

As in Section 13.21 our concern is only with the twist-two contribution and thus we have 
taken the plus-component of the vector current in Eq. (|_b\3() . The definition of the GDA 
leads directly to a sum rule for the time-like form factor 

f dz 

J ^^H^^C^^)=ns), (F.4) 

and hence a means to calculate ^ from the integrand of the time-like form factor. Again 
this is because the effect of the non-local light-cone operator is to fix the plus momentum 
of the active quark. 

Taking the appropriate residues of the five-point function, we arrive at Fig. IF. II for the 
time-like form factor. Keeping only the leading-order piece of the electromagnetic vertex 
F'^, we have 

^-^^Yiz"M-z"p^\M^)G{k,p)d-\k-p) 

X G{P-k,p')-f*{z',k^ -{l-z')p^\M^), (F.5) 

where we have made use of some auxiliary definitions: z = / P"^ , z' = and z" = zjQ. 
Recall 7(x|M^) oc Q\x(\ — x)\. This translates to: < < 1 and < < 1, and hence we 
do not pick up a contribution at zeroth order in the coupling. 

To work at first order, we pick up three terms analogous to those in Appendix for 
the spacelike form factor. We denote these as bJ^ ^ 6Jp and SJ^,. The Born term for the 
three-point electromagnetic vertex 6Jj is quite simple. For the same reason as the zeroth- 
order result, the restriction of 7(x|M^) oc 9[x{l — x)] and momentum conservation force 
the contribution 5Ji^ to vanish. This leaves us to consider only diagrams that arise from 
iteration of the Bethe-Salpeter equation for either of the final-state mesons. 

Considering first the term 6Jp, we have the contribution to the GDA 

^p{z,Cs) = ^PH2z-l) /^^^ (047*(^',k^-(l-^VM-^(A^-p) 

xj*iy",ci^-y"p^)Giq,p)Viq,k)Gik,p)GiP-k,p'), (F.6) 

^And similar to our discussion of GPDs, perturbative QCD corrections lead to renormalization scale 
dependence of the GDAs. For a complete discussion of 0{as) corrections in the crossed channel, see, 
e.g. [12]. 
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a > b 



Figure F.2: Leading-order diagrams for the time-like pion form factor. 



where we have chosen to abbreviate y = q'^/P^ and hence the label y" = y/Q- We have 
customarily omitted the subtracted term containing G, which is zero because there is no 
leading-order term to subtract. Requiring wavefunction vertices mandates < < 1 and 
C<z<l. 

Thus after integration over the loop energy, <I>p produces one contribution to the GDA 

(Ibvr-^j^C J z[\ — z)z'y" [\ — y") z — y 

X D{zM;V,^^\s)r{zM - (1 - /)P^)r(2/",q^ - yV)- (F.7) 



As a contribution to the time- like form factor, we can interpret Eq. ()F.7|) as the time-ordered 
diagram h of Figure IF. 21 



At first order in the weak coupling, we have one term remaining to consider SJ'^,. Again 
omitting the superfluous subtraction of G, we have 

%.{zX,s) = tP^{2z-l) I ^^j*{z",k^ - z"p^\M')G{k,p)d-\p- k) 

X G{P - k,p')V{P -k,P- q)G{P - g,p')7*(y',q^ - (1 - y')v^\M^), (F.8) 

with y' = Y^l^. Both y' and z are restricted: < y' < 1 and < z < Ci and the final 
contribution to the GDA is thus 

9{C-z) f dk^dy'dci-^i2z-l) ^ , ye{y-z) 



(167r3)2( 7 {I- z)z"{l- z")y'{l-y') y-z 

X D{y, q^; z, ]^^\s)r (^", " /'p^)^*(y', + y'p^), (F.9) 

In this form we recognize this contribution as diagram a of Figure IF. 21 

Having found the leading non- vanishing contribution to the GDA namely $ = $p + , 
we observe that the higher Fock components derived in Section f3. 21 fas well as in Appendix 
IE|) do not fit naturally into Eqs. ()F.7|1 or ()F.9|) . One needs a Fock space expansion for the 
photon wavefunction in order to have an expression for the GDA in terms of various Fock 
component overlaps. With the expressions derived for the GDA, we can use Eq. ()F.4|) to 
obtain the time-like form factor. Notice the expressions derived for the GDA and hence 
also for the timelike form factor depend on pieces of the light-front interaction kernel. Thus 
there is no way to directly utilize quark model wavefunctions for these observables. 
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Appendix G 



Complex Conjugate Poles in 
Minkowski Space 

In this Appendix we show that ampUtudes cannot be directly calculated in Minkowski 
space when complex conjugate poles are present in propagators and vertices. We then 
investigate how to calculate amplitudes properly in Minkowski space. Firstly amplitudes 
dependent on the imaginary part of some set of diagrams can be calculated by using a 
straightforward generalization of the cutting rules. We apply this to calculate the quark 
distribution function from the handbag diagram in the Bjorken limit in Section fG. 21 Lastly 
we consider the analytic continuation of space-like amplitudes from Euclidean space. This 
is complicated by the presence of Wick poles and requires their residues to be appropriately 
added when amplitudes are calculated. The details of the GPD calculation appear at the 
end of this Appendix in Section IG.4I Resulting functions calculated in Minkowski space 
after the Wick rotation agree with those obtained from the model defined in Euclidean 
space. 

This work lays the foundation for calculating light-cone dominated amplitudes using 
meromorphic model propagators and vertices constrained by both lattice data^ and Ward- 
Takahashi identities. Distribution functions for such models could be calculated rigorously 
since the pole structure of the propagator is known and relevant integrals converge in the 
complex light-cone energy plane. As far as light-cone phenomenology is concerned, resulting 
expressions would be truly Poincare covariant (as opposed to diagonal with respect to the 
non-interacting operators) and would satisfy field-theoretic identities. Filling these two gaps 
is essential for adequate hadronic phenomenology for processes at large momentum transfer. 
Moreover, such models would help light-cone methods and Dyson-Schwinger studies reach 
complementary standing. 

G.l Model defined 

The model we take is cf)^ theory with electromagnetic interactions. Equivalently we can 
view this model as a bound state of two scalar particles with a trivial Bethe-Salpeter vertex 
r(fe, P) = 1, where the coupling constant is assumed to be absorbed into the overall normal- 
ization, see Section 11.2.21 We make a simple Ansatz for the non-perturbative propagator 

^One must proceed with caution: the lattice calculations employ Landau gauge, while we have tacitly 
used light-cone gauge above. 
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Figure G.l: Triangle diagram at zero momentum transfer used to calculate the quark 
distribution function by projecting onto the hght cone. 



consisting of a pair of complex conjugate poles 

•^(^^ = (A:2_a2 + 52)2 + 4^252' (^-1) 

where a? — h"^ > 0. Defining for ease tit? = a? — iP' and e = 2ab (which is taken to be positive 
without any loss of generality), we can write the propagator as 

S{k) = — 75 75 ; h — 75 '-^—^ ;— . (G.2) 

f^z _ _|_ fc^ — — le 
The hght-cone energy poles of the propagator are thus 

_ k^' + is 

- 2k+ " 2k+ ^^-^^ 
and k'^t, = {k~)* , where * denotes the complex conjugate. Although we use a scalar model, 
results straightforwardly extend to spin-i particles, e.g., since only the pole structure of 
Eq. HG.1|) is relevant. 

Now let us consider calculating this model's quark distribution by projecting onto the 
light cone. The quark distribution can be derived by fixing the plus-momentum of the active 
quark x = k^ /P'^, see Figure IG. 11 and taking the plus-component of the current. Thus up 
to overall normalization, we have the expression 

q{x) oc j (fk S{k+ - xP+) xS{k)S{k - P)S{k). (G.4) 

The k^ integral is then performed by residues. Choosing a frame in which P-*- = 0, the 
spectator propagator has light-cone energy poles 

2 ■ .2 



k + m is 

2{k+ - P+) ~ 2{k+ - P+) 

and k^, = {k^)*. 

Performing the k~ integral in Eq. HG.4|) . we arrive at the quark distribution 



(x) oc 27ri| - e{-x)[Res{k^*) +Res{k^,)] 

+ 0[x{l - x)][Res{k-.) + Res{k^)] + 0(x - 1) [Res(A;- ) + Res(A;^- )] \ (G.6) 
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Figure G.2: Handbag diagram for the forward Compton amplitude. Dashed hne denotes 
the cut which yields the quark distribution in the Bjorken limit. 



This distribution does not have proper support, i.e. it is non-vanishing outside the inter- 
val X £ [0, 1]. Moreover, the distribution is not real valued, whereas it should be positive 
definite. Thus the model based on the propagator in Eq. ()G.1|1 cannot be suitably formu- 
lated in Minkowski space. We will find below that Eq. makes sense as a Minkowski 
space propagator only after analytic continuation from Euclidean space for the amplitude 
in question. In Section [4.2.31 we demonstrated that the model is perfectly well defined in 
Euclidean space and thus the challenge of analytic continuation remains. 

G.2 Cutting rules 

In this Section, we show there is some hope in working with the model propag ator Eq. (fnH) 
in Minkowski space. We demonstrate that the quark distribution can be derived by a 
straightforward generalization of the cutting rules. 

Consider the forward Compton amplitude T^'^[q,P) depicted in Figure lG.21 In the 
Bjorken limit, the imaginary part of this diagram is related to the quark distribution, see 
Section 11.11 For simplicity, we can choose a frame in which q-*- = 0. The minus-plus 
component of the forward Compton amplitude in such frames reads 

iT-+{q, P) = — / d^k {2k- + q-)S{k) S{k - P) S{k) {2k+ + q+)S{k + q). (G.7) 

In the scalar particle case, the minus-plus component of the forward Compton amplitude 
can be used to define the quark distribution in a way analogous to the spin-i case. The 
relation is simply 9m(T '") oc q{x), see [30]. In Eq. ()C.7|) . we have adjusted the overall 
normalization so that equality between 9m(T ^) and q{x) holds. With standard perturba- 
tive propagators, we could calculate the imaginary part of T^'^{q,P) by using the cutting 
rules, whereby one replaces the cut propagators in Figure IG. 21 bv an on-shell prescription, 
namely 

S{k)e=o — > -2m b{k^ - m^). (G.8) 

In Eq. (|G.8|) . we have specified e = for the case of a free particle propagator. Since 
large momentum flows through the handle of the handbag, we may neglect the mass of 
the struck quark and use the standard cutting rule for Sijn -|- q). In the Bjorken limit, we 
define x = —q^/2P ■ q which remains finite as (f' ,P ■ q ^ oo and is kinematically bounded 
between zero and one. Further we orient our frame of reference so that q has a large minus 
component in this limit, and consequently g"*" = (f' /2q~ is finite. Hence we have the familiar 
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replacement 

S{k + q) -— 6{k+ -xP+). (G.9) 
2q 

To complete the cut, we must deal with the spectator particle's complex mass shells. 
We must worry about the propagator Eq. (jG-lj) only where the denominator is zero. Thus 
we are lead to the cutting rule for the propagator Eq. (|G.1I) 

S{k) — > -TTi \5{k^ -m^ + ie) + 5{k^ - - ie)] , (G.IO) 

which puts the intermediate state on its complex mass shells. Furthermore, the limit e — > 
produces the regular cutting rule Eq. HG.8|) . 

Using this cutting rule for the spectator particle along with Eq. ()G.9|) . we can deduce 
the quark distribution from 9m(r ^) in the Bjorken limit^ 

q{x) = — [ dt^k 6{k+ - xP+) \6{k- - k-) + 5{k- - k-)] ^^ML^ (G.ll) 

■K J L ° ° \ \ — X 

where the light-cone energy poles are given in Eq. ()G.5|) . Notice the resulting distribution 
is real and has proper support due to the kinematic constraint x € [0, 1] imposed by the 
Bjorken limit. Evaluation of the two trivial integrals leaves only the transverse momentum 
integration 

^{^) = - E / ^ D^{xM.e,e"\M^)D^{xM.e\e"\M^) (G.12) 

where we have defined the Weinberg propagator generalized for complex masses as 

I)w(x,k^,.,e'|M2)-i = - ^''l + + ^ + (G.13) 

x(l — x) X I — X 

Evaluation of the k-*- integral yields an analytic expression for q{x), which is identical 
to that obtained from the DD Eq. (|4.49)) . Notice Eq. ()G.11|) is equivalent to evaluating 
Eq. ()G.4|) at Res(/i;^)+Res(A;^) and hence the quark distribution is derived as if the complex 
conjugate spectator poles both lie in the upper-half complex plane! We will understand this 
better once we analytically continue from Euclidean space. 



G.3 Analytic continuation 

Above we have seen that modified cutting rules can be used to derive the correct quark 
distribution function in Minkowski space. In essence the result stems from putting the 
spectator particle on its complex mass shells. This must be justified by careful analytic 
continuation of Euclidean space amplitudes. Below we consider the Minkowski space calcu- 
lation of the generalized parton distribution which cannot be derived by cuts. This will force 
us to deal with the underlying Wick rotation necessary to define the model in Minkowski 
space. 

^This may seem obvious. But there are two inconsistencies which will ultimately be resolved in Section 
IG.3I First Eq. ljC7^ is not the Minkowski space amplitude. As already mentioned it requires analytic 
continuation. Second the rule Eq. IG.IOL although plausible, does not produce the imaginary part of 
Eq. 1G.7II . The result Eq. IIG.12L however, leads us to believe these two wrongs make proverbial right. 
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Figure G.3: Complex light-cone energy plane for the propagator. 



In Section I4.2.3[ the model double distributions were calculated in Euclidean space. 
Thus to calculate related amplitudes in Minkowski space, we must Wick rotate in the 
complex energy plane: — > ik^. The analytic continuation can be viewed alternately in 
the complex light-cone energy plane. For k^ = 0, the rotation is from the Qm.{k~) axis to 
the ?R.e{k~) axis. In the general case, such a correspondence can only be made precise by 
considering the Wick rotation in terms of k^ and then boosting to the infinite momentum 
frame. This requires tedious algebraic manipulations and a proliferation of energy poles 
and time-ordered diagrams. Indeed it is easier just to imagine the rotation simply and deal 
with the light-cone singularities. This is the approach we present. 

Before tackling the model's generalized parton distributions in Minkowski space, let us 
imagine a simpler fictitious example. Consider some well-defined Euclidean space amplitude 
having only poles at k~ and k~,, see Eq. ()G.3|) . To calculate the amplitude in Minkowski 
space, we naively integrate along the ^e{k~) axis. In general the correct path on which 
to integrate is one which nears the Ke(fc~) axis except for detours around energy poles in 
the first and third quadrants. Such a path is correct since it can be continuously deformed 
into the Euclidean path. The difference between the naive integration and the correct path 
is a sum of residues of the Wick poles. The energy poles of our fictitious amplitude are 
depicted in Figure IG. 31 Their location depends upon the sign of x = k~^ /P~^. Thus for this 
amplitude the correct continuation from Euclidean space is 

J dQm{k-) — > J d^e{k-) + 2TTi[e{-x) - e{x)] Res{k-). (G.14) 

Notice only k^* is a Wick pole; this is expected because we know the limit e — > can be 
analytically continued in the naive fashion. Closing the contour in the upper-half plane 
to perform the Minkowski space energy integral^ and evaluate our fictitious amplitude, we 
pick up 2TTi[9{—x) Res(A;~) + 9(x) Res(A:~, )]. The net result according to Eq. (|G.14|) is thus 

27ri6'(-3;) [ Res(A;- ) + Res(A:-. )] . 



^One must be careful of zero modes [21] for which k'^ = 0. In such cases, the pole lies on the contour at 
infinity and the integration cannot be performed by residues. Since our fictitious example is only schematic, 
we are neglecting the issue of zero modes and are hence excluding amplitudes of the form 

V / d*k 

4l J (A:2 - m2 + ie)" ' ^ ' 



which must be handled separately. The example 



d^k 



is devoid of zero-mode complications and more closely parallels the expressions encountered for GPDs. 
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Looking back at Figure IG.31 the net result after Wick rotation amounts to both poles 
lying in the same half-plane; or equivalently, we have effectively integrated in either the 
right- or left-half plane. The result that both poles lie in the same half-plane reflects quark 
confinement; the quark wavefunction must vanish as time becomes large [89]. 

The space-like amplitude^ for the generalized parton distribution can now be continued 
to Minkowski space and hence be evaluated by projecting onto the light cone. To do so, 

we refer to Figure E31 and insert the non-local light-cone operator (p{0)iD~^ (f){z~ ) in place 
of the local twist-two operators denoted by a cross in the Figure. Here the plus-component 
picks out the leading-twist contribution according to light-cone power counting. Thus in 
momentum space, we arrive at 

H{x, C, t) = '^^^f [ d^k 6{k+ - xP+) {2k+ + A+) S{k)S{k - P)S{k + A). (G.15) 

Above we have included the ^-dependent pre-factor to normalize the action of between 
non-diagonal states. The overall normalization is then the same as in Eq. (|4.44() . By 
writing Eq. ()0.15|) in Minkowski space, we must also keep in mind the Wick residues 
implicitly necessary so that Eq. HG.15|) is meaningful. In addition to the poles k^ , k^ 
[given in Eqs. (|G.3|) and (|G.5|) . respectively] and their complex conjugates, the integrand of 
Eq. HG.15|) also has the poles 

and k~, = (k^)*. Carrying out the light-cone energy integration in Eq. HG.15|) as well as 
adding relevant residues resulting from the Wick rotation produces (the subtle details of 
this calculation appear at the end of this Appendix in Section [G.4|) 



H{x,C,t) = -2m e{x)e{C-x) Res(A;-) -FRes(fc^,; 

+ 2m9{x-Oe{l-x) Res{k^) + Res{k-) , (G.17) 

where the residue is of the integrand in Eq. HG.15() . As a result of the effective relocation 
of poles to the same half-plane as their complex conjugates, the resulting GPD Eq. HG.17|) 
is real and vanishes outside x from zero to one. 

Using the Weinberg propagator Eq. ()G.13|) . the residues can be compactly written in 
terms of relative momenta. Defining the relative momentum of the final state as 

x' = k'^ = k^ + (1 - x')A^, (G.18) 

and the relative momentum of the photon as 

x" = ^, k"^=k^+x"A^, (G.19) 



^There are additional complications for time-like amplitudes and for amplitudes involving unstable bound 
states. In these cases Wick poles are present even when standard perturbative propagators are used. The 
analysis above must be more carefully considered in these cases where threshold effects are already inherent 
in the analytic continuation to, or from, Euclidean space. 
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the light-cone GPD can be expressed in the form 



(1 - c/2) H{x, c, t) = e{x)e{Q - x) Hi{x, c, t) + e{x - c)e{i - x) H2{x, c, t), (G.20) 

where we have made the abbreviations 



2vr ^ ^f^^^ J x"[l - x"){l - X) 

Dy^{x,k^,e,e"\ M'^)D^ {x' , k'^ , e' , e" \ ) 



H2{x,C,t) = {2x-0^ [ 

£,e',e"=± 



dk' 



x{l — x)x' 

(G.22) 



Firstly one can see analytically that the correct quark distribution results at zero momen- 
tum transfer, namely H2{x,^,Q) = q{x), where q{x) is given by Eq. ()G.12|) . As remarked 
in Section IG.21 this function is identical to that obtained from the double distribution 
Eq. ()4.49|) . Secondly, the resulting light-cone GPD Eq. HG.20|) agrees numerically with that 
found from the double distribution, via Eq. 1)4. 4U|) . which is plotted in Figure 1^31 Finally 
the electromagnetic form factor found from the sum rule 



Jo 



dxH{x,C,t) (G.23) 



agrees numerically with the result of Eq. ()4.41|1 (which is plotted in Figure 14.4(1 . The 
^-independence of Eq. ()G.23|) stems from Lorentz invariance which is present, however, 
not manifest in Eq. (|G.20() . Thus with the calculation of Eq. ()G.20() from analytically 
continuing to Minkowski space, space-like amplitudes now agree with those calculated from 
the Euclidean space double distribution. 



G.4 Calculation of the GPD 

Below we derive Eq. ()G.17|) for the GPD in Minkowski space. The details have been 
relegated here since there is some subtlety. In order to evaluate Eq. ()(t.15|) , which implicitly 
needs analytic continuation, we must shift the energy integration variable and define a 
prescription for dealing with vanishing real parts. Let us see how these difficulties arise. 

In considering the Wick rotation, one is usually only concerned with the single denomi- 
nator that results from combining propagators via Feynman parameters. For the moment, 
let us ignore the complication of complex conjugate pairs of poles. In this case, combining 
the denominators of Eq. (|4.44() using Feynman parameters results in [P — Do{x , y; t) + ie]~^ , 
where Do is given by Eq. ()4.45|) . Since the bound state is stable and t is space-like, Do{x, y; t) 
is always positive and hence there are no Wick poles. Moreover, we have shifted the variable 
k'^ to arrive at P. 

In analytically continuing the expression with uncombined propagators (and again no 
complex conjugate poles), we are confronted with a problem. The pole , for example, is 
shifted by the energy P~ . Thus the location of the singularity in the complex plane will be 
shifted parallel to the real axis depending upon the relative magnitude of the spectator's 
kinetic energy and . In the schematic example shown in Section fG. 31 the pole k~ does 
not have such a shift. Thus for the b and c poles, the location of the singularities depicted 
in Figure Rt.3I will shift along the real axis (depending on P~ and A~) and there will be 
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Figure G.4: Complex light-cone energy plane for the shifted poles Eq. HG.26|) for the gen- 
eralized parton distribution. 



threshold Wick poles [the threshold is defined when 5Re(A;^^ = 0)]. This is unphysical: we 
just demonstrated the Wick rotation can be done for the combined denominators without 
crossing any poles. To perform the same Wick rotation at the level of separate propagators, 
we must use the freedom to shift the energy variable as well as the stability of the bound 
state. 

On the light-cone, the bound state stability condition can be expressed as 

where k-'-jei = k-*- — xP-*" is the relative transverse momentum of the two constituents. 
Because we consider the elastic electromagnetic form factor, there is an analogous relation 
for the final state P' . Since each propagator contains the kinetic energy of a single particle, 
the bound state stability condition can never be utilized without shifting k~ . Yet in order 
to perform such a shift, the real part of one pole must be zero and hence we must invent a 
prescription for moving this pole off the Euclidean contour. 

We choose the translation: ^ + 6k^ + rj, where is a positive infinitesimal and 



k^^ + 



m 



2 



The resulting poles of the integrand in Eq. (jG.lSf) we denote 

k +m le 







-6k~ 


= -P- 


K 


= K 


-5k- 


-v = - 




= k- 


-5k~ 


= -P'~ 



+ 2P+x{l-x) 2P+X 

V - wf^ ' (G.26) 

I k'^ +m'^ iE 

2P'+x'{l-x') 2P+(x-0 

and similarly for their complex conjugate partners. Notice imaginary parts of the poles 
are unaffected by the energy translation. The rj prescription has displaced the resulting 
spectator pole away from the Euclidean path independent of x. Moreover the a and c poles 
have non-zero real parts; so r/ has been set to zero for these poles above. 
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Using the expressions for the new poles Eq. HG.26|) and the bound state stability con- 
dition Eq. (|G.24() . we can determine the quadrant location of the singularities independent 
of P~ and (see footnote 2). These quadrant locations are depicted for the full range of 
X in Figure IG.4I Accordingly poles in the first and third quadrants are Wick poles. As we 
saw in Section fG. 31 the net result of the analytic continuation is to evaluate the integral by 
effectively closing the contour in the right- or left-half plane. Hence the GPD Eq. HG.15|) is 



H{x, C,t) = -2Tri (9(x)(9(C - x) Res{k^ ) + Res(fc~. ; 

+ 2TTi e{x - C)0{1 - x) Res{k-) + Res(fc^-. 



(G.27) 



Evaluating the residues in Eq. (|G.27|) yields the result of Section IG.31 namely Eq. HG.20|) 
which is algebraically equivalent to Eq. HG.17() . 

Notice from Figure IG. 41 other rj prescriptions for the shift, such as + 2p+(L-i) ' ^^^^ 
an incorrect result for the case when there are no complex conjugate pairs. The figure shows 
that the infinitesimal prescription must be positive and independent of the sign of x, x — C, 
etc, in order to reproduce the familiar result. It is interesting to note that for x > 1 in the 
case without conjugate pairs, all poles of the integrand are Wick poles. Though since the 
integral is convergent, the sum of these Wick residues vanishes. 

The interested reader can verify that the alternate shifts which use the same pole pre- 
scription 



' ^ j^- -i- k^^+m^ 



+ 



k- 



k- 



2P+X 



+ r] 



2P+(x-C,) 

also yield the correct results provided t is space-like. 



+ -q 



(G.28) 
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Appendix H 

Ambiguity of the Double 
Distributions 



The decomposition into F{f],a;t), G{(3,a;t) and K{(3,a;t) DDs is ambiguous. Experi- 
mentally one cannot access the DDs directly, only the H{x,(,,t) and E[x,^,t) GPDs. In 
Ref. [75], the ambiguity inherent in defining DDs for the pion was likened to the gauge 
ambiguity of the vector potential of a two-dimensional magnetic field. In this Appendix, 
we extend this analysis to the proton case. A full treatment has been given in [119]. 

Utilizing the Gordon identities (which appear in Appendix^), we can rewrite the bilocal 
matrix element in Eq. as 



(P',A>,(-z/2) {z/2)\P,\) 



1 - 

AM- 



1 /-l-l/^l - 

1 J-i+m 



P-z 



M 



GE{(3,a-t) 



4M2 



G{P,a;t) 



uxiP), (H.l) 



where we have defined new double distributions 



GEiP, a; t) = F{P, a; t) + -^K{(3, a; t) 
Gm{I3, a; t) = F{p, a; t) + K{P, a; t) 



(H.2) 
(H.3) 



in analogy with the Sachs electric and magnetic form factors. After integrating by parts 
and neglecting the surface terms, we have 



(P', (-^/2) M (^/2) 1^' ^) = M y ^ dp 



da e-^/5^-+-^-/2 



i+L 



X -^—^nyiP'n 



AM 



N{p, a; t) - ie^-°^z^A,P„7^75 Gm{P, a; t) 



ux{P), (H.4) 



where 



NiP,a;t) = ^^GM»;t) + l[l ^ 



t \ d ^.^ , 
— G(/3,a;t). 



(H.5) 
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Now we define an arbitrary potential function xiP^ ol\ ^) which has the property: 
x(/3, = — — The above expression Eq. (|H.4() is invariant under the trans- 
formation 



GEi/3,a;t) -^GEiP,a;t) - 1 



d 



4M2 / da 



d 



G{(3, a; t) ^ G{(5, a; t) + 2— xW, a; t) 



(H.6) 



which imphcitly leaves GmW, ol\ t) unchanged. Notice the a symmetry of the distributions is 
preserved because the potential «; is an odd function. Translating the transformation 
Eq. (|H.6|) back into the original DDs, we have invariance under 

a; t) ^ a-t)-^ x(/5, a; 

G(/3, a; t) ^ a; t) + 2^ x(/3, «; (H.7) 

One may try and choose the potential so that a; t) can be completely removed. 
The analysis is now the same as in [75]. To choose such a x(/?)Q?;i), we must maintain the 
zero boundary conditions so that surface terms produced by partial integration vanish. The 
following choice of potential preserves the boundary conditions 



Xo(/3,a;t) 



/3 



where the D-term is given by 



dp' G{p',a;t) 



l-\a\ 



dp' G{p', a; t) - sign(/?)D(a; t) 



(H. 



D{a;t) 



dpG{P,a;t). 



(H.9) 



Under the transformation generated by Xo{P,0!;t), the function G{P,a;t) is changed ac- 
cordingly to 

fo, if G(-/3,a;t) = G(/3,a;t), 
1 6{P)D{a; t), otherwise. 



G{P,a;t) 



(H.IO) 
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Appendix I 

Basic Identities 



For reference we list identities used in computing DDs above. To calculate numerator 
structures so that symmetric and traceless tensors can be calculated, we used generalized 
Gordon identities [130] of the form 

uyiP')Tux{P) = 4^^A'(^') (2mr} + [A,r]) UxiP), (i.i) 

where T is any Dirac matrix, = + Py and A'^ = {P' - Py. The usual Gordon 

identity is a special case of Eq. (|I.1|) . namely for F = 

ux'{P')rux{P) = ^^A'(^') (2P^ + i'T^'^u) uxiP). (F2) 
We also require the following two cases of the general identity with 75 = i7''7^7^7^ and 

ux'{P'h^l5 uxiP) = i^ux'{P') (a^5 - e^^'^^cj^pP,) ux{P) (F3) 
ux'{P')cJ^''ux{P) = ^^x'iP') (t'^A^ - 7'^A- + 2ie^'^"^P,7^75) ux{P). (F4) 
Also of use is the gamma triple product 

^M^-^P = S^^-P-^^ + ^£'^-'"^7,75, (F5) 

where 
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Appendix J 

Calculations on the Light-Cone 



Here we include our conventions for projecting quantities onto the light cone. While the 
development and derivations above rely exclusively on manifest Lorentz invariance, the 
light-cone Fock representation provides transparent physical intuition about our model. 
The light-cone wavefunctions for our model are admittedly simple and thus provide a useful 
guide to understanding the DDs constructed above. 

The light-cone spinor u\{k,m) satisfies the Dirac equation — m)u\{k,m) = and is 
explicitly given by (e.g., see [129]) 

ux{k,m) = -^^-j={y2k^ +pm + cx^ -li^^Xx, (J.l) 

where /3 = and ex = 7^7 and the unit spinors Xx are given by 

x| = i=(i,o,i,o) 

Using Eq. (|J.1|) . we derive the useful product of spinors of different mass and momentum 



ux'{Km)ux{P,M) = [6x',x{k+M + P+m) - X6x',-x{k-^Px - P'-kx)] , (J.2) 

\/k+P+ 

where the notation ax = + iXa? has been employed with the understanding that spins 
correspond to the signs t= -|-1 and |= — 1. 

Using the light-cone spinors, one can find the Hq{x,^,t) and Eq{x,^^t) GPDs in terms 
of the light-cone, non-diagonal matrix element '^(x,^, t), namely 

A^^'^(x, i, t) = J— [(1 - e)Hq{x, t) - fEgix, t)] (J.3) 

yi - C 

M^q'-Hx,^,t) = - J't^ Eq{x,^,t). (J.4) 
21V1 Y 1 — 

These expressions can alternately be used to derive GPDs for the model considered in Sec- 
tion l5.2l bv directly inserting the quark bilocal operator between non-diagonal proton states. 
Expressions are derived in this manner by integration over the relative light-cone energy 
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at the cost of sacrificing manifest Lorentz invariance. This description in terms of hght- 
cone Fock components, however, is more intuitive to work with than the DD formulation. 
Moreover, the light-cone energy integration obviates the positivity properties of our model 
GPDs. We shall not present complete expressions for the GPDs on the light-cone, however, 
the diagonal overlap region 2; > ^ is particularly simple to consider. 

For our model, we can find the lowest Fock component's light-cone wavefunction by 
projecting the covariant amplitude onto the surface x"*" = 0, namely 



, / ,± u 1 f dk Us^{P - k,m) f ^ \_ T 
VLc{x,KehSi,X) = ^2P+)^ J ^/T^ yY-i->iC'jUs.,[P -k,m) 

^-j+^{k,P)ux{P,M), (J.5) 



Us^{k,m) 



where x is the fraction of the proton's plus momentum carried by the residual quark [x = 
k^ /P^), and the relative transverse momentum is k^j = k"*" — xP"*". Above, A labels the 
spin of the proton, whereas the Sj label spins of the three quarks. We have omitted the color 
and isospin parts of the wavefunction which are trivial: tx eci,c2,c3(<5i,u'^2,d — ^i,d^2,u)S3,u for 
the scalar diquark. In the scalar channel, the quark-diquark Bethe-Salpeter wavefunction 
is 

and above we have chosen a single Dirac component F^*) = —ig^^^ and the coupling constant 
g has been absorbed into the overall normalization above. Our choice of vertex functions 
corresponds to only modeling a subset of the possible three quark light-cone wavefunctions 
of the proton [132]. Carrying out the projection in Eq. ((XsJ, yields 

lpLC[X,k^^l] Si, X) OC — -TyT- (J-7) 

* X{l — X) X 1 — x 

Above a = 1 for the valence light-cone wavefunction. The effective wavefunction 
ipf^ix^k.-^^^] Si, X) generated by the analytic regularization is given by Eq. (|J.7|) with a = 2 
(or in general a > 1). There are no true higher Fock components in this model since the 
interaction kernel is instantaneous in light-cone time [47]. Focusing on the spin structure, 
the scalar diquark channel consists of two states: a state where the proton spin is aligned 
with residual quark's spin and a state where the quark and diquark are in a relative p-wave. 
The quark distribution function can be obtained from 

fuix) = Y, J dk^\^Pll{x,k^;s„X)\\ (J.8) 

and agrees with the covariant calculation of Quix) from its moments in Section [5.2l Similarly 
the Dirac and Pauli form factors can be expressed in terms of the effective wavefunction 
since the light-cone contour integration only receives contribution from the diquark pole. 
The expression are 

Flit) = J^y"(ik^(ix^£f(x,k^ + (l-x)A^;s„AM^(x,k^;s„A) (J.9) 
-^F2{t) = ^y"dk^dxV'rc(x,k^ + (l-x)A^;s„-A)^£ff(x,k^;s„A)(J.10) 
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Also for this reason, we can express the GPDs as simple convolutions in the diagonal overlap 
region x > Let xi = and X2 = yE| a-nd then we have 

e{x-i)Ml^\x,i,t) = Y, I tZk^^LC (^2,k^ + (l-X2)^;s„A) 

xV£c(^xi,k^-(l-xi)^;s„A^ (J.ll) 

with a very similar expression holding for the spin-flip amplitude. In the above form the 
positivity bound is manifest. 

In the axial-vector diquark channel we have the orthogonal amplitude 

i± ^^ 1 f dk- Us^{P - k,m) ( + \_ t 
V'(3;,k,,i;si,A) = ^2P+Y J ^JY^ \^^lp,iCy^Us^[P -k,m) 

Above, we have omitted the color and isospin parts of the wavefunction which are: oc 
£ci,c2,c3[iSi,uS2,d + Si^dS2,u)S3,u - 25i,u(52,M(53,d] for the axial diquark. In the axial-vector 
channel, the quark-diquark Bethe-Salpeter vector wavefunction is 

= J '—-(- i9^''hu75)Df^''iP - k), (J.13) 

where we have chosen a single Dirac component F^") = —ig^"-^ and the coupling constant 
above. Carrying out the projection in Eq. (|.T.12|) yields the axial diquark 
piece of the light-cone wavefunction. Analogous formulas to the above hold for the axial 
diquark contribution to the quark distribution, form factors and GPDs for x > ^. 
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Appendix K 

Derivation of the Double 
Distributions 



In this Appendix, we detail the calculation of the DDs in the scalar diquark channel and 
comment on the calculation in the axial-vector channel. The crucial steps in the derivation 
hinge upon reducing the numerator by factors present in the denominator or by use of the 
simple identity: 

-^^i'^-^f'^'i = 0. (K.l) 

In the scalar diquark channel, let us take Eq. H5.31() as our starting point. Denote the 
propagators simply by 2t = (A; — P + A/2)^ — m|) + ie, *B = (fc + A)^ — w? + ie, and €. = 
k"^ — m? +ie. The DDs in the scalar channel can be deduced without reducing factors in the 
numerator. We merely introduce two Feynman parameters {x, y} to cast the denominator 
specifically in the form [x'^ + y!B + {1 — x — One then translates A;^ to render the 

integral hyper spherically symmetric via the definition k'^ = + f3P^ — (a + l)A^/2. Here 
P = x and a = X + 2y — 1. Carrying out this procedure on Eq. 1)5. 31() produces 



2 70 7-1+/3 J 4-^T{ay 



where the numerator Dirac structure is given by 

TV^ = (f+/57^-(a-l)A/2)7^f+/3F-(« + l)A/2) 

+mV + ima^'^A^ + 2m{l + (5P - aA/2)^. (K.3) 

Using Eq. (|K.H) . as well as the identities in Appendix ^ we can cast Eq. ()K.2|) in the 
form 



/ da d^^ ,2 Hl-Do{P,a;t)] 

^ Jo J~l+f3 J ^ '- \'^) 



x,J.(P"')A^W.(P)g^j^--^P«...P--.(--) ^^^(--) , (K.4) 
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where 




+icj^^A^(l -P){m + I3M) - a{m + /3M)A''. 



(K.5) 



The / integration is then standard and the DDs can be read off simply by using Eqs. (|5.1|) 
- H5.4|) and we arrive at Eq. H5.35() . 

Calculation in the axial-vector diquark channel is similar, however, the numerator is 
more complicated. In the same units as Eq. (j5.I^H) . we have the contribution from Fig. 15.11 
for the axial diquark 



where the vector propag ator Z?"^(/c - P) is given in Eq. (|5.25|) . The terms which result 
from the g"^ structure in the vector propagator can be dealt with straightforwardly after 
evaluating the contracted gamma matrices. The integrals encountered are then similar to 
those in the scalar channel. 



9 



nA'(^^^)757a(^ + A + m)r^^i-/^"(^ + m)7^75nA(P) 
-m? + ieY[{k + A)2 - m? + ieY 



(K.6) 
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